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SECTION  I 


INTRODUCTION 


Understanding  the  performance  of  high-powcr  lasers  requires 
a  model  of  the  system  response  to  external  perturbations.  In  the 
case  of  small  perturbations  about  some  equilibriam  value  the  sys¬ 
tem  response  can  be  characterized  by  a  (frequency-dependent)  trans¬ 
fer  function. 

As  the  simplest  example  of  such  a  transfer  function,  consider 
the  problem  of  relaxation  oscillations  in  a  conventional  laser.  The 
system  can  be  described  by  the  following  set  of  equations: 

It  =  c  (g  -  (1) 


_  2  (1  +  ^)  +  p  (2) 

dt  T  'tig 


where 


())  =  system  flux 
g  =  system  gain/unit  length 
d)  =  saturation  flux 
T  =  upper  state  deactivation  time 
c  =  speed  of  light 
g^  =  cavity  coupling 
P  =  external  pumping 

Let  us  now  assume  steady-state  pumping  with  a  small  per¬ 
turbation  ,  i .  e .  , 

P  =  P  +  P, ,  (3) 

o  1 
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where  is  the  steady-state  component  and  is  the  perturbation. 
Similarly,  if 


<{>  =  4>o  +  'J’l 


(4) 


g 


g  + 


(5) 


r  =  1  +  (}>  /<}> 
o  ^s 


(6) 


we  have  the  following  equations  (after  Laplace  transforming,  de¬ 
noted  by  4)  (s)  ,  etc  .  ) 


In  this  case 


If  any  singularities  occur  with  Rea]  s  >  O,  tlic  system  has 
an  instability. 
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In  this  case,  the  roots  occur  at 


s 


3go(r-i: 


(10) 


i.e.,  at  a  frequency 


w 


o 


(r-1) 

T 


(11) 


with  a  damping  factor  of 


r  =  2r/T 


(12) 


and  no  instability  obtains. 

In  this  report  we  will  show  that  a  realistic  transfer  func¬ 
tion  for  high-power  lasers  must  include  the  acoustic,  or  thermo¬ 
dynamic  degrees  of  freedom  as  well.  The  reason  for  this  inclusion 
car.  be  seen  in  schematic  form  in  Figure  1. 

Gain  and  flux  perturbations  load  to  density  perturbations 
via  local  heating.  This,  in  turn,  causes  refractive  index  fluc¬ 
tuations  and  flux  perturbations  which  close  the  cycle. 

We  v/ill  show  below  three  separate  examples  of  this  e.xpandod 
approach  to  the  system  transfer  function  and  system  stability.  The 
first  two  deal  with  CO2  EDLs,  although  the  form.alism  can  bo  gen¬ 
eralized  to  other  types  of  lasers.  The  third  deals  explicitly  v.’ith 
the  chemical  laser. 

We  shall  see  that  when  the  tliermodynamic  degrees  of  freedom 
are  included,  the  possibility  of  a  system  instability  exists.  In 
fact,  such  an  instability  has  actually  been  found  experimentally 
in  the  case  of  a  cw  CO2  electron-discharge  laser  (EDL) . 

VJe  shall  now  briefly  outline  the  remainder  of  tliis  report. 
Section  II  contains  a  discussion  of  an  instability  arising  from 
the  transfer  of  energy  from  an  initially  oscillating  cavity  mode 
to  a  different  transverse  mode  by  an  acoustic  'wave.  Tl".c  approach 
to  this  problem  is  the  solution  of  a  picrturbed  propagator  equa¬ 
tion,  i.e.,  an  integral  equation  for  the  change  in  the  electro¬ 
magnetic  wave  after  a  round-trip  traversal  of  the  acoustically 
perturbed  cavity.  Tv;o  simplifying  assumptions  made  in  this 
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Figure  1.  Schematic  Indication  of  the  Coupling  of  Flux,  Gain, 
and  Thermodynamic  Degrees  of  Freedom 
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section  arc  that  (a)  the  initially  oscillating  mode  has  a  uni¬ 
form  transverse  flucncc  distribution  and  (b)  that  both  electro¬ 
magnetic  modes  have  the  same  natural  frequency  of  oscillation. 

It  is  shown  that  the  oscillatory  heating  v/hich  occurs  can  result 
in  an  unstable  interaction  between  the  two  electromagnetic  waves 
and  the  acoustic  wave,  i.e.,  a  perturbation  whose  amplitude  in¬ 
creases  exponentially. 

In  Section  III  we  discuss  a  different  type  of  instability: 
alignment  sensitivity.  Variations  in  the  gas  density  along  the 
flow  direction  lead  to  both  gain  variations  and  slight  changes  in 
the  propagation  direction  of  a  light  mode.  We  show  that  this  can 
result  in  extreme  cavity  alignment  sensitivity. 

In  Seacion  IV  we  return  to  the  mode-medium  interaction 
problem  of  Section  II,  which  is  now  reformulated  in  the  language 
of  Brillouin  scattering  (differential  treatment).  We  show  the 
equivalence  of  the  differential  and  integral  formulations  of  the 
mode-medium  interaction.  We  focus,  however,  on  an  alternative 
instability  mechanism.  We  show  that  an  acoustic  wave  of  fre¬ 
quency  tends  to  resonantly  couple  an  initially  oscillating 
mode  to  a  new  transverse  mode  whose  natural  frequency  is  off¬ 
set  by  from  the  natural  frequency  of  the  intial  mode.  This 
is  just  the  ordinary  Brillouin  effect.  Wo  derive  equations  of 
motion  for  the  amplitude  of  the  now  mode  and  for  the  acoustic 
wave.  The  secular  equation  for  the  system  is  derived.  It  is 
shown  to  be  unstable  in  the  limit  of  negligible  gain-flux  oscil¬ 
lations  (that  is,  with  a  ratio  of  cavity  flux  to  saturation  flux 
<<  1,  so  that  gain  saturation  may  be  ignored).  A  brief  discus¬ 
sion  is  given  of  the  nonlinear  saturation  of  this  process. 

An  idealization  of  the  mode-medium  interaction  theory  of 
Sections  II  and  IV  is  the  implicit  assumption  that  the  not  loss 
(=  loss  minus  gain)  of  two  cavity  modes  may  be  very  nearly  equal. 
(See  Appendix  A.)  At  the  beginning  of  this  program  it  v/as  realized 
that  no  self-consistent  calculation  of  the  eigcnmodcs  of  a  loaded 
(i.e.,  gain-saturated),  unstable  resonator  of  moderate  Fresnel 
number,  such  as  a  cw  electron-discharge  laser,  existed.  Such  a 
calculation,  which  clearly  is  quite  relevant  to  this  program,  is 
described  for  the  first  time  in  Section  V.  Using  a  new  computa¬ 
tional  technique  developed  at  AERL,  we  have  solved  for  the  self- 
consistent  eigcnmodcs  of  an  unstable  resonator  whose  gain  is 
saturated  by  the  presence  of  a  strong  oscillating  mode  (the 
"fundamental"  mode).  We  find  that  the  fluence  distributions  of 
the  various  eigcnmodcs  arc  modified  by  the  saturation,  and  that 
the  eigenvalues  are  shifted  from  their  empty- resonator  values. 

In  particular,  the  empty-resonator  loss  degeneracy  which  exists 
at  certain  hal f- in tegra 1  values  of  Fresnel  number  disappears  in 
the  presence  of  a  saturated  gain,  i.e.,  different  modes  are  found 
to  have  different  losses.  This  is  a  significant  new  result  which 
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suggests  that  it  is  not  advantageous  to  avoid  ha  1 1- integra 1  values 
of  Fresnel  numbers  in  the  design  of  an  unstable  resonator  which 
will  be  operated  in  tlie  siiturated-gain  regime.  Fx'om  the  point  of 
view  of  our  program,  the  results  of  this  section  suggest  that  th.: 
model  of  Sections  II  and  IV,  which  assumes  that  at  least  some  modes 
have  nearly  equal  net  loss,  is  somewhat  simplistic.  However,  the 
nonuniformity  of  the  gas  density  has  been  neglected  in  this  sec¬ 
tion;  it  is  possible  that  when  the  nonuniformity  is  taken  into 
account  the  net  losses  of  several  modes  may  be  very  nearly  equal. 

In  Section  VI  wo  extend  our  instabi 1 ity- theory  techniques 
to  a  relatively  simple  model  of  a  supersonic  chemical  laser 
(HF/DF).  Tile  model  corresponds  to  a  one-dimonsicna  1  flowfield. 

The  fluid-mechanic,  chemical  and  optical  coupling  equations  are 
obtained  and  Laplace  transformed.  The  stability  of  the  system 
is  studied  by  use  of  the  Routh  Criterion,  which  allows  the  sign 
of  the  roots  of  a  polynomial  equation  to  be  determined.  Numerical 
calculations  for  typical  operating  conditions  of  the  model  laser 
show  that  it  is  stable.  The  important  result  of  this  model  is 
that  it  determines  the  transfer  function  of  the  system.  The 
transfer  function  can  be  used  to  predict  the  magnitude  of  the 
laser  response  to  any  external  perturbation. 

Finally,  in  Appendix  A,  v;c  discuss  the  acoustically  per¬ 
turbed  optical  resonator  in  a  more  rigorous  fashion.  Using  the 
integral-operator  formalism,  wo  derive  tl’ic  equation  satisfied 
by  the  amplitudes  of  the  various  eigenmodes  in  the  presence  of 
the  acoustic  perturbation.  Wo  reformulate,  in  the  language  of 
quantum  mechanical  perturbation  theory,  the  difference  between 
the  instability  mechanisms  discussed  in  Sections  II  and  IV;  both 
follow  from  the  present  treatment.  We  suggest  an  alternative 
model  to  explain  the  experimentally  observed  flux  modulation  in 
the  electron-discharge  laser:  one  in  which  only  a  single  eigen¬ 
mode  oscillates,  but  an  acoustic  standing  v/avo  both  modulates  the 
cavity  loss  and  is  driven  by  nonuniform  heating  caused  by  the 
nonuniformity  of  the  eigenmode  and  changes  of  the  cavity  flux. 
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SECTION  II 


A  THREE-WAVE  LASER  INSTABILITY 


A.  INTRODUCTION 

It  has  been  known  for  some  time^^^  that  high-power  CO_  cv/ 
Electric  Discharge  Lasers  (EDL)  exhibit  output  instabilities. 

This  is  evidenced  via  a  spontaneous  oscillation  of  the  output 
flux  usually  resulting  in  full  modulation. 

An  example  of  such  behavior  for  a  typical  CO2  cw  EDL  is 
shown  in  Figure  2.  To  the  left  of  point  "a"  the  output  is  stable, 
while  to  the  right  of  point  "b"  the  instability  has  already  built 
up.  Between  "a"  and  "b"  the  output  may  be  considered  ] inoar ly 
unstable ,  and  periodic  in  nature.  This  suggests  the  possible 
success  of  a  linearized  model,  with  the  (linear)  instability 
evidencing  itself  via  an  exponentially  grov;ing  root  of  the  char¬ 
acteristic  equation  of  the  system. 

Corroboration  of  such  a  model  may  be  obtained  by  reviewing 
interf erograms  of  the  laser  taken  with  flux  extraction  (Figure  3). 
Also  shown  is  the  laser  output  flux  versus  time  with  the  points 
corresponding  to  the  three  interf  erograms  as  indicated.  \<c  note 
the  presence  of  low  amplitude  density  (or  acoustic)  standing  waves 
between  anode  and  cathode  whore  the  instability  in  the  output  flux 
is  already  evident.  It  is  as  if  the  laser  were  being  acoustically 
Q-switched. 

Further  evidence  of  an  acoustic  connection  can  bo  found  in 
the  period  of  the  output  oscillations.  Typical  oscillations  of  a 
gain-flux  system,  such  as  are  found  in  relaxation  of  a  gain- 
switched  spike,  are  of  the  order  of  a  few  microst'conds .  As  can 
be  seen  from  Figures  2  and  3,  the  period  of  the  oscillations  is 
of  the  order  of  50  gsec.  For  typical  anodc-cathexio  separation 
of  about  5  cm,  with  an  acoustic  mode  number  ~  5  (Cf .  Figure  3) , 
the  acoustic  time  is  typically 


(5  cm/5) 

4 

3  X  10  cm/scc 


~  35  ijsec 
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Figure  2.  Evidence  of  Linear  Instability.  The  upper  trace  in  the 
top  figure  is  the  flux  and  the  bottom  is  sustainer 
current.  The  first  three  msec  show  a  dc  output  where¬ 
upon  the  flux  begins  to  oscill.itc.  At  this  loint  there 
is  a  floating  zero  level  duo  to  detector  sensitivity. 
The  lower  figures  are  expanded  time  scales  at  incipient 
breakup  of  flux. 
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INTERFER06RAMS  OF  CW  EOL  CAVITY  WITH  FLUX  EXTRACTION,  20.000fps 
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OUTPUT  FLUX  TRACES  FOR  MOVIE  INTERFEROGRAmS 
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Figure  3.  Mach  Zonder  Interf erograms  of  the  Flov;.  The  output  flux 
trace  is  shown  as  the  lower  curve  in  the  upper  right  with 
positive  downward.  The  flux  is  stable  up  to  about  2  msec 
whereupon  the  Q-switching  occurs.  The  three  interf ero¬ 
grams  are  before,  at,  and  after  onset  of  full  scale 
oscillators.  Each  picture  shows  two  exposures  separated 
by  25  psec.  The  large  scale  waves  evieJent  at  4.8  msec 
are  a  result  of  the  oscillations  and  not  the  caur.e  as 
they  do  not  appear  earlier.  At  2  msec  notice  the 
presence  of  a  standing  transverse  wave  of  a  few  friiini':'. 
and  about  3  or  4  cycles.  It  is  this  v/ave  ttiat  proiiuces 
the  oscillations  in  the  flux. 
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In  this  section,  v/o  will  describe  a  model  of  flux-gain- 
acoustic  (density)  coupling  which  admits  an  instability  and 
derive  an  analytic  expression  for  the  unstable  root  in  certain 
limiting  cases.  A  subsequent  paper  will  present  a  more  rigorous 
formulation,  comparison  with  experiment,  and  experimental  methods 
of  removing  the  instability. 

B.  THE  MODEL 

A  schematic  model  of  the  CO2  CV-JEDL  is  shown  in  Figure  4. 

The  flow,  anode-cathode,  and  optical  axis  directions  are  along 
the  X,  y,  z  axes,  as  shown  in  Figure  4. 

The  anode  is  below  and  the  cathode  above,  the  plane  of  the 
paper.  Furthermore,  we  take  the  plane  z  =  0  as  lying  midway  be¬ 
tween  the  mirrors.  Assuming  that  the  change  of  the  medium  is 
negligible  in  a  round  trip  photon  time  (2  L/c)  we  may  write  the 
following  integral  equation  for  the  complex  amplitude 


u  (x,  y,  z  =  0,  t^  +  — )  = 

f^op  '  y*  y"'  Z  =  0,  s,  t^)  U  (x',  y',  z  =  0,  t^)  dx'dy' 

(13) 

In  Eq .  (13)  u  represents  the  amplitude  for  either  the  left 

or  right  going  flux;  the  propagation  operator  K^p  is  a  function 
of  the  state  of  the  medium;  and  s  denotes  the  path  taken  by  a  ray 
from  (x',y',z'=0)to(x,y,2=0). 

In  the  presence  of  small  perturbations,  we  make  the  decompo¬ 
sition 


K 


op 


=  K 


op  (o) 


+  K 


op  (1) 


(t) 


for  the  propagator  and 


(14) 


u  =  exp  {i(k^z  -  w^t))[u^(x,  y,  z)  +  u^(x,  y,  z,  t) ]  (15) 

for  the  complex  amplitude,  whore  k^,  ii>o  ‘ire  tlic  wavenumber  and 
angular  frequency  of  the  light.  both  Uq  and  u,  arc  ay.r.umed  to 
vary  slowly  over  a  distance  comparable  to  a  wavelength.  By  defi¬ 
nition  of  Uq  ,  Kqp  (q)  and  of  what  constitutes  a  laser  mode  wo  liavt’ 
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u  (x,  y,  z  =  0,  t  +  — )  = 
o  '  '  j  '  '  o  c 


'op  (o)  z  =  0)  (x' ,  y',  z'  =  0)  dx’dy' 


Keeping  terms  to  first  order  in  Eq.  (13)  and  using  Eqs.  (14)  -  (16) 
we  have 


['‘o*'=o  “1  "=o  *  f’’]  ""P  ■■ 

+  /  V  (1)  “o  '^>'■<5'"  *  j'^op  (o)  ' 


K  /  1  u  dx'dy' 
op  (o)  o 


.  .  u,  dx'dy’  (17) 
op  (o)  1 


If  2Lu)o/c  =  nu  where  n  is  an  integer  (i.e.,  we  have  a  longitudinal 
mode)  then 


Uj^(x,  y,  2  =  0,  t^  + 


^)  =  Jki(VUq  dx'dy'  +|k^  Ujlt^)dx'dy 


Expanding  the  left-hand  side  of  Eq.  (18)  in  a  Taylor  series,  we 
obtain 


u^lx,  y,  z  =  0,  t)  +  — 


«  y* dx'dy '  +  J  dx’dy' 


2h  _ 
c  9 


[K  -  6(x-x')  5  (y-y')l  u,  (x' ,  y',  z'  =  0,  t  )dx'dy' 


ir  “  jiKo  - 


y,  y ' )  dx’dy' 


In  physical  terms  the  second  term  on  the  right-hand  side  of 
Eq.  (29)  represents  scattering  "into"  u^,  whereas  the  first  term 
is  leakage  "out  of"  uj •  If  is  an  approximate  eigenmode  of  the 
loaded  resonator  with  eigenvalue  close  to  unity,  the  first  term 
on  the  right-hand  side  will  bo  second  order  and  can  be  dropped. 

We  investigate  the  meaning  of  this  condition  in  greater  detail  in 
Appendix  A.  We  have  then 
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J  K^(x,  X' ,  y. 


y  '  )  u  dx  '  dy ' 
o 


In  most  applications  it  is  the  flux  perturbation  which  is 
experimentally  of  interest  and  yhich  couples  to  gain  and  density 
perturbations  in  the  medium.  To  obtain  the  flux  equation  corre¬ 
sponding  to  Eq.  (21)  we  multiply  both  sides  by  u^*  and  take  the 
real  part  of  both  sides.  We  obtain 


^  (Uo*Ui  +  u^u^M  =  ^ 


We  note  that 


l^u^*  ^  u^  dx'dy'  +  c.c.  (22) 


=  <<>0  +  ‘<’1  1^0  ^ll^  "  ^o^  *“o*^l  ^0^1*^  ^  °  (lu^l^) 


so  that 


<Pl  =  +  “o“l* 

We  now  assume  that  the  lowest  mode  Uq  is  independent  of 
x,y  (i.e.,  spatially  uniform),  so  that  u©  may  be  taken  outside 
the  integral  on  the  right  hand  side  of  Eq.  (22).  This  will  be 
most  valid  for  small  equivalent  Fresnel  numbers,  but  is  a 
reasonable  first  approximation  in  any  case. 


o  [  f  +  c.c.j 


34'  _  "'1 
9t  3t 


As  we  have  noted  will  be  directly  related  to  the  perturbations 
of  the  medium.  We  now  turn  to  expressing  in  terms  of  density 
and  gain  perturbations. 

Assuming  the  medium  reacts  slowly  in  a  time  2E/c,  the 
propagator  Kqp  can  be  written  in  the  form 

K^p  (X,  y;  X-,  y';  z  =  0,  s,  t)  = 


^x,y 

exp  |i^/  n(x,  y;  x',  y';  t,  z'  =  0,  s)  d  s| 

‘'x',y'  # 
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Where  the  path  is  the  line  segmented  geometric  path  from  (x',  y') 
to  lx,y)  and  n  is  the  complex  refractive  index.  An  illustration 
for  the  typical  points  is  given  below  in  Figure  5. 

For  the  unperturbed  medium,  we  have,  simply 

/•(x,y) 


op  (o)  2XL 


rix,yj 

^  exp  {  i  ^  as  I 


(x',y') 


Hence,  if 


n  =  n  +  n, 
o  1 


op  (1) 


K  ,  .  =  K  , 
op  (o)  op  ( 


For  sufficiently  small  values  of  (less  than  one  fringe  dis¬ 
placement  on  an  interf erogram)  we  may  expand  the  exponential  in 
Eq.  (29)  to  obtain 


K  ~  i  IZ 

op  (1)  X 


P  p(x,y) 

I  ds 

^  «  I 


(x* ,y ' ) 


'op  (o) 


Substituting  Eg.  (30)  into  Eq.  (25)  there  results  for  either  the 
left  or  right  running  flux 


^■^1  _  '^‘^o  2Tri 
3t  “  2L  X 


.(x,y) 


'op  (o) 


I  /  n^  ds  I  dx ' dy ' 

(  J(x'y')  ’ 


■  /V*  (o)  I  [ 


(x,y) 


ds|  dx'dy' 


(x-,y') 
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Figure  5.  Illustrating  the  Propagator  via  which  the  Medium  at 

(x- ,y' )  on  the  Nth  Pass  Influences  tlic  Light  Amplitude 
at  (x,y)  on  the  (N  +  l)st  Pass 
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where  Bg  is  the  Glcrdstone-Dale  constant  tor  the  nictliu:-.i  and  , 
is  the  density  of  the  gas  mixture  at  one  atmosphere  pressure. 
Noting  that,  by  definition 


n 


R 


and 


we  have,  finally 


(37) 


(38) 


x';  y,  y')  (x',  y',  t)  d.x’dy' 

Y>  y’)  9]^  (x',  y',  t)  dx'dy'J  (3  9) 

To  appreciate  the  physical  content  of  Eq.  (39),  wo  trensfor 
spatially  as  follows: 

(a)  In  the  y  (anode-cathode)  direction  we  take  a  finite 
Fourier  cosine  transform. 


TT  “  -IT 


/ 


aj(x,  x'; 


(b)  In  the  x  direction  we  perform  an  exponential  Fourier 
transform.  Since  the  steady-state  flux,  to'  nontero  only  over 
a  finite  region,  wo  get  convolutions.  To  finesro  this  proi^lem 
we  smooth  in  kj,-spaco  over  a  region  2/b,  where  b  is  the  length 
of  the  flow  channel.  That  is,  we  define 


-K  +l/b 
X 


IK  2/b 


<fi^  (K^)  dx 


'k  -1/b 

X 


We  may  now  cast  Eq.  (39)  in  the  form 


(40) 


^»lm,k 

dt 


c^o 

2L  ^“Rjnm',KK'  ^Im’K'  '''  ^ ^  * 

m'  ,K  ' 
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We  now  postulate  that  one  acoustic  inode  will  be  favored  above  all 
others,  and  that  this  mode  has  no  variation  on  a  scale  loss  than 
b  (in  the  flow  direction)  and  has  m  =  I,  say,  in  the  y  direction. 
The  reason  for  this  assumption,  as  well  as  a  prescription  for 
estimating  I,  will  be  given  below. 

Equation  (41)  now  becomes  the  pair  of  equations. 


d(() 


1  C(^ 

0,0  _  _ o 

2L 


dt 


1,0 


I  R  ^ 

(_  ox.,oo 


+  2;.^,  g^ 


0,0, 


(42) 


(where  is  the  medium  thickness) 
M 


d4) 


1 0  c<i' 

£  ,  o  _  o 


1 


£,o 


^££,oo  ^o 


dt  2L 

In  Eqs.  (42)  and  (43)  we  have  tacitly  assumed 


+  2£j^j 


£ ,  o 


(43) 


a 


Imm' ;oo 


mm ' 


(44) 


since  the  main  effect  of  a  gain  variation  is  to  change  the  flu.v 
in  the  same  spatial  manner.  We  note  further  that  Eq.  (42)  is  an 
equation  for  the  rate  of  change  of  the  output  flu.x  and  as  such 


‘R 


o£,oo  1,0 


(45) 


can  be  considered  to  be  a  time  varying  modulation  of  the  cavity 
coupling.  Equation  (40),  on  the  other  hand,  describes  the  varia¬ 
tion  of  the  "£"th  Fourier  component  of  the  flux  distribution. 

The  following  physical  interpretation  can  be  given  to  Eq .  (40) . 

When  the  laser  is  turned  on  there  will  be  some  acoustic  noise 
(density  fluctuations)  between  anode  and  cathode  which  can  be 
expanded  in  a  Fourier  series.  Each  acoustic  mode  constitutes 
a  pliase  diffraction  grating  for  the  incident  flux  Iq.  In  Appendix 
B  we  will  show  that  the  follov/ing  statements  can  be  made  about 
the  scattered  light: 

a.  The  acoustic  gain  .md  optical  modes  will  have  a  trans¬ 
verse  variation  with  period  a/n,  where  "a"  is  the  anode-cathode 
separation,  and  n  is  given  by 


n^  =  a^/XL  (46) 
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b.  A  plane  wave  scattered  by  an  acoustic  wave  will,  after 
one  trip,  produce  a  flux  oscillation  with  the  same  transverse 
wavelength  as  the  original  acoustic  wave. 

c.  The  variation  of  the  acoustic  perturbation  in  the  optical 
propagation  direction  can  be  neglected. 

Furthermore,  wo  note  that  if  we  take  Uq  as  a  constant  for 
the  lowest  mode,  then  the  higher  Fourier  components  uj^r  (recall 
4>lJ.  =  Uo*  u;)^£  +  Uq  uif,*)  are  orthogonal  to  Uq  and  it  is  reasonable 
to  assume  that  the  higher-order  modes  are  representable  by  sines 
and  consines  to  the  same  extent  that  the  lov;cst  order  mode  is 
representable  by  a  constant.  This  simple  picture,  incidentally, 
allows  us  to  make  a  realistic  estimate  of  the  influence  coefficient 
aRj),ji,oo  (henceforth  referred  to  as  a)  in  Eg.  (43). 

Technically,  the  quantity  a  enters  into  the  model  cejuations 
as  the  influence  coefficient  from  a  sinusoidal  density  variation  on 
one  pass  to  a  sinusoidal  flux  variation  on  the  next  pass  of  the 
laser  radiation. 

The  radiation  diffracts  through  the  standing  acoustic  wave, 
which  acts  as  a  phase  grating,  and  completes  a  round  trip,  rein¬ 
forcing  the  density  wave.  Thus  reinforcement  arises  by  having 
the  perturbed  radiation  in  the  same  physical  shape  (same  spatial 
Fourier  component)  as  the  density  wave.  The  reinforcement  is 
accomplished  via  heating. 

The  perturbation  in  the  amplitude  of  the  scattered  light 
(due  to  diffraction)  is 


u 


1 


1)  u  ~  i  u  5 
o  o 


(47) 


Where  Uq  is  the  incident  amplitude  and  6  is  the  (sinusoidally, 
spatially  varying)  change  in  the  optical  path.  The  change  in  the 
flux,  <{),  is  given  by 


u  *Ui  +  u  u  *  ~ 
o  1  o  1 


4,o6 


(48) 


Hence,  in  one  pass 


(49) 
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Whore  f-m  is  the  optical  length  through  the  medium,  Fg  is  the 
Gladstone-Dalo  constant,  and  is  the  density  of  the  medium 
at  one  atmosphere 


The  equation  in  which  a  enters  is  of  the  form 


^IS.  c  ^1£ 

From  this  equation,  the  change  in  one  pass  due  to  a  density  per¬ 
turbation  is  given  by 


2L/C 


By  comparing  Eqs.  (49)  and  (51)  we  obtain 

2'irR  p  £, 

~  q  o  m 


We  now  proceed  to  a  consideration  of  the  gain  and  heating 
equations  necessary  to  complete  the  model.  The  gain  equation  for 
the  flowing  CO2  medium  is  (fs  ^  saturation  flux) 


/l2  +  V 
yat  3x/ 


_  2  _  ^  +  p 

T  1 


Where  t  is  the  effective  deactivation  time  from  the  upper  lasing 
level  of  the  CO2  molecule,  including  the  effects  of  the  nitrogen 
in  the  gas  mix,  and  P  is  the  pumping. 

Once  again,  we  consider  the  perturbed  form  of  Eq.  (53) 

(v  is  the  mean  flow  velocity) 


+  V 


^I'^’o 


T)io  effect  of  the  flowing  term  is  '  lamp  the  growth  of  any  per¬ 
turbations  with  a  damping  time  of  t  .lu-  order  of  the  convection 
time  through  the  cavity.  This  is  of  t)ie  order  of  b/v  wlierc  b  is 
the  beam  dimension  in  the  direction  of  flow.  Hence,  Eq.  (54) 
becomes 
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(55) 


rgi 

vgi 

go(r-l) 

Where  r 

0 

s 

at  ~ 

T 

b 

^  ^0 

The  J2.th 

coiT^ponent 

of  Eq. 

(55) 

will 

be  of  the  form 

^^IJ-  ^  V 

rgi^  go(r-i) 

at  b 

T  T 

(56) 


To  complete  the  coupled  set  of  equations  we  use  the  fact  that 
perturbations  in  the  gain  and  flux  amount  to  a  perturbation  in 
the  heating  profile  in  the  medium.  This,  in  turn,  will  affect 
the  density  perturabations ,  completing  the  loop. 

The  acoustic  equations  for  a  flowing  medium  and  non- 
isentropic  flow  can  be  obtained  from  th''>  equations  of  momentum, 
continuity,  and  from  the  first  law  of  thermodynamics.  In 
Appendix  C  we  show  that  the  heating  equation  is  of  the  form 


Where  Pj  is  the  perturbation  in  the  power  per  unit  volume  given 
by 


^1  =  T  <Vl  ^1^0^ 


^l^s 


(58) 


where 


V  =  V-,  -  V 


1' 


(59) 


and  V3,  V2  represent  the  upper  and  lower  frequencies  of  the  lasing 
transition.  Combining  Eqs.  (57)  -  (59)  with  the  prescription  (Cf. 
Appendix  C) , 


9  V 
— —  — 

9x  b 


(60) 


(where  b  is  the  dimension  of  the  beam  in  the  direction  of  tlie 
flow)  we  obtain 


27  • 
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For  the  tth  component  we  have 


2  2 
~  3 

y 


Where  "a"  is  the  anode-cathode  spacing.  Hence, 

=  (-t)  V  <Vi!i  +  -V  ^rJ’ 


Equations  (43) ,  (56) ,  and  (63)  constitute  one  set  of  coupled 

equations.  To  find  the  normal  modes  of  the  system  we  use  Laplace 
transforms  of  our  variables 


e  dt 


etc.  This  allows  us  to  mal;e  the  replace.ment 


<})^(t)  ->  S  <1)^(5)  ,  etc, 


In  matrix  form  our  coupled  equations  become 


if  ♦ 
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Where 


-  eg  -p- 
L 


r  V 

S  +  -  +  ^ 
T  b 


(S+g) 


■/  \  2  2  /  ,\  2 


Equation  (66)  contains  both  the  normal  modes  and  the  instability 
criteria  for  the  system.  The  normal  modes  can  be  found  from  the 
secular  equation 


A  (S)  I  =0 


Equation  (68)  can  be  rewritten  in  the  form 


s"  =  0 


where  y  arc  functions  of  the  system  parameters.  The  roots  of 
Eqs.  (68)  and  (69)  arc  the  normal  modes.  An  instabi]ity  obtains 
if  one  of  the  roots  has  a  real  part  greater  than  r.cro.  The  con¬ 
dition  that  no  root  has  a  real  part  greater  than  zero  can  bo 
formulated  as  a  condition  on  the  coefficients  Vn  (the  Routh- 
Hurwitz  conditions).  We  now  proceed  to  a  deter  nu  nati on  of  the 
stability  boundary  (Kq.  (98)). 

C.  STABILITY  BOUNDARY 

Equations  (67)  -  (69)  yield  the  following  values  for 


*  =Ra'>ISb'  "  =Rb'>  -  fr  ” 


Avea  f  Vi  HL  1 1 
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Tho  Routh-Uurwitz  conditions  require  that  for  stability  the 
following  conditions  obtain; 


>  0  (n  =  1...5) 


Yo  -  Y, 


V4^3  -  ^2^5/ 


^0^5 


V  I  ""2^^ 


To  nake  the  algebra  tractable  wo  take  tho  leading  terns  in 
£qs.  (70)  -  (75).  In  determining  the  leading  terms,  we  note  that 
for  typical  systems 


Then  we  take 


V  q  q  2  „  2  P  q 

^0  ~  ^RA  ^RB  T 


9  P  2  2  ^  T 

Y,  w  +  2S  +  Pa  — 

‘  1  RA  RB  V  RB  V 


^2  ^Rb' 


>3  Sb 


.  ;A\/i::o  c  vtnrTi 


where 


/,  \2  q  1 

^  Ie  ^ 


Now  Eq.  (84)  is  equivalent  to  t)ie  condition 


2  2  2  2  3 

^0^2^3'^4  ^  '  o’’’!'*  4  “  '''  *  ' t'' "  ‘''i  >  -  ^ 


0 ' 3  '4 


1 ' 2 ' 4  '  '  3 ' 4 


1  '  4 


Inserting  Eqs.  v86)  -  (91)  in  E(j .  (93)  and  using  the  con¬ 

ditions  of  Eq.  (85)  we  obtain,  after  some  tediou  s  Q 1  y  t'  L)  r  ti  f 


{  0  r  ^pu  '1 

^  2S^^)  4S^  S3  ^  - ^_Kl^  (1^3-1  S^.) 


'' 


Or,  if  a  <  0, 


Pa  < 


2'=v’'  -"kb' 


1 


1  .  3(S^1)  -- 


The  Eqs.  (81)  -  (83)  arc  trivially  satisfied  excegt  for  the  co 
dition  that  y ^  >0.  If  a  >  0,  this  iinpos.es  the  condition  tiiat 


Pa  <  (S^i) 


KB 


2 


( 


Combining  Eqs.  (95)  and  (96),  we  have 


ll 

P 


S  1 
V 


1  +  3S  1 

V 


V 


< 


a 


S  1 
-Jii- 


e  2  ,,  2 

IIkA _ ^  Kh_ 

P 


as  our  stability  criterion. 
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utilizing  the  definitions  contained  in  Eqs.  (76) 
and  Eq.  (92),  v/e  have,  alternatively 


-  2 


2L  V 

b  • 


where 


’A 


1  +  3  —  ■ 

V 


2tiB  P_ 

-  <±— ^  —  £  (Y-1)  <  , 

VT  a  S 


-  (80) 


2L  V 
b  ’  C 


S 

(98) 


“=  f^o  S' 


(99) 


D.  DISCUSSION 

Equation  (98)  can  bo  put  in  a  more  useful  form.  Wo  note 

that 


'^SA  Ip 


(100) 


T  = 


'’o 


(102) 


o 


AA  p. 


(102) 


Where  t;^,  l',YA  ‘Jenotc  the  values  of  i,  /c;,  Fn-c.ssure 

of  one  atmosphere.  For  a  value  of  <  0  (for  wh.ich  t  l.e  loft  side 
of  the  inequality  of  Eq .  (98)  obtains)  tlie  limiting  value  of  ,  /. 

which  we  denote  by  (Pq/Pa^  cr  will  bo  given  by  °  ^ 


(103) 
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Let  us  consider  a  typical  system  v/ith  the  fol  1  ov.iny 
paramters 


P 


A 


0,6  kg  m 


v/C^  =  0.5 


4’S;^ 

X 


2 

1.6  X  JO  I'*  /  m 
10  ^  m 


(Y-1)  =  0.5 

e  =  1.7  X  10“ 

g 

Cg  =  360  m/soc 

b  =  7  cm 

T  =10  sec 

o 


(105 


Under  those  conditions  the  maximirm  a  1  lowaij  1  e  d-'n;  ity 

(pressure)  is  given  by  Eq.  (103)  as  ( ,■  /■.■)  -  0.38. 

o  A  cv 

It  is  interesting  to  note  the.t  (93)  drvs  not  contain 

the  ilux  explicitly  and  i..s  thus  i nd.ei. endont  of  than  hv'.']  of 
pumping.  V.’e  must  be  carc-ful,  liowevor,  to  iKite  tlsit  hi.  (98) 
holds  only  in  Eej .  (85)  liolds.  'J’liis  imjdjos  a  co:niiti<n  on  tiV' 

flux  exceeding  some  minimum  amourit  .  Also,  i'q .  ('tf)  a.';'un,i's  ti'.i.' 

substitution  of  Kgs.  (86)  -  (91)  fcjr  Kqs.  (70)  -  (75)  .  If  the 

latter  are  used,  Jq  will  enter  tlie  stability  conditi'nn  explicit! 

The  linear  instability,  if  it  exists,  will  as  we  lia'.'e  rann 
involve  an  exponent  i  a  1  1  y  erowinri  aconrtic  wave.  This  w  s/e  nee) 

have  e”'''.-  a-a’l  ^ . 1  M  ■  'e  I'avlt  in  1  .  s  w  :l;n  •'  s  ^  '  ,t-  •  . 

Ke  via\'  see  this  .sa a  1  1 1  a  t  .  v  i  y  :  :  e  1  m  i  ,  )  .  .  : . e  s  ’.  ;  :  t 

will  te..g  la  h:'  ales  ■  t'  t  '  ( i^nn- vt  a  11  -  ..-i'  :  '  m  ' 

tliey  liavo  a  small  cas  :  i  n  ;  tis.e  come,  am;  ro  t  :.<■  :ais.-;  lax  m  .  t.- 

(i.e.,  Sv  <<  Ey) .  Hence, 


1 1 
t 


■  1 


(in 


M  • 
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and 


IZ 


4> 


o 


a 

2L 


_1_ 

Ztt 


[21  g, „ 


(107) 


and  the  amplification  can  be  quite  significant  (in  practice 
~  103  -  104)  . 


Equation  (98)  suggests  difficulties  in  trying  to  run  a 
high-pov/erod  CO2  EDL  in  a  cw  mode.  A  solution  suggested  by  this 
analysis  is  to  damp  the  acoustic  mode  further  by  means  of  an 
acoustic  absorber:  for  example,  an  acoustic  muffler  on  the  cathode. 

This  has  been  tried  and  found  to  be  successful .  (2 , 3 )  a 

subsequent  paper  we  will  investigate  the  effects  of  an  acoustic 
absorber,  as  well  as  more  realistic  kinetics  and  compare  our 
results  with  experiments  undertaken  at  AERL. 


2.  Kellcn,  P.F.,  Mattson,  A.C.,  Ahous.c',  D.R.,  and  VodL'r,  M.J., 

Optical  Engineering  3^.0  (1979). 

3.  Korff,  D.,  Glickler,  S.l,.,  and  Dau'ih.orty,  J.P.,  "Acovistic 
Instability  Model  for  High  Power  ew  EDI.  L.trers,:  30lti  Annual 
Gaseous  Electronics  Conference,  Palo  Alto,  Cal.,  Oct.  18-2-1, 
1977. 


.  ’AVCO  1  V!  Ml  n 


SECTION  III 


ALIGNMENT  SENSITIVITY 


In  tho  previous  section  v.'o  have  ignored  c.iiri  \':i  r  i  a  t  i  or. ;;  in 
the  direction  of  flov.’.  Such  viiriation.s  can  load  to  extreire  alicn- 
ment  sensitiv  i  ties ,  as  we  shall  now  siiow. 


There  arc  two  effects  involved.  The  first  i ;;  th.e  \’nriatior, 
of  output  coupling  v.-ith  position  of  the  optical  axis  on  tl'.e  fc;‘ed- 
back  mirror.  This  is  basically  a  diffraction  effect  (“1)  and  is 
illustrated  in  Figure  6. 


In  Figure  6  N^q  is  the  system  equivalent  Fresnel  number, 
defined  by 


N  =  J'l 

eq  2.\L  M 


(108) 


'..•licrc  a  is  the  colllmatina  mirror  diam.'ator,  M  is  tiie  systi,  :n 
magnification,  and  t'nc  other  symbols  ha\-o  be(_n  ■ 'Xi.)!  a  i  ned  aiiove. 

£  is  defined  as  the  normalized  position  (one  dimonsional)  of  ti.e 
optical  axis  on  the  feedback  mirror  ('■  ==  ■»■  1  ,  0,  -  1  correspond 
to  top,  middle  and  bottom  of  the  mirror,  respectively). 

The  second  effect  is  the  variation  of  the  syste.m.  cain  with 
position  of  the  optical  axis.  This  effect  arises  boca'.ise  of  the: 
variation  of  the  aain  with  temperature  for  a  particule.r  \’alue  of 
J,  the  initial  annular  momontu,."  of  t.he  j-adiat.i.ng  molcculC'.  This 
variation  arises  through  the  partition  function.  In  fact. 


ga  in 


2J  +  1  -BJ (J  +  1 )  /T 
e 


(109) 


where  B  is  the  vibrational  constant  for  the  transition  and  T  is 
tho  temperature. 

Our  flux  cqvintion  (for  idu?  variati(,>n  of  tli"  flux, 

not  a  Fourier  coniponcnt)  is  (x  is  tlie  flow  direct  ion,  aiuj  x  =  0 
is  the  optical  axis) 


94’ J 

TtT 


c? 


(g 


Ui) 


(110) 


4.  Kell  on,  r.  ancj  .‘’'.mitli,  M.  ,  Oj  »t .  Eng .  ,  18,  pp.  1’>7-]G0 

(March/April  1979). 
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where  Pi,  proportional  to  the  tilt  angle,  is  defined  by 


U 


1 


? 

3x  ^1 


(X) 


0 


=  dimensionless  coupling  constant 


(111) 


We  shall  derive  below. 

The  gain  equation  along  the  optical  axis  (or  the  average 
system  gain  across  the  beam,  which  can  be  shown  to  the  equivalent) 
satisfies  the  equation 


_  2  _  ^Vi  ^  ^I'o*  + 

9t  T  T(f)^  dy,  dt 

s  1 


(112) 


In  Eq .  (112)  the  last  term  on  the  right-hand  side  repre¬ 

sents  the  effect  of  the  gain  variation  with  optical  axis  position. 
We  rewrite  Eq .  (112)  in  the  form 


Sg,  (g  d-T  +  gif'  )  c«. ^bg  ?y, 

^1  _  _  3  _  ^o  1  ^1  o  2  -’o 

3t  "  T  Kfg  at 


(113) 


where  a2  is  a  dimensionless  constant  whose  value  will  be  found 
below. 


The  heating  equation  talces  the  form 
3t  a^^)  p^(x)  =  (Y  -  1)  a^f(x)  (g^4.^  +  g^P^)  (114) 


Where  f (x)  denotes  the  variation  in  gain  along  the  flow  direction 
(g^  and  g^^  are  assumed  to  have  the  same  spatial  dependence).  The 
equation  for  yi  ("  p]^  '  (0))  then  becomes 


at^t 


-  c. 


a,  )  P-, '  (X) !, 


=  (Y  -  1)  9^  ffx) 


X  =  0 


(g. 


A 

i"  1 


where  the  prime  denotes  derivatives  with  respect  to  x. 
Henceforth,  we  take 


(115) 


(116) 
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to  simplify  tho  analysis.  If  a  standinq  v.mvo  oxist:;  in  thi-  f  1  o-..- 
channel  (analoqous  to  an  open-ended  oryan  pip'e)  tlK-n  I'.q.  (llfd 
follows  naturally.  In  any  event,  it  serves  as  a  reasonable  ap¬ 
proximation.  Similarly  wo  take 


(Y  -  1)  3^  f(x)|^  ^  0 


<-  (Y  -  1) 


(117) 


where  03  is  derived  below. 

Our  three  coupled  equation  for  this  case  are  given  by 

/^1 

A  (  gi  (s)  )  =  0 

(S)  ' 


(118) 


where  (r  =  1  +  -t—  and  v  is  the  gas  flow  velocity) 


-  c-f. 


V  r 

s  +  T-  +  — 
b  T 


a  jb 

eg  (f - 

°  °  Pq 


bg 

■'o 

-  a^s 

P  2 


a^q 
3  -’o 


(119) 


(Y  -  1) 


3  ^o 


(Y  -  1)  -  (s  +  f)  (s  +  ^)^  + 


A  Routh-Hurwi tz  staljility  criterion,  using  leading-order  terms 
as  in  the  previous  chapter,  leads  to  the  stability  requirement 


where 


,  ,  ,  2  1 

“l«3  ^  ^  ’  T-  ■  ^g~ 
^s  c  s^o 


^  3  ,  b 

0.  -  p  C  and  =  — 
’^A  o  s  c  V 


(120) 


(121  ) 


V'Je  shall  now  derive  expressions  for  and  (wliich,  as 

wo  shall  see,  is  ~  n2) . 
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To  dcteriuine  we  must  relate  chancjes  in  c  to  chancres  in 
p*  (x  =  0) .  We  have 


2iL  *  =  lii  ft  J_  0 

A  ^tilt  A  p 
^  a 


(122) 


where  is  the  thickness  of  the  optical  medinm  and  !tilt  is  tlae  tilt 
of  a  ray  due  to  refraction.  The  optical  axis  will  then  tilt  by 


°trlt 


(123) 


Furthermore  , 


Ac  - 


(124) 


where  are  the  diameter  and  radius  of  curvature  of  the 

feedback  mirror,  respectively.  Hence, 


-6  • 

d.  M-1  ’^g  M  p. 


(125) 


Now,  let  us  denote  the  change  in  g^.  when  c  changes  by  (l/N'eq)  as 
ngo-  Then, 

Ag  =  II  Ac 


=  (ngo)  Ac 

=  (ng„)  •  - 


o  eq  d,  M  M  -  1  p. 


Ap'  (0) 


(126) 


«igob 


Ap' (0) 


Therefore , 


^  „  _JM _  ^  ^ 

““l  ^eq  M-1  b  d^  "g 


For  a  typical  ca.se 


~  n  1  M  in  M  n  1  20  =  200 

n  0.1,  10,  M  =  2,  Q-j- 


(127) 


=  0.1, 


2  X  10 

d  2  cm 
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C.)  < 


so  that 


«  a 

We  note,  furthermore  that  nj  will  be  much  less  than  its 
alculated  value  if  the  system  is  aligned  at,  e.g.,  points  A  or 
'  in  Figure  6.  (The  calculation  in  Eqs.  (126)  and  (127)  is 
for  point  B  in  that  figure.) 

We  may  schematically  denote  the  variation  of  system  gain 
with  optical  axis  by  the  diagram  in  Figure  7. 

An  analysis  identical  to  that  above  leads  to 


-  n'  M  4M  ^  r  !!o 

^2  ^eq  M  -  1  b  d,  9  P-, 

1  a 


(128) 


Now,  0-2  arises  from  the  temperature  variation  of  the  gain 
as  given  in  Eq .  (109). 

Talcing 


T 

inlet 


=  173°K, 


=  380° 


we  obtain 


02  *  0.5 


for  a  typical  case. 

Since  02,  013  arc  both  dimensionless  constants  due  to  the 
same  effect  we  take  them  to  bo  equal.  Wo  then  obtain  the  follow¬ 
ing  stability  boundary. 


<^eq’  [(M  -  Dbd  J  ‘‘g  \  j 


2n  C' 
o  s 


C  g 
s  -^o 


(129) 
^  c 


For  a  numerical  example,  tak 


M 

=  2 

«1 

n 

0 

II 

n’ 

=  0.1 

^1 

^o 

0.5  m“^ 

^s 

20 

N 

cq 

20  m 

b 

1  0  cm 

0.05 

- 

1  m 

-  4 

2  cm 

2  X  10 

300  m/sec 

T  = 

c 

10  ^  sec 

SYSTEM  GAIN 


In  such  a  case,  we  have 

^  4  (unstable) 

2  b  J 

An  instability  in  this  system  w'ould  result  in  larc?e  wander¬ 
ings  or  the  pointing  angle  and  possibly  to  lino  ji-imping  if  the 
optical  axis  moved  sufficiently  far  so  that  the  gain  were  greater 
for  a  different  J  number. 

Even  in  the  absence  of  an  instability,  of  course,  the 
matrix  A  determines  the  system  response  to  external  perturbations 
In  particular,  if  v;e  have  mirror  vibrations  (n^  perturijation)  th.e 
coupling  ensures  that  we  would  obtain  a  response  at  the  acoustic 
frequencies.  Such  a  response  v/ould  not  obtain  in  standard  treat¬ 
ments  . 


Finally,  we  stress  that  if  the  system  is  aligned  at  a 
stable  point  (A  or  A'  in  Figure  6)  no  iiistaLility  obtains,  .i 
(or  n)  is  negligable  and  Eq.  (129)  is  a utomat ica 1 ly  satisfied. 


■  AVCO  i  VLML  r  I 


SECTION  IV 


MODE  MEDIUM  INTERACTIONS:  DIFFERENTIAL  FORMULATION 


A.  INTRODUCTION 

The  analysis  in  this  section  will  deal  mainly  with  an  attempt 
to  recast  the  instability  discussed  in  Section  II  in  terms  of  more 
conventional  medium  instabilities  such  as  are  found  in  plasma 
physics.  Some  of  these,  such  as  the  free  electron  laser,  harmonic 
generation,  and  optical  m.odulation  by  sound  waves  are  beneficial. 
Others,  like  thermal  blooming,  self-focusing  and  mode-medium 
interaction  (MMI)  are  deleterious. 

The  analysis  in  the  previous  two  sections  emphasized  the 
Huygcns-Fresnel-Kirchhof f  integral  solution  to  the  v.’avo  equation. 

In  this  section,  instead,  the  differential  equations  wiil  bo 
analyzed  and  comparison  with  the  other  derivation  will  be  made 
along  the  way.  The  differential  formulation  has  the  advantage 
that  it  can  be  easily  extended  into  the  nonlinear  regime,  in  vdiich 
the  amplitude  of  the  perturbation  has  built  up  to  a  large  value. 

It  also  has  the  advantage  that  it  demonstrates  how  the  instability 
discussed  in  Section  II  (hereafter  referred  to  as  MMT-1)  is  related 
to  more  conventional  instabilities.  To  do  this.  Maxwell's  equations 
including  a  time-varying  conductivity  (to  account  for  cavity  and 
medium  losses)  and  a  time-varying  dielectric  function  will  be  used. 

It  will  be  seen  that  MMI-1  results  when  the  conductivity  is  modulated 
and  when  all  the  electromagnetic  v;aves  have  the  sam.e  unperturbed 
frequency.  This  is  a  form  of  nonresonant  amplitude  modulation  and 
instalji  1  i  ty ,  which  is  different  from  the  convt>ntional.  m.ethod  of 
amplitude  m.odulation  in  which  resonant  coupling  of  electromagnetic 
modes  which  are  separted  by  the  modulation  frequency  occurs.  On 
the  other  hand,  Brillouin  scatterinci,  which  is  primarily  a  phase 
modulation,  results  when  the  dielectric  function  is  m.odulated. 

The  coupling  is  strongest  wh.en  the  electromagnetic  mocLes  are 
separated  by  the  modulation  frequency.  This  is  the  mechanism  that 
will  bo  studied  in  this  section.  The  cq'aalions  of  motion  for  tl;e 
electroinagnc't  ic  and  acoustic  waves  are  derivi'd.  The  secular 
determiiiant  is  s.olvod.  It  is  shown  to  have  unshrihle  rerg  in  the 
limit  in  which  the  ratio  of  cavity  flux  to  the  r.at  urat  i  on  flux  is 
small,  so  that  cjain-flux  o:;ci  1  lat  ions  can  be  necjl  ectc'd .  Finally, 
the  nonlinear  regime  is  brie^'ly  examined. 

D.  MODE-HF-DIUM  INTERACTION.^ 

Our  starting  point  v/i  1  1  be  Maxwell  's  e<!ual  ieui;;  v.’ i  {  h  liu' 
appropri,ate  boundary  conditions  (cxivity  lonc.ru;)  arul  ;■  a.’  i  i  u;-  'i.iin 

m 
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and  losses.  Since  the  laser  is  in  an  unstable-resonator  config¬ 
uration,  the  boundary  conditions  arc  rather  complicated  (sec 
Appendix  D) .  Thus  to  simplify  matters,  we  will  lump  all  the  losses 
in  a  phenomenological  conductivity  function.  The  medium  gain  will 
be  computed  semiclassically  from  the  imaginary  part  of  the  suscep¬ 
tibility.  A  simple  model  of  this  is  discussed  in  Appendix  E. 

Thus,  the  governing  equations  are,  in  mks  units. 


V  •  D  =  0 


(130) 


B  =  0 


(131) 


V  X  E  =  - 


(132) 


V  X  H  =  J  + 


(133) 


where  D  =  cE,  J  =  oE,  c  contains  the  effects  of  dielectric  polari¬ 
zation  and  0  contains  both  the  medium  losses  and  the  cavity  loss. 
From>  Eqs .  (130)  -  (133)  t)io  following  wave  equation  can  be  derived 

for  transverse  electromagnetic  waves; 


V  E  -  u 


0-at2 


(134) 


(Mq  is  the  froe-space  permeability) .  Assume  that 

c  =  (r)  +  c  (r,t) 

o  =  0  (r)  +  o  ^^Nr  ,t) 


(135) 

(136) 


The  entire  scattering  calculation  can  actually  be  done  in  terms  of 
a  or  1,.  The  only  reason  for  keeping  botli  is  to  stick  to  the  con¬ 
vention  which  associates  the  gain  of  the  system  with  the  im, aginary 
part  of  c  and  the  lossc-s  with  the  real  part  of  :  .  Similarly, 
scattering  is  associated  with  the  real  part  of  >  and  not  the 
imaginary  part  of  a  ev'^n  though  formally  (with  tlie  inter¬ 

pretation)  it  can  be  calculated  in  terms  of  eitlu-r  function.  In 
what  follows  we  maintain  the  convention.  Substi  tviting  Eejs.  (135)  - 
(137)  into  Eq .  (134)  leads  to. 


( 0 )  ,  ,  ^  E 

(J  fv  e  '  '  ( r )  —  „ 

3tr 


(0)  fdl  _ 
at 
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This  equation  will  now  be  solved  in  a  perturbai  i  ve  f  .ish  i '  i; . 
However  rather  than  expanding  in  normal  moder. ,  v/.  ■  '.-.■ill  vir.'^  tii'^ 
multiplc-time-scalc  expansion  of  secular  equat.i  ■ 'P.:; .  ^  t )  hxiandii 
E  =  L'O  +  e(1)  +  realizing  that  E^O)  ;  0  since  there  is  no 

steady-state  space-charge  field,  leads  to 
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^0^  2  5t  ^  ° 
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(139) 

If  and  0^^^  in  turn  depend  on  the  electric  fields,  tlu  p  third, 

and  higher-order  equations  will  result.  Equation  (13S)  dc'scnljcs 
the  equilibrium  eigenmodes  of  a  loaded  resonator  including  all  tlie 
losses.  This  equation  for  the  unstable-resonator  aoomet  I'v  is 
discussed  in  Appendix  D  and  in  more  detail  in  anotlu'r  rut  of  thi:; 
report.  'Die  normal  modes  of  e(2)  arc  tlio  same  as  tliose  for  i:(l) 
(since  the  left-hand  sides  of  the  eejuations  are  i  dicnt  i  r  a  ]  )  . 
However,  since  e(2)  ig  driven  by  products  of  fi  rr.  t -ord.or  quan'itic 
it  is  conceivable  that  the  left-hand  side  is  driven  at  a  natural 
frequency.  If  this  should  happen  to  be  t!ie  case  tlu-n  our  p'.  rtur- 
bation  expansion  will  broa!<  down.  The  reason  is  that  i;^^)  ronta;:, 
a  secular  response  v;hich  will  make  e(2)  >  e(1).  g-Q  cc>t  arc  in. d.  t 'n  i 

problem  those  components  of  the  right-hand  siche  of  t:(3)  pj-c.  inc-  p- 
porated  into  the  equation  for  eH),  v;hich  is  now  allowed  to  liave 
a  slowly  varying  amplitude.  First  the  eiaenmodc'n  of  i:u.  (lib)  liav 
to  be  found.  Towards  this  end  lot. 


E^^^  (r,t)  = 


(1)  -iai^t+ikz 
Ej^  o  +  c  •  c 


(140) 


whore  c-c  stands  for  complex  conjugate  and  F 


;i ) 


sat  i 
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whore 


1,2^2  (0 
,  k  C  .  q 

^  ~  2  ~  ‘^‘rO^O  ^  k' 


(142) 


g  ,  are  the  equilibrium  gain  and  loss  terms.  The  loss 

includes  medium  and  output-coupling  losses;  it  is  related  to 
o(0)  by,  g^(0)  =  (r)cjg^/k.  The  re  lat  i  onsli  i  betv;een  the 

dielectric  function  and  the  medium  susceptibility  is  given  in 
Appendix  E.  Tlie  dipole  moment  per  molecule,  (liven  by 
p  =  (aE  +  C"C),  results  in  a  macroscopic  polarisation  P  =  '-qP- 
Finally,  since  d(1)  =  +  F,  we  fiiid  tlint  in  ov 

simple  semi-classical  picture  t(0)  =  cod  +  ‘KQld  “  iu(0)/ki. 
The  equation  for  the  equilibrium  gain  is  given  by 


V  .  g(0)v(0)  = 


+  9_- 


(0). (0) 


(143) 


where  V'  '  is  the  equilibrium  flow  velocity  of  the  c)as ,  •  is  the 

collisional  deactivation  time  of  the  upper  state  a.nd  P  is  the 
pumping  rate.  The  solution  to  Fqs.  (142)  and  (143)  conr,tituteB 
the  eejuilibrium  eigcnmodcs  of  the  resonator,  since  thi'  resonator 
boundary  is  included  in  c(0).  The  details  are  given  in  Ap;  i.-ndix 
D.  Any  perturbation  in  the  equilibrium  fields  can  bo  c-xt  andod  in 
the  complete  or thonorir.a  1  set  of  cigonfunct  ions  (v.diich  arc-  as.sura  d 
to  exist).  h’ov.-GVor ,  since  a  secular  rospon.se  is  ant  i  c  i  t  ate: ,  tin 
Fourier  coefficients  will  be  assum'^d  to  bo  tiiao  dependent  anui  tin' 
appropriate  set  of  equations  dcscriljincj  the  time  evolution  of  tlie 
coefficients  will  now  be  derived.  That  is,  let 


E(r,t)  =  i;  u^.  (t)E^_^^'  (r)o 


-ii.  ),  t  +  i  kz 


(144) 


and  sul'S,  t  i  tout  e  t i  s  into  I’q .  (137),  keeping  in  rd.n.l  1!m‘  '.:i-(t)  is 

slov.'ly  varying  and  i:p(l)  fnitis.fies  i;q.  (141).  Ti\'‘  }'')tui!  at  ion 

expansion  in  Eqr..  (138)  ajid  (139)  thus  i;;  r*  don'-  to  take  ir.to 
account  antici[iatcd  secularity  in  t)io  re.sponnv  of  e(2)  .  Avnuimo 
also  that 


r<^Nr,t)  =  Ladt)./^Nr)e"'"^^^"'"^ 
k 


(14'.) 
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a  (r,t)  -  ^  b.  (t)  '  (r)t 
^  K  K 


(146) 


where  'y  (r)  and  <' y  (r)  arc  the’  aj'jjroj  rieitc’  siKitial,  dieleetrii 
and  loss-modulation  e  i  aon  fund  ions  to  be  tii't  ermi  ned  sel  f-cc  r.:' i  s- 
tontly.  Note  that  in  F.qs.  (144  )  -  (1  46)  k  is  a  sum^-iat  i  on  i!'.';  ; 

there  is  no  implication  that  tt\c  k's  are  f  h.e  same  in.  there  c  '.pia- 
tions.  Substituting  Egs.  (  144)  -  (146)  into  I’q .  (137)  yield.;;, 


-iuut+ikz  2i.a',  du 

k  ^ 

g2  -im^.t-fik'z  -io^„t  +  ik"z 

*^'0  2  ^  ^k'^k'  ^  '^k^'^k" 

k*  k" 


3  y  ,  (1)  .t+ik'z„  -i(i'j^„t+ik"z 

^0  3t  ^  ^ 
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Not'"'  ♦hat  the  secular  i  ty-oroduci  ng  terms  lead  to  slow  variations 
in  the  mode  am.plitudcs  as  expected.  To  pick  out  a  ]  articular 
eigenm^ode  on  the  left,  m.ultiply  the  left  and  riab.t  sides  by  idie 
complex  conjugate  of  the  left  and  integrate  over  the  resonator 
volume  and  time: 


3^k  ic^  f 

^  /  dvdt  exp[i  (Wj^-Ui^^.-Cj^,,)  t-i  (k-k'-k")  z] 

^  k  '  ,  k  " 


X  (wj^,  +  o;j^„)  ko^k'^k"^k  ^  ^k"  ^““k' 


+  S -  S  /dvdt  exp  [  i  (o,  -U3.  , -01.  „  )  t- i  (k-k  ' -k  "  )  z  ; 
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2  ( 1  )  *  ( 1 ) 

where  t)ic  cigenmodcs  are  normalized  such  tliat  /d 


.  ■  Avea  L  VI  ( is.  r  T 


There  are  several  possible  cases  of  ititere.st  : 

(i)  Consider  first  the  case  in  which  there  is  no  dielectric 
modulation  (cj^i  (r)  =;  0)  . 

Taking  the  limit  that  (.:p.  e  (the  slow  modulation  of  the 
gain  and  loss  term)  compared  to  optical  f  req  uenc  i  o;'  is  neglicjiblc 
recjuircs  tJiat  That  is,  the  two  el ect  roneuintut  i  c  '.-.■aves 

have  the  same  natural  (unperturbed)  frcquejioy  of  oscillation. 

One  way  this  could  haj-’P^-n  would  be  if  k"  k,  i.e.  ,  the  lasina 
mode  is  "coupled"  to  itself  by  loss  modulation  arisinn  from  the 
acoustic  wave.  Alternatively,  the  tv.-o  niodes  may  he  dis,tinct, 
although  degenerate  in  frequency;  this  case  is  treated  in  detail 
in  Section  II.  In  tlic  sel f -coupl ing  case  we  obtain 


3u^.  2 

^=C  UQ(bj..O  )Uj 


V. 


with  a  similar  equation  for  Uj.„.  Making  the  follov.-ing  identifica¬ 


tion 


aR^p^o 


(r)  '  (r)  (r) 


results  in  the  equation  used  in  M.MI-1  witli  being  the  density 
pcrturbatioii  by  the  acoustic  wave.  The  gain  term  can  be  obtained 
from  the  polarization  term  as  is  done  in  Appendix  D.  The  above 
identification  of  fn.i  (1)  with  density  f  1  \ictviat  i  or.s  is  not  unrea¬ 
sonable  since  the  ca\'ity  and  medium  Ic-ssc's  for  sir-.plo  systems  can 
be  related  to  the  quality  factor  of  the  system  which  in  turn 
depends  on  the  index  of  refraction.  Besides  the  case  '/v.  =  :.]^n  it 
is  also  possit:le  to  consider  the  resonant  triplet  .  =  .qq’  +  ' 

k  -  k'  4-  k "  where  now  the  two  eloctromafpK'ti  c  "'.odes  are  scuiarated 
by  the  acoustic  frequency.  The  only  difference  is  that  now  the 
-iWk ,  t 

e  term  in  I’q.  (149)  is  rci''lacc'd  by  unity.  The  right-hand 

side  is  again  related  to  the  density  fluctuation.  This  is  the 
case  of  conventional  amplitude  modulation  by  a  resonant  triplet. 

(ii)  Consider  nov;  separatf  ly  the  case  of  no  cavity  and  medium 
loss  variat  icjn  ('r'  0)  but  nfr’  inclvuio  a  dielectric  :re;li;l  at  ion 

fp'  9'  0.  I'c;  procur'd  furtlic)'  v.-e  dr'Vi'lop  a  sirqjli.'  '‘'U'  ^  v  "  . 

The  sp'r’c:ific  modulation  of  the  medium  that  is  to  )ic  studied  is 
that  due  to  an  acoustic  wave. 
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From  microscoiiic  considvrcitions  (soc  Appcnidix  E)  it 
follov.'s  that  tho  dipole  mor’.ent  of  a  molcculo  is  given  by 
P  ~  (''i’cff  +  C'C)  v.'hc-TO  a  is  a  f  reciuency-dc'pc'ndcnt  p>ro!jor  t  iona  1  i  ty 
constant  and  =  F.  for  a  qas.  Tho  m^icroscopjc  po  lar  i  za  t  ion  (  8  ) 

of  the  medium  is  then  given  by  (for  N’  molcculos/cubic  centimeter) 

ff  =  Np  (150) 

with  the  dipole  moment  given  by  p  =  (ctE  +  c-c).  Note  that  a, 
both  depend  upon  the  intensity  of  the  field.  Since  the  secular 
part  of  the  response  comes  from  products  of  first-order  quantities, 

using  -  c<Ue<1)  =  -  i  ^  ‘  +  cc  leads 

to 


3t  . 


+  co,„prg,<l'(t,r)E/<l'E,,.'l> 


(151) 


where  it  was  assumed  that  the  fast  acoustic  exponential  variation 
has  been  extracted  via  the  requirement  that 

u  =  (152) 

k  =  k'  +  k"  (153) 

and  the  fast  transverse  variation  in  and  is  taken  care 

of  by  the  electric  field  variations  which  have  different  transverse 
mode  structure 

k,;  =  k^-  +  k^"  (154) 

where  k  is  fixed  by  the  acovistic-rosonator  goom.otry.  In  the 
plane  wave  case  the  spatial  integrals  can  actually  be  done  and 
lead  to 


3u 


k 


at 


(1) 


(155) 
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^’g  ~  tono-Da  1  c’  coef  f  ic  i  en  L  ,  ia  ^  tt:a  atr.'.o- 

sphoric  dci.sity  and  *  dc'noti's  complex  coniniate  nt  i  {. :  ex  . 
Equations  for  qy  <  ,  Nj^,  will  be  derived  in  ti.e  next  xc  ctiox. 

The  only  term  whicli  can  drive  tJ'.c  cain  equation  rcnonantl' 
[that  is,  so  that  Kes.  (152)  -  (154  )  are  satisfied]  is  the  flux- 
term.  Thus,  the  secular  part  of  the  gain  equation  is 


Yg^k- 


^0  T 


(r-1)  ^k^k” 


(157) 


where  r  =  1  +  -0  -  cto'Jk^/2  and  (-,0  =  1/-  +  ('.■.■i'ert 
Vq  is  the  flow  velocity  and  b,  the  anode-cathode  distaiicc)  in,  i'u 
total  damping  of  the  gain  term  including  convection  and  deactiva¬ 
tion  of  the  upper  level. 

The  equations  dcscribiiig  the  medium  arc  the  continuity, 
momentum  and  heat-transpert  equations. 


-f  V  •  pV  =  0 


(158) 


|v  V  •  vv .  -  -':e 
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+  V  .  VP  =  t  V  .  Vp  +  (-, 


(160) 


where  p,  V,  P,  dQ/dt  are  the  ma.ss  density,  fluid  velocity,  prcBsuix 
and  heat  deposition  rate;  y  is  the  ratio  of  specific  heats. 


Linearizing  Eqs.  (1  58)  -  (160)  would  produce  tin. 
acoustic-wave  equation.  How-evor,  since  the-  heat  detxisit- 
cari  have  cnmjionc'nts  at  the  acoustic  f  rcxiu-vicy ,  111'.:',  lead: 

secular  ref,p(insc  in  the  acous.t.ic  equation.  M'h  i  s  c.ecular 
once  aejain  taken  earn  of  by  Icttincj  the  am.pl  i  tvuler;  In'  s,  1 ' 
varying.  The  linearized  acoust i c-wavc  equation  (fer  tlu' 
density)  ir;  given  by 
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ui  tL'i'm 
V  o  a 


it  "  ^Oi: 
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where  Vq  is  the  transverse  flow  velocity  and  C,i“  =  I RT  ^  the 
acoustic  velocity.  The  heatiny  teria  h.as  beem  incUuied  in  antici¬ 
pation  of  the  sccularity.  The  heatinn  term  is  obtciiii'-d  few  a 
laser  by  consorvation-of-eneryy  arguments;  it  takes  the  following 
form : 


(1C2) 


whore  is  the  ratio  first  excited  state  to  ground  state  and 

c*c  stands  for  comprex  conjxigatc.  Since  P  docs  indeed  '.’ary  at 
the  acoustic  frequency  (note  that  use  of  Eqs .  (152)  -  (154)  was 
made)  it  v.'ill  produce  a  secular  response  in  Eq .  (161).  The  slov.-ly 

varying  piece  of  the  acoustic  wave  is  produced  by, 


3t 
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v/hcrc  3  a  is  the  convective  damping  of  the  acoustic  waves,  is 

given  by  Eq.  (162)  and  p;.^  >  =  mNj.;<  where  m  is  the  mass. 


The  linear  instalaility  is  obtained  by  assuming  that  one 
of  the  electromagnetic  waves  has  largo  amplitude  and  hence  decays 
into  the  second  ei  cctromagnetic  w’a'/o  and  the  acoustic  oscillations. 
Thus  assuming  the  (q  constant  and  using  Eqs.  (  156)  ,  (157)  and 

(163)  leads  to  the  determinantal  equation  for  the  instability. 
Assuming  that  Uj^„  -  e^^t,  ,  ~  e^*t,  g^  ,  ,  e-*^  loads  to. 


Setting  the  determinant  of  Eq .  (164)  etjual  to  ;nTO,  Icadr,  to  the 

polynomial  equation  for  the  instability.  Uefinim 
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1)  ,  a2 


(y  -  1) 


eg 


0  V, 
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Pq  u  -  1)  y  -  1  / 

c? 


where  -  Gladstone-Dalo  constant  and  g  is  atmospheric  density, 
we  obtain  ^ 


(fi)  (6  +  Yg)  (6  +  Y^)  +  a^(6  +  Ya)  "  ^  ''g^  +  iS  =  0  (165) 

The  acoustic  damping  is  given  by  Ya  ==  whore  is  the  trans¬ 

verse  gas  flow  velocity  and  b  is  the  anode-cathode  distance.  The 
damping  in  the  gain  equation  is  given  by  Yg  =  ^Q/b  +  1/t  where  t 
contains  all  the  collisionai  relaxation  rates. 


C.  NONLINEAR  SATURATION 


If  now,  it  is  assumed  that  u^  docs  not  stay  constant  but 
decreases  significantly  (as  energy  is  transferred  to  the  other 
v;aves)  ,  it  is  necessary  to  include  an  evolution  equation  for  u^. 
That  is  Eq .  (155)  must  be  included  together  v;ith  Eqs.  (156)  , 

(157),  and  (163).  Unlike  the  thre.o-wavc  equations  of  plasma 
physics,  which  can  be  solved  in  terms  of  Jacobi  elliptic  functions, 
this  set  of  equations  can  only  be  solved  numerically.  Some  insight 
into  these  solutions  can  bo  obtained  by  oxamininej  the  conserved 
quantities.  Using  Eqs.  (155)  -  (157)  and  (163)  leads  to. 
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This  can  further  be  reduced  by  using  the  equation  for  g,  ,  . 
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The  convective  damping  has  been  included  fjimply  by  replacing 
VqI/'X  in  the  acoustic  equation  by  Vo/b  where  b  is  the  anode- 
cathode  spacing. 

In  order  to  understand  Eq.  (167)  it  is  necossarv  to  examine 
Figure  8.  During  the  1 inear- lasing  stage  the  second  electromag¬ 
netic  wave  is  negligible  and  the  acoustic  noise  has  not  made  a 
transition  in  the  anode-cathode  direction.  Once  the  primary 
electromagnetic  wave,  uj^  has  boon  saturated  and  reaches  stocidy 
state  and  the  acoustic  noise  has  had  time  to  i^rojjtigatc  acioss  the 
anode-cathode  separation,  the  linear  regime  of  tlie  instability 
sets  in,  whereby  the  primary  electromagnetic  v.-avo  feeds  the 
secondary  electromagnetic  (em)  wave  U}-"  and  the  acoustic  noise. 
Eventually  tlie  primary  cm  w'ave  is  totally  depleted.  Tlnis  results 
in  the  secondary  em  wave  also  dying  out  in  one  or  two  round  trips 
in  the  resonator  because  of  the  output  cout)ling.  However  Eq .  (167) 

implies  that  the  nonlinear  process  of  supplying  energy  back  to  the 
primary  wave  will  not  take  place  because  of  the  dissipative  nature 
of  the  acoustic  waves.  The  energy  can  be  m.ade  to  flow  into  the 
secondary  em  wave  and  the  acoustic  w'ave  but  will  not  flow  back  to 
the  pump.  Thus,  once  the  pump  is  depleted  the  instability  stops. 
The  entire  process  begins  again,  howev'cr,  because  ug  builds  up 
due  to  the  laser  properties  of  the  medium.  The  acoustic  v.’avo  may 
rcm">in  into  this  region  because  of  its  relatively  long  decay  time. 
This  feature  may  load  to  a  degradation  of  beam  quality  of  pulsed 
lasers,  as  well. 

D.  SUMMARY 

The  physical  processes  involved  in  the  instability  are 
illustrated  in  Figure  9. 

Consider  a  large-amplitude  w'avo  at  frequency  o  ^ind  wave 
number  kg.  In  the  presence  of  a  density  fluctuation  at  frequency 
u)  and  wave  number  k,  a  second  el cctromagnct  j  c  wave  (  .'g  -  , 

kg  -  k)  is  excited  through  tlic  nonlinear  polarization.  This  wa'/o 
is  a  higher-order  mode  of  the  unstable  resonator  and  iicnce  could 
have  a  higher  loss  rate  than  the  fundamental  mode.  Therefore,  it 
may  not  bo  excited  in  the  absence  of  the  density  fluctuation.  Tlie 
two  electromagnetic  wav^cs  in  turn  amplify  the  density  fluctuation 
through  the  heating  term.  In  addition  they  excite  fluctuations 
in  the  equilibrium  gain  throinh  tlio  flux-saturation  term.  Tli  i  s 
completes  the  feedback  "looji."  .nit'.ce  the  h-'at  depen'.  i  I  i  on  is  a 
dissi[jative  process,  this  instability  is  Tioncoi'isc' rva  t  i  \'e  and.  turns 
itself  off.  I’ho  reason  for  tliis  is  that  once  .all  tho  energy  from 
the  primary  wave  (i‘'gr  k.Q)  has  been  deposited  into  the'  seconrlary 
clectromarjnct  i  c  wave  and  th(?  acoustic  wave-,  lliere  is  no  v;ay  of 
returning  the  eiicrrty  to  the  prim<ary  wavi'  as.  would  be  ihe  c.is.i'  of 
ordinary  conservat  i\'e  nonlinear  inst.ibilit  ies.  H('cau.se  of  t  lie 
additional  loss  of  the  secondary  wave  (highr-r-order  mock-)  it,  too. 
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Figure  9.  Physical  Processes  Involved  in  a  Three-Wave  InstabiliL 
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leaks  out  once  the  primary  wave  has  decayed.  The  whole  process 
restarts  because  the  primary  w'ave  is  once  atjain  built  up  by  the 
lasinq  medium.  The  acoustic  v;avos  may  remain  present  throufhiout 
this  phase  because  of  their  relatively  loncj  dam^^iing  time  constant. 

The  equations  we  have  derived  in  our  differential  formulation 
reduce  to  those  of  Section  II  when  the  two  elcctromaynctic  modes 
are  degenerate  in  frequency,  as  they  must. 


The  process  of  dielectr 
shown  to  be  unstable.  This 
including  the  effects  of  ga 


ic  modulation  was  studied  and  also 
is  a  form  of  Brillouin  scatterirg 
in-flux  oscillations. 
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SECTION  V 


UNSTABLE  RESONATORS  'NITH  SATUllABLE  GAIN 


A.  INTRODUCTION 

Unstable  rosonatorr;  liavc  fonnd  widesinaMd  vim  in  lii  n!.- 
gain  lasers.  As  early  as  19&[)^''')  it  was  la'ce.in  i  zed.  t;,iV  t  liLsa 
resonators  arc  desirable  for  a  number  of  reasons.  Tlit-se 
include : 


1)  The  large  mode  volume  that  can  be  obtained. 

2)  The  uniformity  of  illumination  {and  h.ence  enemy 
extraction)  soon  by  the  gain  medium. 

3)  The  relative  freedom  (compared  to  stable  resonators) 
from  undcsired  liiglier-order  transverse  modes. 


In  thc'  realm  of  geometric  optics,  the  oiie'ration  of  an 
unstable  resonator  such  as  the  one  shown  in  Finure  10  is 
obvious:  The  output  is  a  plane  wave.  However,  whenv  diffrac¬ 

tive  effects  arc  included,  mode  properties  are  not  ol'vions  vind 
the  calculation  is  complicated  in  the  extreme.  The  uncL- rst  a;id- 
ing  of  unstable  resonator  properties  that  lias  evolved  ovew  tlr'’ 
past  15  years  has  been  derived  from  the  -rradual  devc- 1  op.men t  of 
ever  more  sophi  £;t  icated  calculation  teclin  iepaes  .  ^  ^  ^ )  Prior 

to  this  report,  thc  study  of  unstaljlc  i.'csonators  had  advancc‘d 


9  .  S iegman  ,  A .  E .  ,  Proc.  IEEE  5_3  ,  p .  277  (1  065). 

10.  Fox,  A.,  and  Li,  T.,  Bell  .Syst  .  pp .  453-488. 

11.  Stroiffer,  W.  ,  IEEE  J.  Quant.  El-'Ct.  OE-4  ,  pp.  220-230 
(April  1968). 

12.  Bergstcin,  L.,  Appl  .  Opt  .  1_,  p.  495  (1968). 

13.  Chen,  L.,  and  Felsen,  T,.,  IEEE  J.  Quant.  Elect.  QE-9, 
pp.  1102-1113  (Novemiitr  1  973). 

14.  Hon-.-i  t7. ,  P.,  J.  Onl  .  Eci.  Am.  63,  inv.  1  528-1  54  3  (December 
1  973)  . 

15.  Moore,  G.,  and  McCarthy,  R.,  H A_m_.  Gy,  pp.  2  2  8- 
241  (February  1977). 
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rc  10.  Tcloscopic  (.Jnf.tablo  Resonator,  Showing 
Col lima Led  Output  Beam 
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to  tho  point  where  the  effects  of  diffraction  and  dotcnninistic 
intra-cavity  qain  distributions  were  under;;  tood  .( 1  )  In  thi;; 
reijort  v/c  carry  tho  lc;vel  of  understandinci  one  step  furth.or, 
and  include  the  self-consistency  problem  associated  with  gain 
saturation  as  found  in  actual  lasers. 

B.  NOMENCLATURE  AND  STATEMENT  OF  PROBLEM 

In  this  section,  we  express  tho  coupled  diffraction  and 
saturated  gain  problem  in  terms  of  the  ono-dimonsional  strip 
confocal  unstable  resonator,  as  shown  in  Figure  11.  Coninarablc 
analysis  can  be  made  using  two-dimensional  rectangular  symmerry 
or  circitlar  symmetry’  resonators,  but  wo  wi;;h  to  avoid  those 
additional  mathematical  and  computational  complexities  h.ere. 

It  should  be  noted  that,  for  the  case  of  empty  (no  ge.in)  reso¬ 
nators,  tho  diffractive  properties  differ  only  in  detail  among 
the  several  geometries. 

In  the  absence  of  either  a  saturable  or  non-uniform  gain, 
the  optical  properties  of  the  strip  resonator  arc  described  by 
the  well-known  integral  equation (16) 


Xf  (X) 


_ /  2  i 

^  V  -  1 


eff 


+ 


+  e 


c 


exp 


cf  f 


.  f{x’)  dx' 


(168) 


In  this  equation 


f (x)  =  the  complex  optical  amplitude  in  tho  plane  of  the 
feedback  mirror 


X 

M 


F 


eff 


tho  eigenvalue  corresponding  to  f(.x) 
geometric  magnification 

the  so-called  effective  Fresnel  number 
(a^/X’L) ( (M  -  1) /21 


X'  =  wavelength  of  light  (not  to  lu  confused  wi  t  li  ) 


16.  Weiner,  M.,  Arnd  .  Oet  .  18,  pp.  1  8  ^  x- 1  h  (Uune  I'^'V'’).  Tlu' 

basic  rc.sonator  t-rp:  it  ioii  is.  d.-sciilid  in  nuo' a '  "is  i.ipers, 
but  thi;;  work  i  tls'  only  on<.  Ui  less  Jt  I'Xplicitly  in 
terms  of  the  effective  Fresivl  nusi.'/cr. 
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2a  -  foodback  mirror  size 

L  =  the  normalized  optical  axis  offset  or  tilt.  r 
takes  on  values  from  0  to  1 . 

Equation  (168)  is  just  a  formulation  of  the  requirement 
that  the  optical  iimplitude  distribution  reproduce  itsi'lf,  to 
within  the  constant  1,  after  makinq  a  round  trip  back  and  forth 
in  the  resonator.  For  example,  in  the  limit  Ff.ff  difirac- 

tion  effects  vanish  and  Eq.  (168)  takes  on  the  form 

Xf(x)  =  f(x/M)  (169) 


wliich  has  as  a  solution  f(x)  =  constant  and  \  =  1. 
readily  identified  with  the  geometric  optics  result 
wave  scaled  in  size  by  a  factor  M  every  round  trip, 
output  coupling  seen  to  be 


This  is 
of  a  plane 
with  the 


=  1  -  (1/M)  (170) 

which  is  just  the  fraction  of  the  mode  not  blocked  by  the  feed¬ 
back  mirror. 

In  many  cases  of  practical  interest,  liowovor  (e.a.,  infra¬ 
red  lasers),  a  qcomotric  description  is  not  sufficient  and  dif¬ 
fraction  effects  must  bo  included,  even  though  the  optical 
properties  may  bo  largely  aeometric  in  character.  This  corre¬ 
sponds  to  solutions  of  the  integral  equation  for  which 


1 


eff 


100 


1 . 5  <  M  r  4 


For  a  given  combination  of  M,  Fpff  and  i  ,  t  tuu’e  are  several 
solutions,  each  with  its  eiuenvalue  and  ci uenf unct i on ,  corre¬ 
sponding  to  the  several  transverse  mochrs  of  t  1k'  resonator. 
With  tliG  eigenvalue  equation  defined  as.  ’ tv;.  (  1 1  fi )  ,  tdie 
(one-dimensional)  output  coupling  l,^;  ( i  )  ar.scK- i  a  t  twi  urith  each 
mode  is  given  by 


L,(i)  -  1 


(171  ) 


Wt;  may  as 
output  coii[jling 


CSV.  the  effC'Ct.  of  d  i  f  f  r  a  e  t  j  ,  ,  1  ,y 

or  equivalently,  pi,  as  F,,.)  J 


nr i  1  ■  ■  :a  ng 
\M  r  i  ■  'ci . 
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Figure  3,  fro;n  Ref.  13,  shows  such  a  curve  for  .M  -  2.9,  ;  =  0, 

and  Feff  boLwcon  5  and  20.  An  overall  per  ioflic  i  t  y  :n  F,.ff  ir, 
apparent,  as  v.'oll  as  dc'qeneratc  crosslnq  points,  (e.q.,  Fj.r.-  cr  R.9) 
where  two  niodc-s  share  the  same  value  of  |'  l-  Iti  an  actual  las.'r, 
such  points  may  cause  uridcsirable  mode  hop.pinn.  A  very  simila.r 
pattern  of  diffraction-induced  output  couplinci  variation  c::ists 
if  F(2ff  is  kept  constant,  but  the  position  of  tiie  or-tical  axis 
(c)  is  var  ied .  {  i  7  ,  1  8 )  Even  if  the  mode'  arc  not  .ictiially  dLM:>;-n- 
crate,  the  magnitude  of  the  output  couplint;  ripplt'  i;:  importai'.t 
for  assessinu  mode-medium  interaction  cff^'C'ts.  'Itnis  latter  s-iti;- 
ation  is  perhaps  the  more  likely  in  an  orieratiiui  device,  as  it 
can  readily  l^e  induced  liy ,  say,  a  mirror  vibration,  or  by  a 
density  perturbation,  leading  to  MMI  instabilities. 

In  any  realistic  high-efficiency  laser  tfio  gain  will  alvays 
bo  saturated.  The  effect  of  the  presences  of  ciain  s.at  ur.it  ion  on. 
theoretical  curves  such  as  depicted  in  Fiauro  12  h.is.  n.en:  tc;  o:'. 
not  been  assessed,  and  provides  the  purp'oso  for  tiiir.  nt.niy.  In 
the  main  theoretical  analysis,  we  consider  only  ali  r'c  n  rui- 

tors  (c  =  0)  ,  but  in  .Section  E  we  si'.ow  liow  tilted  r..';v  .r;  a  t  c  r 
behavior  can  be  related  to  the  aligned  case. 

C .  METHODOLOGY 

For  the  present  analysis  v.'e  l-avo  developed  a  w  A.  n  I  •  :  , 

combining  a  num.ber  of  tochinicjues  wliicli  ii.ivi,'-  not  : 

before.  The  result  is  an  iterative  procedn.re  •..naio:-;  -nl 

not  only  the  lowest  loss  mode,  but  also  i  he  h.i  vh.er-:':  :•  n 
tor  modes,  simultaneously.  Even  for  a  hiclily  sa t  ■ ; inC  i,;  n,.  ii  nn, 
convergence  is  obtained  in  five  iterations  or  I'-'ss,  a:'.,!  t!,.  r.^  i  ■ 
no  difficulty  obtaining  a  solution  in  tlic  presence  of  m.eae  dt  s. 
eracy . 


We  outline  here  the  t  oclin  i  ejues  whicli  c:ompri:se  t':...  c'-re.r- 
nents  of  the  algorithm.  Howcvc^r,  since  tl,.'  th.i.'ort  -  icil  as:', 
of  these  arc  in  places  quite  intricate,  th"  iS'  sK  r  i rr'ft  tua  1 
to  the  original  references  for  complett;  details.. 

1.  Initial _ Empty  Resonator  Calcul  iti'On 

For  a  given  i'lot  of  resonator  via  rar,'.'’ t  e  r  s.  >c  ,  with 

as.sunied  zero  for  the  ics.t  of  tlie  d  i  scuss,  i .  .n )  ,  t  >..  •  is.iv'y  r.'scriatcr 
eigenvalues  '' i  and  cor  res.pond  i  nu  mode  d;  i  :•  t  r  i  lu  t.  i  c  v, :  :  ^  (;■■,)  isi' 

detc'rmincHi .  The  a.symiUotic  thf'ory  of  lesamat  i.e:;  js  us.ivi.  (11) 


17.  I’erkitis,  J.  ,  and  Cason,  C.,  Appl  .  ld>y  1,' ■  t  t  .  H  ,  ]  p-  I'bi- 
200  (Au<p;s,t  1  977  )  . 

IB.  K'^'lli'ii,  r.,  .'ind  Fi'.iith,  M.,  Ov>t  .  Env.  18,  ;.p.  I'i7--lf0 

(March/Apr  i.  1  1  979)  .  '  " 
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This  theory  can  be  suiiraiarizod  ,  in  fairly  uliyr,  lea  I  t  e  r;:'.:;  ,  ^  ^ ^ 
as  follows.  The  amplit\ide  function  fi  (x)  is  t  ;-; 


series  of  slowly  modulated  cylinder  waves,  plu;: 
as  follows: 


.ein'U/L;  a 

a  core  term 


f^(x)  = 


N 

E 


i  k  r^^(x) 


i  k  r^^(-x) 


f  ■  (x) 
n  ,  1 


+  f  .  (-X) 

-  n,  1 


-r  (->;) 
n 


(172) 


The  origins  of  the  cylinder  waves  arc  successivclv’  furtlier  rcsnc'/ed 
virtvial  images  of  the  feedback  mirror  edacs  (Finuro  13)  ,  corre¬ 
sponding  to  light  that  lias  made  n  trips  around  the  resonator. 


Then  is  the  distance  from  tlio  n^h  virtual  imano  to  t!ic  obser- 


n  1 


feedback  mirror 
(172)  becomes 

plane.  In  the 

exp 

2  ik(x  -  x^^) 

w 

\  n  / 

n 

n ,  1 


(X) 


(173) 


+  C;xp 


r  1  21 

2  ik(x  +  x,^)2 

f,(-x) 

W' 

n 

/*s 

+  (X) 


The  origins  of  the  cylinder  waves  are  at  the  points 


.2n 


w  =  d  (M  -  1 ) 
n 


Xn  =  aH 


( 1 7  i  a  ) 

(171b) 


Use  of  Fq  .  (173)  in  Fq .  ( Hi  8 )  loads  to  a  series  of  infenrals  wi-.ie 

arc  solved  l^y  tlie  n’.otiiod  of  stionai'v  pliare.  (H')  Tiiis  is  v.ilui  for 


1  1  (:  s  1  C'  r 


Ffifj-  > '■  1,  but  ha;;,  bio'^n  s.  iu'wn  to  uiv'  reasonable  las 
keff  ~  !•  Wlion  thi;:  is  dcuie,  tlie  modulatirui  fimctio 

the  corr'  term  i  ^  are  four.d  tc.i  10"“  expreuss  ili  1  e  in  t  err; 
resonator  na  r  a.r.i' a,  , ;  rs  ,  Tile  final  result  i;;  a  crargle;-:  polynomi 
in  a  param-rter  1:  : 


•  n  ,  1  '  ■  ■  ■ 
r.  f  ('  t  bi.  '  1' 


19.  [lorn,  M .  ,  and  Wolf,  F.,  Prinoiejes  <>;  0:>;  les,,  fgh  j.d. 
( Pe  rita.mon  ,  Pew  York,  l')7a),  pp'.  '!  \>2-'i . 
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Figure  13.  First  Virtual  Imago  of  Fec'dback  Mirror  Fdaos. 

Coherent  addition  of  contrihution.a  f rora  numerous 
such  virtual  images  oives  rise  to  the  complicated 
mode  structure  of  the  resonator. 
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1  =  - 


2^^  \  2  F 


N  +  1 


•  E 


n  =  1 


exp  2  TT  i  F  ..  B 


ef  f  ‘  n  /  ’  n 


+  exp  I  2  r  1  F 


-1  V.  +  1/2  I  ,n 


ef f  ' n  /  ■ n 


for  symmetric  modes  v.’hcn  n  =  N  +  1 


(175) 


The  plus  and  minus  signs  denote  symmetric  aiid  antisymmetric 
modes.  Here 


e„  = 


1  -  M 
1  +  m' 


(176) 


and  the  original  eigenvalue  X  is  found  from 


X  = 


(177) 


The  lowest-loss  mode  corresponds  to  tlio  larges  value  of  j  X 
The  corresponding  eigenfunction  is  found  from  using 


f  .  (x) 
n ,  1 


=  (VM)'^(I  -  M 

(1  -  (a  +  m'^^x) 


(178) 


^  ^  ^  -0—1 
fi(x)  =  [ad  -  M  ^)] 


2  miF 


(179) 


in  Eq .  (172)  .  Since  the  expression  for  f  lias  a  facLoi-  1  -  n  in 

the  denominator,  and  since  |xj  =  ]  / 1  ii  [  ~  1  for  tlie  ]  ov.'e:;t- less 
mode,  it  dominates  Eq .  (172)  for  regions  in'md'  tlie  r.hadow  boun¬ 

daries  (x  -  +  M<i)  and  represents  tlie  "geometric"  part  of  the 
solution , 
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The  number  of  virtual  imaqcr. ,  N,  v.’hich  are  ncK''de(.1  to  fix 
the  order  of  the  polynomial,  can  be  e;;tiniated  in  the  follov/inq 
way.  (20)  Consid<->r  t’ne  confocal  unr.tai^lo  rouonator  shov.’n  in 
Figure  14  v/ith  a  sphei'ical  wavefront  emanating  from  the  common 
focal  point.  Radiation  within  the  angle 


0 


1 


D/M 


(180) 


which  makoc  one  round  trip  tlirouqh  the  cavity,  will  exit  magni¬ 
fied  to  the  whole  mode  size  D.  Similarly,  radiation  within  the 
small  core  angle 


0 


n 


(181) 


which  makes  n  round  trips,  is  magnified  M  tinvas,  and  also  exits 
with  the  full  mode  size  D.  [This  core  region  is  v/hcre  (x)  in 
Eq.  (172)  arises.]  If  we  desire  the  cavity  to  bo  "diffraction 
limited"  in  the  usual  sense,  the  exiting  radiation  must  have  a 
cavity  angular  spread  of 


0 

c 


(182) 


By  requiring  the  angular  spread  of  the  core  region  to  match  that 
of  the  whole  cavity,  i.e.,  0^^  =  0^,  wo  have 


whence 


D 


M^f. 


X 

D 


(183) 


N 


In  (8  M 
In  M 


(184) 


In  practice,  it  turns  out  to  bo  useful  to  add  two  or  thrr’o  to 
N  as  calculated  by  Kq.  (184),  to  guarantf'c  th.it  N  will  be  larcie 
enougli.  A  similar  formula  used  by  Uorwitz(l'^)  in 


N 


In  (250 
In  M 


(185) 


20.  Anan'ev,  Y.,  Sov.  J.  Quant.  Elect.  5,  j.p.  615-0  17  (1075). 
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2. 


Intensity  Ronorma 1 i zati o n 


The  lov.’cst- loss  mode  intensitv  distr  i  livit  ion 


fo(x,d)  I  2 


,  corresponciina  to  >q,  is  norma  1  i  x.cni  to 
ated  intensity,  Ig.  To  do  this,  wo  co!',r.'  to  a  u<'o; 
intensity,  Ig,  usincj  the  theory  of  qeomo'tric  ;;tri; 
with  saturable  aain.  (21)  v.'e  then  require  that  i.‘i. 


I  ( )  - 
t'ne  satur- 
le  trie 

5  innajiK! t  ors 
Jtvexnn-.'c 


intensity  over  the  feedback  mirror  be 
out  diffraction: 

+a 


the  same  with  and  with- 


-a 


l(x)  dx 


(186) 


This  is  to  a  certain  extent  heuristic  and  is  strictly  true  only 
in  the  limit  F^ff  “b  but  is  a  rosoniiblo  requirement  for  the 
range  of  Faff's  of  interest , 

In  the  geometric  theory,,  we  begin  with  a.n  afu'ur.m'd  fi'-ld 
form  in  the  cavity;  namely,  a  modulated  cylindrical  wave  ema¬ 
nating  from  the  confocal  point,  and  a  modulatc.-d  plane  wave 
exiting  the  laser: 


E 


ik  ( z  +  L  +  d) 
e 


f (z,r) 


g  ( r  ,  '* ) 


(187) 


(Nomenclature  is  defined  in  Figure  15.  The.'  funct.ion  f(z,r)  hero 
is  similar  to  but  not  the  same  as  the  f  function  oi  F'.  ction  1.) 
The  boundary  conditions  are  that  the  amplitudes  vanish  oil  tl'.e 
mirrors : 


f (0,r) 


1 


(L  +  d) 
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L  +  d. 


_ 

L  + 


and 


f (d,r) 


1 

1/2 
d  ' 


c 


2ikL 


By  starting  with  the  time-dcfjcndent  wave  x'quation 


c  V  E  = 


fit 


2 


2  c  G 


t 


21  . 


Monro,  C , ,  and  McCarthy, 
227  (February  1977). 
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Figure  15.  Noineucl  aturo  for  GeomoLric  Strip  RosoMator 
(after  Ref.  21) 
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and  neglecting  all  second-order  derivatives,  wc  obtain,  using 
the  assumed  form  for  E,  the  rate  equations 


Using  Eqs ,  (192,  (193),  and  (199)  in  Uq .  (191)  yields 


dH 

G 

0 

dz  1 

2H  ^  1 

f  L  4  d 

\  -2H 

^2 

1  + 

0  +  I 

^d  t  L  -  z 

)‘'  J 

f 

g 

(19 


Determining  a  solution  consists  of  findin<!  a.  value  of  such 
that  the  boundary  conditions  [Eq.  (IBH)!  ere  obtai;u.'Cl  vdien 
Eq .  (196)  is  inteqratod.  Then  Iq  =  i  M  f '  j  j  .s  used  in  Eq .  (186) 

to  normalize  the  diffractive  result. 

3.  Propagation  Througli  the  Gain _ Ee^gjon 


The  renormalized  fo(>-)  is  propanaiod  U'.rcnich  the  cjain 
region,  using  simultaneous  forward  and  ijach.v.'a  rd  propauat  ion  .  ^  ^ 
This  process  is  illustrated  in  Fitnii-o  16.  Ccuisidcr  a  point 
(>:,w)  in  the  roaion  bctv.’oon  mirrors.  Here,  v;  is  a  coordinate 
measured  from  rig'it  to  left  av.’ay  from  th.e  feedL'>ack  mirror. 


The  "  f  on'.’ard "  part  of  tlie  pror'aepat  ion.  consists  of  th.e 
expanding  protion  of  fo(:-0  reflected  into  the  cavity  l.vy  tlie 
feedback  lairror,  and  t'nus,  propaea.tc;;  ‘‘rum  r'icht  to  left  in  ttr' 
figure.  In  the.’  abscnc<>  of  ejain  we  iiavo ,  from  the  scaliir  tlicory 
of  diffraction,  the'  com[.)K'X  amplitude  at  a  point  (x,w): 


^cx 


K(x,x^)  fy(>:2)  ^xp 


ik.  2 

2d  ■■'1 


(19 


The  exponent  i  .a  1  term  i  ;;  tli^  j’h.ise  curviture  resultin':  from  th.e 
curved  fcedba.ck  mii-ror.  in  tl!e  Eres.nol  app.roxi  mat  ion  K  is  tlu’ 
propagation  fujiction  defined  );y 


K(x, 


i 

ib. 


1/2 


exp 


i  k. 
2], 


(19 


The  "Ijackv.’ard "  perl  ion  i  r:  th.e  0:11  i  re  <  xtent  (C'  fr,  (:d  , 
prior  to  liittin-i  tl.e  miri’or,  evolvin':  I  ickv.M  rci  in  1  1 :  .  I'his 
propariation  is  effected  using  tlie  coiiqil  ■ 'x -con  j  uq.i  t  c  ’  proi.Mca  t  i  f  n 
f  unct i on : 


24  . 


I, cm  i  se  11,  W  .  ,  e  t_  a  1  .  , 
August ) 
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•4-o> 

‘  ^ 

—  W 

K*  (x,x^) 

fo<"l> 

dx  , 

(199) 

Thus  it,  too,  propagates  from 

right  to 

left. 

In  this  way, 

bot)i 

waves  move  into  the  gain  region  in  the  same  spatial  sense. 

Locally,  the  saturated  gain  is  computed  using,  again,  the 
Rigrod  formula: 


G  ( X ,  w )  = 


o 


(200) 


Of  course,  light  reaching  a  particular  (x,w')  must  have.’  t  I'aversed. 
gain  in  tlie  roaion  w  <  v.’'.  To  account  for  this,  we  insert  in 
Eq.  (198)  a  line  integral  term 


Ggx^^'^"'-1^  =  exp 


G(C 


(201  ) 


where  C.  is  a  vector  connectinci  tlie  obsfu'ving  point  (x,w)  witli 
each  point  (xi,0)  on  tlio  fccdbaclc  mirror. 

For  the  backward  wave,  gain  becomes  attenuation,  so  Diat 


^bkw 


(x,w,x^) 


=  exp 


G(.') 


(202) 


is  used  in  Eq .  (199).  This  scliome  permits  s.ystcin.it  ic  c.ilcula- 

tion  of  G(x,w)  in  the  direction  of  inerL-asing  w. 

4.  Modes  in  the  Presence  of  Gain 

The  analysis  of  Monie  f,  MeCartliy  ^  ^  ^  ac  ts  asyr,r.)tot  i  c 

calculation  of  unstai.ile  rc'sonator  mod.i's  •  .a  in  fui-'.ctit'u 

can  be  .sfiecified.  In  tlieir  work,  G(>;,wd  v;a,;  a  [h  ir  t  u  1 , 1 1  ••d,  f-anc- 

tion,  Init  this  need  not  b’-  so  in  uianoral  ,  aani ,  as.  r.liown  a!;o\'i', 
we  calculiiLe  it  from  tlie  m<Hh'  ]'atl.ern  itself. 

To  include  the  C'frefds  of  qain  in  the  law.cnni  1  oi;  ar.e.lys.i:' 
the  eigenvalue  polynomial  [  Eq .  (17'i))  is  nujdifi'd  t')  I'ead 
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+  r'  (-a)  exp  (2  x  i  F  p 
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■(r-lxi 


(20  3) 


for  symmetric  modes  when  n  =  M  +  1 


The  nev;  functions,  rj!,(x),  are  integrated  gain  functioi'is  verv 
similar  to  Eqs .  (201)  and  (202)  ,  except  tliat  tlie  ’.'uth  of  inte¬ 
gration  forms  n  round  tri{:)s  tlirougli  tlie  cavity,  startiii  !  at  tlic' 
mirror  edge  (x  =  a)  and  terminating  at  tine  point  x.  Sp’.'c  1  f  i  ca  I  1 , 


r;(x)  ^  exp 


p 

E 

r.'  -  1 


G[X^(v;'),w-]  dx- 


o 


(20-5) 
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G[X2(w'),w'1  dw' 


In  J  G(0,w' )dv.-'| 
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where  the  integration  paVh.s  are  alone;  rays  with 


X3^(2') 
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In  this  analysis,  wo  always  assumo  that  the  1  owost- 1  os,;-,  riffle 
dctcrntines-  the  intensity  in  the  qain  region,  so  that  T'  func¬ 
tions  have  an  implicit  i  =  0  subscx'ipt. 

^  •  Iteration  to  Convo rciencc 

The  updated  f  (x)  amplitude  distribution  is  rirfq.  '.fnw  oj  into 
the  gain  region  as  in  step  3,  and  a  new  f(x)  calculated  as  in 
stop  4.  This  process  is  repeated  until  conveT'ceriC"  is  ote^  a  i  n.(  d  , 
by  which  we  mean  that  the  change  in  the  lowest-lcnn;  eigenva  1  xu • , 
Xq,  is  less  than  some  values 


<  6  (20G) 

Since  \q  ~  1 ,  this  is  the  same  as  having  the  fractional  change 
into  Xq  bo  less  than  li . 

For  a  GqL  of  5  (i.e.,  round  trip  unsaturated  intcn.sity 
gain  is  ,  a  convcracncc'  criterion  of  '  =  0.002  rc>(.]uires 

~  5  iterations;  for  =  0.0005,  the  nunilier  is  ciclit  tet  ten. 

There  appear  to  be  two  reasons  why  convergence  is  rapid  '••.•ith. 
this  algorithm.  First,  tl.o  saturated  ce.in  ca  Irul  at  i  on  always 
heains  with  an  dcriv''’d  from  ari  cm:  ty  rojw.'nator  calcula.- 

tion.  Thus  in  a  si'iist  it  tr'sats  nain  a;,;  a  C'crturl  .;t  i  on .  A 
second  reason  is  that  t);-'  intenratod  gain  functiot'.s,  1'  {:■:),  are 
fairly  smooth,  I'ecausc  t  I'.ey  intecrale  out  fine'  :aca]e  1  one  1 1  ud  i  :‘'.a  1 
(w-di  rot.:t  ion)  gain  ripnl'.,  and  do  not  chance  Ixy  much  betva'-en 
iterations . 


It  mur;t  ho  emt'ha.si.u  1  that  th.e  iterative  a  1  nor  it  h.;::  ■..".•e- 
sonted  in  t.h.is  papi-r  is.  '  u:'.damcnt  a  1  1  y  d  i  f  f  • f  ro:n  tie  Ft 
itcrat.  i  on  ( 1  c.-ne;'ally  to  .simulate  l.'un.'r.s  wild:  sa.te.r:'.:'.! 

gain.  in  the'  Fox-J.i  ,  thu  rr'r'-ona  t  tu'  i  :it  .;r  ■;  1  1 1 ;  n 

itself  is  iterat^-cl  to  fn.xe.’eT'tn.nico ,  .'‘.tart  inn  v.'ith  com''  ^  u  ri 'U 
inititil  mode  rdiape  .sucl:  .la.  a  pi. me  Vvavc'  or  C'>ai.\:;s  i  a;i .  In  s',;ch 
schenu'f;  ,  the  only  k:K3v;l(d  c'f  tinv  mocio  s.u rv i  i  nc  '  p.iss  tt' 

pass  is  contiinf'd  i  rx  t;v'  nu.merical  array  whlcli  'c  i  f  i  ■  - f, 

The  inath.rmiat  ica  ]  man  i  pu  1  ,i !.  i  on  nc't  ded  to  p,rei:  no  1 1  •  ■  f,,!:-:)  .iround 
the  ctivity  of  nocc:s,'-.  i  t  nccunralati  s  c'lux'r,  (s;!  h  fro:;:  nn-'!.  ;uv,;;‘c.  . 
as  roLincloff  error,  aixd  f  roni  tir’pro;-:  i  n;  a ;  ions  '..iso  i  is.  ore: 

tion  schem.'.  Tho  necessity  for  m  i :  ti  t  a  mi  i  n  ■:  a  hichly  iccuintf 
copy  of  f  (:■:]  implies  hich  .s,imp !  j  f?-e<u;'-ncy  ct.-;p.;;ad  ;  >  t-.-p;c.i] 
structure  scale  sire's  coi.miincd  in  f(x),  with  tie'  i.ssi;  thm. 
comrniter  slcainu.'o  g:u  i  le  anon  *_  s>  diotate  a  pa'.ictic.!!  up:  <-  ]  nu  t.  mi 
f'c'ff»  visually  around  unaty  ajid  certainly  li_s:;  t  h.ux  ID. 


25.  .Cn  i 1  ,10  ,  F.,  and  Fir-cmui,  A.,  Aptil  .  O:  .  It,  i  i;.- lo  i 

(Au'iunt  1''7").  'I'ii’e  oomrut.-r  •:‘imd  in  thin 

})a;ier  in  t  Fie  nrchc-tyiM'  of  a  immbor  ul  s.  i  mu  1  a  t  1 1 ;  n  des  in 
use  t-otlay. 
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In  contrast,  tlic  method  presented  hiere  .nec'd  only  yarn  the 
ciqcnvaluc  Xq  from  iteration  to  iteration.  Fur  tlv'rr.iore ,  r;inci.' 
the  r'  (x)  functions  are  .sen.nitivc  only  to  tl;c  latcjent  n.patia] 
vnriatio’  s  in  gain,  a  highly  accurate  propagation  through  the 
gain  region  is  not  essential. 

D.  RESULTS 

The  basic  purpose  of  this  calculation  is  to  asr.ess  the 
effect  of  cfain  saturation  on  well-known  diffractive  mcjde  prop¬ 
erties.  To  this  end,  v/e  consider  the  particularly  il  lustrati\’o 
case  of  a  resonator  with 

M  =  2. 9 

F  =  8.892,  9.390,  and  9.803 
ef  f 

Referring  again  to  Figure  12  those  F^,ff  values  correspond  to 
(in  Horwitz's  terminology)  a  so-called  crossing  point,  a  point 
of  maximum  mode  separation,  and  a  cusping  point,  respectively. 

To  choose  approi^riatc  small -signal  gains,  wo  consider  tlv? 
steady-state  oscillation  condition  for  a  strip  resonator  in  the 
geometric  limit: 


1  2  G^L/(1  +  (I/I3)] 

M  °  1 


GqL  =  ^  [1  +  (I/Ig) ]  In  M 


(207) 


(208) 


We  will  consider  the  three  regimes  of  interest: 

I/Is  «  1 

3/1,  -  1 
I/I^  >>  1 

A  very  low  value  of  l/l>.  is  tlm  same  as  th.i,'  cxnpl  y-re;.(uiator  cal 
culation  and  i r,)]cu]ated  as  a  natter  eif  r'ourS''  lo’"  any  cit 
case.  As  an  example  of  tlio  interngdi  at  <a  ca:;e  we  tain-  Gy,!. 

In  ^  l.OGf,  so  that  1/1.-.  =  1.  For  3/1,-  >  1,  v;o  liave  chonaui 

GqT,  -5.0. 
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A  co:i!;  ).:i  r  i  ;;o:i  of  prf5o;.^  nirt:.c>l  v;itAi  otin-r  i 

frustrated  ty  th"  y^aunty  of  otS':r  ri'-sulfs,  a;ui  ry  tls-  fact  ti.it 
since  other  si'e.ha;?!,-;  ft. 1  :t  •.■ration  rcv:uirin(:  ]  o-,v  V,..'.-,  n! 

results  are  available'  at  Live  cou'ditions  listed  at  the  leC  a:  nn  i  n 
of  this  section. 

Nonet  ho  1 '.!ss  ,  ’.ve  havt'  ca  1  cu  1 t.eJ  a  car-.c-  foi'  coripa  r  i  sc 'n  , 
with  results  obt.iiiiiHi  by  K'.'n-.-.cl;  atal  Ch.est'  (26)  ] »,  t  ir  '.■.'c.  rh.  , 

a  Fox-Li  strip  le^r-onator  calculation  w.i.';  ir.c.^d  i  f  i  eci  to  i  n.c  1  .'.dt  ■  a 
satiH'able  oain  slu'ot  positioned  rifuit  at  the  eririarv  :nir!ej;  : 
the  amplitude  v.’ave  as  prop.iaato'l  to  tt'.e  p.ri.tr'.ry  ;airr'''r,  an.-pli- 
fied  and  saturated  usine  the  Rieinxi  fornula  (v;ith  local  inf-  n- 
sity  doublc'd  to  accov.nt  for  incotdnc:  and  outeoi:!':!  '.■.•■•.••■s)  acn.; 
propagated  back  to  tlio  feedbvick  lairror.  Ti;C  re.sonator  parari- 
eters  arc: 

M  =  2.5 


F 


Off 


0.64 


G  L  -  4.583 
o 


this  value  of  F^ff-  tlio  asymptotic  nrr-tlv  d  is.  bx'lcw  its  rr:ni;;-,.'' 
of  proven  validity.  Also,  Rensch  iinc;  Ci.  -stv.';'  dc  ra'd.  :  rc-.- i 
eigenvalues  but  only  intensity  d  i  s  t  r  i  ha;t  i  on.: .  ’bn.us  a  ccm.p.'ir  i  sct; 
of  tlu'ir  results  v;itli  ours  must  of  necessity  be  soas.-v.-iiat  srccu- 
lativc . 


Figure  1.7  compares  bare  resonator  i  n  t  r'ns  i  ty  patt.  rns,  a:-;:! 
gives  a  clue  as  to  hov;  th.e  asymptotic  nt.-tlicni  brr'ak-'  dev.'s,  a.t  vr'ry 
iov;  F(^ff;  the  fine  scale  detail  is  lost,  tliou  hi  Ih"  ■:rc.‘-"'  rh,.v.  •: 
i.s  basically  tho  same.  The  situation  v.’itl:  ''niin  is  ecru c  d:  !'■ 
(Figure  18)  .  It  will  be  noted  that  th.ein-  is  a  r-liei'.t  di;  ftu'.  ncr 
in  tlte  location  of  the  side  maxiraa.  .In  our  ca  1  cu  1  i  t  -p,  ^  t  h  ’p 

occur  at  slightly  hichor  value.*-,  of  x/a.  su rnu.-r-t  th't  tl’.'lr- 

is  due  to  an  imprec  i  s  i  oj',  on  Rc-nsch  and  Cl."Stor’s>  put  in  calc-r.- 
latinq  the  effective'  Frc'sne-l  numl'.' r  ,  ;;  i  nc-e*  witdi  *d.  i  r*  t  >  ■  in.i''- 

qration  sch.enie.s  there^  is  always  an  uncc'rtainty  in  just  ■chcix'  Ida.’ 
feedback  mirror  ends. 

For  thf'  ;>rincipil  c:a  1  cu  1 , 'it  i  rm  s  ,  ;  la  cain  v.-.!S  C'.::‘“i:r  d  t  . ' 

a  rec'ion  from  x  --  0.25  L  to  x  =  0.75  b,  si;:ce  in  r,  us.-  lu  ui  :  r 

devices  (wh''!'e  unsl.dib'  r*  soiia  ‘  '  r  s  i’ i  n  d  ids  or-’O'  u  !.  a:  ■  )  ike 
mirrors  are  si'':-i,-.v;ha  t  T'em-jve'd  fic'm  the  lasitu!  r.'-liusi.  M  *•■.■,■•■•.■  e ;  , 

thi.'-  i.s  not.  r'f,:u -iit  i  a  1  \'i  l  lie  a  1  nor  i  t  I'.n  ,  .ukI  t  h*,-  aa  i  n  ran  be 
plac'd  arbttrarily  b''tv.'een  the:  mirror;;. 


IMTENSITY  (ARBITRARY  UNITS) 
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Rosiilts  for  Ff.ff  =  8.892  and  GqL  =-*  1.065  cire  .shov.'n  in 
Figures  19  througli  24.  Figures  19  and  20  shov/  intensity  a;vd 
phase  for  the  first  pass,  i.o.,  before  ai'iy  uain  saturation, 
but  after  renormalization.  The  units  of  interisi'ey  are  s-jch 
theut  Icj  =  1.  The  results  after  5  p'asses  appotir  in  Ficyirc^s  21 
through  23.  The  first  depicts  r^(x)  for  veii'ious  n.  As  ri 
increases,  the  curves  aprjroacli  a  limiting  form.  As  or.a  v.’oiild 
expect,  the  have  maxima  and  minima  invein'ely  related  to 

those  scon  in  the  intensity  plot.  Inspection  of  th.e  bare  and 
saturated  intensity  profiles  shov;s  that  the  fine  scale  diffrac¬ 
tion  features  are  v'ory  nearly  tlna  same  in  both  cases,  Fuivia  tint 
the  overall  envelope  of  tlio  intensity  is  made  more  uniform. 
note  that  it  would  be  iiard  to  clai.m  t.hat  the  djffracti\’e  cem- 
ponent  of  the  pattern  is  exactly  suppressed  by  gain  saturation. 
Figure  15  depicts  the  saturated  gain  in  th.o  resonator.  Th.e 
optical  axis  coordinate  i.s  in  percontago  of  w/L;  i.e.,  0  is  the 
feedback  mirror  location  and  100  the  primary  mirror.  The  aain 
scale  is  in  units  of  gain  per  optical  axis  increment.  Thus  the 
unsaturated  gain  is  seen  to  be  1. 065/50  ~  0.02. 

Figures  25  through  28  arc  for  Ff.fr  =  8.89.2  Ivat  v.-ith  G.jh 
increased  to  5.0.  The  most  noticeable  effc'Ct  of  tlio  in.cr'  ac.od 
gain  is  that  the  very  small  amount  of  "stray  liul'.t"  well  beyond 
the  geometric  mode  size  becomes  signi  f  icar.t  ly  nm.r-lified.  Sin.cc' 
real  devices  often  have  scrat'or  aperture'?;  ;.:.ieoi  f  i  ra  1  1  y  lo  (•;•:- 
elude  such  radiation,  its  ap[)earanco  in  the  piv.T-ent  ca  1  cu  1  ?if  i  or-', 
is  probably  not  important.  One  can  nl;;o  discern  th.it  the  ranuo 
of  intensity  excursions  from  a  mean  value  has  !>ceii  still  furth.' 
flattened.  Intercstinaly  enough,  the  pliaso  profiles  are  nearly 
identical  among  the  three  cases. 

Figures  29  through  38  depict  rc.su]  ts  for  Fp^f  -■  9.39,  and 
those  for  F(.,ff  =  9.863  are  shown  in  ricyares  39  through  48.  It 
v;ill  bo  scon  that  the  same  qualitative  ol.)sorvat ions  can  be  made 
for  these  cases  as  well. 

The  prcccdincj  plots  do  not  hint  .at  thcr  dramatic  effect  of 
gain  saturation  on  cigonvalno,  and  Iv-nec',  output  coutdinc:.  ]n 
Figure  4  9  a  portion  of  the  empty  cavity  plot  of  Fieure  12  i 
reproduced  for  F^ff  betv/ocn  8.5  and  10.5,  alone  with  result:; 
from  the  pre.sont  calculations.  It  will  i)'-'  r.een  that  the  iriole 
degeneracy  at  ~  8.892  is  lifted,  even  for  tlie  valvn- 

of  GqI,  --  1.065  correspondi  nu  to  I  /  I ’  1.  Ca  1  cm  1  a  1 1  ,'.t 

closc’ly  r-gsaced  intcuval.s  ne.ar  Fp^e  8.8''2  Irav  ii  .t  the 

the  dc'generacy  is  truly  lifti'd,  and  net  jur;t  laove.:!  to  a  r.>  eiy 
Fresnel  nunrljor.  It  i.s  a  1  .so  cavidenl  that  even  at  1h,"  eus.pire; 
{Joint  th'-'  j'irr'ple  in  mode  r.ep.arat  ican  is.  redum;!,  1  I'.e'.;.  :ii  ,  suri 
singly,  tdu'i'c'  is  not  ar;  mucrh  (iiffi  renec-  1  e  |  u  th>'  (/g.,  l.tue. 

and  Gpl,  ■  5.0  cases  a;;  there  is  betwcian  IIkt-.c  sikI  tlu'  ''nii'ty- 
resonataar  cvis.e. 
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Figure  32.  Intensity  Profile  for  G^L  =  1.065. 
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Figure  49.  Effect  of  Saturable  Gain  on  Mode  Degeneracy.  The 

solid  curve  is  a  portion  of  Figure  12,  rof>rcsonL  ing 
an  "empty"  resonator  in  the  absence  of  gain  satvira- 
tion . 
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In  a  qualitative  sense,  saturable  gain  seems  to  produce 
an  effect  not  unlike  that  of  mirrors  with  rounded  cdqes^^^)  or 
tapered  ref ] ectivity ,  ( 28)  which  reduce  the  influence  of  diffrac¬ 
tion  on  mode  properties.  This  is  in  spite  of  the  fact  that,  as 
we  have  seen,  the  intensity  profiles  themselves  are  not  really 
smoothed  on  a  fine  scale. 


Horwitz  has  pointed  out  that  empty  resonator  mode  degen¬ 
eracy  occurs  for  Feff's  loss  that  some  critical  value,  with 
that  critical  Fgff  increasing  quite  rapidly  as  M  approaclics  1. 

In  other  words,  the  propensity  for  mode  degeneracy  is  emphasized 
for  low  Fcff  a^d  low  M.  We  hypothesize  that  for  each  valvie  of 
Fj2^f  and  M  corresponding  to  a  mode  degeneracy,  there  is  also  a 
critical  GqL  which  will  lift  that  degeneracy.  Thougl'i  wo  have 
not  explored  this  thesis  exhaustively,  we  have  performed  a  num¬ 
ber  of  further  calculations  using  G^L  =  2.  For  M  =  2.9  there 
turn  out  to  be  no  degeneracies  for  F^ff  '  2.  When  M  is  reduced 
to  2.0,  degenerate  mode  behavior  is  suppressed  for  Fcff  13. 

For  the  empty  resonator  case,  the  corresponding  critical  F^jff's 
are  17  and  43,  respectively. 

This  result  is  of  considerable  importance  for  device 
design,  as  it  frees  F,^£f  in  many  cases  of  interest  from  being 
constrained  to  half-integer  values,  and,  by  inference,  permits 
a  more  relaxed  view  of  resonator  alignment  tolerances. 

In  light  of  the  present  results,  one  is  tempted  to  specu¬ 
late  how  the  mode  structure  of  circular  mirror  resonators  would 
be  affected  by  gain  saturation.  For  sucli  resonators,  it  has 
been  shown  that,  in  the  absence  of  gain,  the  in-phase  diffrac¬ 
tive  contribution  from  the  mirror  perimeter  leads  to  strong 
mode  degeneracy,  persisting  well  into  what  miglit  bo  thought  of 
as  the  regime  of  geometric  optics.  ^^9)  certainly  iiot 

intuitive  whether  gain  saturation  could  bo  expected  to  lift  such 
a  strong  degeneracy.  Clearly  this  is  an  area  which  warrants 
further  study. 


E.  IMPLICATIONS  FOR  MMI 


In  an  empty  resonator,  the  poak-to-poak  excursions  in 
the  value  of  the  \q  eincnvaluc  are,  in  the  region  whore  mode 
degeneracy  exists,  'about  equal  to  the  aeomotric  optics  limit 
symmetric  mode  separation.  (See  Figvire  12  or  the  figures  of 
Refs.  14  and  16.)  Since  as  F^^f  -►  ■«,  '  1  and  j  »■  ~  1/M---J. 

This  amounts  to  a  fractional  variation  in  the  output  coupling  of 


27.  Santana,  C.,  and  Folscn,  L. ,  Appl .  Opt.  17 ,  pp.  2239-2243 
(15  July  1978) . 

28.  McAllister,  G  .  ,  vt  ji  ,  IKF.F.  J.  0\iant  .  F.lcct.  Qj  >  l_q , 
pp.  346-355  (March  1974). 
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which,  for  M  -  2.9  is  about  25%. If  saturable  gain  effects 
are  included,  on  the  other  hand  this  viiriation  is  clearly  cut 
at  least  in  half,  depending  on  the  value  of  GqL. 

Even  though  these  results  are  derived  for  varying  Fgff 
with  e  =  0 ,  we  now  present  an  argument  that  they  apply  equally 
well  to  the  varying-e  case,  with  constant,  because  the 

diffractive  processes  involved  are  essentially  the  same. (31) 
Consider  the  tilted  resonator  shown  in  Fiqure  50,  with 
£  =  Cq  ^  0.  We  want  to  consider  two  effective  P’resnel  numbers, 
derived  from  the  portion  of  the  resonator  on  either  side  of  the 
optical  axis: 


(210a) 


p(2)  ^  p(0) 

^eff  eff 


(1  -  e) 


(210b) 


if  e  is  small  compared  to  unity,  then 


(211a) 


(1  -  2c) 


(211b) 


Thus  as  c  increases  from  zero,  the  two  half-resonator  Fresnel 
numbers  move  in  opposite  directions  av;ay  from  F^Qf.  In  parti¬ 
cular,  the  output  couplings  associated  with  each  half-resonator 
follow  the  usual  ooriodic  curves,  as  indicated  in  Figure  51,  so 
that  when  2  c  f(01  is.  an  integer,  f(3)  and  F  ( J  t  are  eacli  an 


that  when  2  c  fJ^^  is 
integer  away  from  F^^^ 


The  observation  tliat  the  next-to- lowest- 1  or 
a  geometric  optics  asymptole  of  ~  ^0-25 
previously,  but  is  apparent  from  a  close  e: 
1  X  I  vs  F(,f  f  plots  . 

A  more  rigorous  exposition  of  tlie  contril)';1 
effects  of  tlie  two  edges  can  be  fovind  in  C 
L.  Felscn,  Apnl._Opt^  11,  t'p.  23'''?- 2  357  (1 

Curiously  enougji,  that  i^aper  conr'.idors  a  f: 
(c)  with  varying  F^f-j-,  Ijut  does  not  take  uj 
ary  case,  altliougli  tlie  formalism  for  doin'! 
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As  justification  for  this  argument.  Figure  5?.  shov.'s  nn 
empty-resonator  calculation  of  i  j  vs  tilt  in  unit.s  of  2  r  F,', j  f 
for  F^p^  =  5.37  3  and  H  =  2.9.  Tlie  validity  of  this  way  of  look 
ing  at  tilted  resonators  can  also  be  discerned  from  an  examina¬ 
tion  of  the  mode  patterns  themselves.  In  Figure  53  we  compare 
an  exact  calculation  of  =  9.3,  M  =  1.9,  2  l  F^ff  =  1  with 

a  composite  formed  by  the  left  half  of  a  mode  pattern  for  F^^^j. 
8.3,  M  =  1.9,."  =  0  juxtaposed  with  the  right  half  of  a  mode- 
pattern  for  F^.p}-  =  10.3,  M  =  1.9,  c  =  0.  The  resemblance  is 
quite  closed  indeed. 
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SECTION  VI 


STABILITY  ANALYSIS  OF  A  SUPERSONIC  CHEMICAL  LASER 


A.  INTRODUCTION 

The  objective  of  this  study  is  to  develop  an  analytic  r'.odfl 
of  the  stability  characteristics  of  a  supersonic  laser;  in  parti¬ 
cular  that  of  a  chemical  laser.  The  flowfield  and  optical  inhomo- 
geneity  of  the  chemical  laser  make  it  an  extremely  difficult 
problem  to  analyze,  and  it  becomes  necessary  to  simplify  tiie  laser 
characteristics  in  order  to  obtain  a  first-order  model  of  mode- 
medium  interaction  inside  tlie  cavity. 

The  chemical  laser  generically  consists  of  an  array  of  flov; 
nozzles  from  v;hich  alternate  layers  of  fuel  and  oxidizer  are 
ejected  (Figure  54).  This  loads  to  a  series  of  mixing  and  reac¬ 
tion  zones  across  the  laser  cavity.  The  subsonic  las(.^r,  as 
described  previously,  can  couple  to  transverse  acoustic  waves  and 
result  in  instability.  The  supersonic  laser  v.’ill  generally  not 
have  such  a  mechanism,  as  the  lateral  transit  time,  D/cs^  is 
much  longer  than  the  axial  flow  time  1/u,  such  that  acoustic  dis¬ 
turbances  will  be  swept  out  of  the  cavity  (bc?o  the  glossary  of 
terms  at  the  end  of  this  section) .  For  the  purposes  of  this  work 
the  flowfield  is  treated  as  being  one  dimensional,  and  therefor 
any  coupling  of  disturbances  in  the  flowfield  to  the  lasing  wil 
be  primarily  axial  in  form.  This  ignores  t))c  possible  instability 
mechanisms  introduced  by  the  series  of  planar  mixing  layers  acting 
as  diffraction  gratings  or  prisms  v.’hich  could  then  interact  with, 
the  laser  flux.  This  mechanism  is  considered  to  be  outside  the 
scope  of  the  present  analysis. 

B.  MODEL  ASSUMPTIONS 

As  stated  above,  our  model  treats  the  laser  as  :f  it  has 
a  one-dimensional  flov/field.  This  is  obviously  not  true,  but  it 
leads  to  an  overall  simplification  of  the  problem  tliat  may  none¬ 
theless  )3e  reasonably  valid.  The  consequence  of  mixinej  is  to 
reduce  the  overall  rate  of  chemical  reaction  in  the  laser.  Tlius, 
a  reaction  of  the  form 


is  described  by 


A  +  B  ->-  0 


dX 

-dF  = 


where  X  denotes  mass  fraction,  and  p  is  gas  dens.ity. 
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Figure  54,  Schematic  of  Laser  Gcometr 


Under  turbulent  mixinq  condition;'.  th('  local  valuen  of  the 
concentrations  [A]  and  [D]  may  vary  su!).sLantia  11  y  v.’ith  time,  and 
their  product  more  so.  Altlioucjh  the  (.lyr.amics  ot  the  mix  iiuj  pro¬ 
cess  are  very  complicated,  and  many  v;oi:ker;;  have  i  nvc;;t  i  qa  t  ed 
such  flowfields,  it  would  appear  feasible  to  use  a  simple  model 
to  estimate  the  reaction  rates  for  a  quasi  one-dimensional  treat¬ 
ment.  In  the  model,  the  concentra tions  X^y,  Xg  of  species  A  and 
B  (averaged  over  a  plane  normal  to  the  flow  direction)  are  giver, 
by  their  mass-averaged  flow  rates  and  an  effective  reaction  rate 
as  follows: 


dX 
_ c 

dt 


kfAX^X^P 


where  A  is  a  mixing  parameter  which  is  typically  much  loss  than 
unity,  and  kfA  is  the  effective  reaction  rate.  The  rate  of 
disappearance  of  f uel /ox id izer  species  can  bo  related  to  a  lasing 
medium  length  in  the  flow  direction.  Thus  for  the  C?  laser,  with 
the  D2  concentration  significantly  greater  than  F2 ,  we  can  write 

X„  =  X  exp  (-k^pXX  t) 
o  2,0 

v.’l'.ere  Xp^  is  the  initial  fluorine  concentration.  Thus,  we  have 
an  e-folding  time  for  disappearance  of  fuel  given  by 


T 


k^pXXp 


2,0 


The  convective  distance  i  traversed  during  this  time  is 
thus  given  by 


£.  =  UT 


where  u  is  the  convective  gas  velocity.  For  a  given  mixing  rate 
we  can  thus  relate  the  length  of  the  lasing  mecium  in  the  flow 
direction  to  the  kinetic  rate  and  the  flow  velocity,  such  that 


u _ 

kj  Ji.pX^^ 


2,0 


In  the  following  we  denote  the  effective  reaction  rate  by 
kf  =  kf  A . 

A  chemically  pumped  DF  laser  exhibit;',  multiline  operation 
sucli  that  laser  action  can  occur  on  an  many  as  th.i'f'o  v  i  lira  t  i  cxia  1 
transitions  and  a  total  of  15-20  vibra t  i  ona 1 - rota t i ona  1  linos.  (12) 


32.  Gro.s;;,  R.W.F.  aiui  Spencer,  P.J.,  ”CW  llyd rog(.n-Ua  1  i do  Laser;'.," 
Handbook  oj  Cliem  i  c^ij_  La  :'.ers ,  R.W.F.  Gros.s  and  J.F.  Dott,  cd .  , 
Wiley- Intersci  once ,  New  York,  NY,  1976,  Chapter  ■) . 
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In  order  to  ^^ccuratcly  model  such  a  system,  the  time  evolution 
of  all  lasing  states  and  the  flux  for  caci>  lasing  transition 
should  be  followed.  Tliis  leads  to  a  system  of  approximately 
20-30  coupled  differential  equations  which  must  be  ijitograted 
numerically.  Such  a  largo  sot  of  equations  is  too  cumbersome  to 
be  of  use  for  the  present  application;  therefore,  simp'lifying 
assumptions  must  be  made  in  order  to  provide  a  manageable  kine¬ 
tics  scheme. 

The  first  assumption  that  we  make  is  that  three  vibrational 
levels  are  sufficient  to  model  the  behavior  of  the  lasing  species 
These  levels  are  the  upper  lasing  level,  OFy,  the  lower  lasing 
level,  DFi ,  and  the  vibrational  ground  state,  DFg.  We  assume 
each  level  is  in  rotational  thermal  equilibrium. 

In  order  to  model  the  lasing  action  in  the  system,  it  is 
assumed  that  the  multiline  operation  of  the  DF  laser  can  bo  de¬ 
scribed  by  a  single  transition  from  DF^  to  DFi.  Implicit  in  this 
assumption  is  the  additional  assumption  that  no  lasing  occurs  to 
the  ground  state  of  DF.  This  allows  us  to  use  a  steady-state 
approximation  for  describing  the  unperturbed  lasing  medium,  since 
with  this  assumption,  the  ground  state  of  DF,  whose  population 
increases  with  time,  is  decoupled  from  the  remaining  kinetic 
equations.  Since  only  ~  20%  of  the  laser  energy  in  a  DF  laser 
arises  from  the  v  =  1  ^  v  =  0  transition ,(  32  ,  3 3)  this  assumption 
is  not  unreasonable. 

The  reactions  necessary  to  describe  the  DF  laser  are  as 
follows : 


(i) Formation: 


F  +  D2  DF^  +  D 


-*■  DF^  +  D 

No  ground-state  DF  is  formed  by  this  reaction, 
(ii) Deactivation: 

^D 

DF  +  M  DF,  +  M 

u  1 


-D 

DF,  +  M  -»  DF  +  M 
1  9 


33.  ”NAC1,  Program  I’inal  P^gToJ^l"  (U)  ,  LTM-2‘J2,  THW, 
Classified ;  CONFIDKMTIAh. 


July  R,  1970. 
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M  is  the  combined  dcnisity  of  DP,  r)2 ,  D  and  F;  deactivation 
by  the  carrier  gas  is  assumed  negligible. 


( i i i ) Lasing: 


DF^  +  hv 


g 

->  DF^  +  2  hv 


Three-body  recombination  is  assunu^d  ncKjlibible  and  only 
single-quantum  dcacti vations  are  considered.  Estimates  for  tlie 
formation  £ind  de>:ic  t  iva  t  i  on-ra  te  constants  wei'C’  obtained  from  th.c 
recommended  values  of  Co!u'n{34)  by  assuminn  that  the  n2/F  ratio 
was  2-3  and  20-30  o  of  the  F  atoms  were  reacted  in  the  lasincj 
region.  Our  estimated  rate  constants  and  tiie  entluilpies  for  the 
various  reactions  are  listed  in  Table  1. 


One  more  quantity  is  necessary  to  model  the  behavior  of 
the  lasing  medium  and  that  is  the  gain  per  unit  length,  which 
can  bo  written 


where  is  the  stimulated  emission  cross  section  and  .‘Nv/j 
the  population  inversion  for  the  v,j-l  v-l,j  transition. 
Writing  this  equation  in  terms  of  the  species  discussed  above 
yields 


g  =  o ( [DF^l  -  0 [DF^] ) 

where  the  brackets  denote  species  concentration  and  o  is  the 
product  of  the  stimulated  emission  cross  section  and  t!ie  frac¬ 
tion  of  the  total  DFu  population  which  is  in  tlie  lasing  rotation 
al  state.  For  a  P-branch  transition,  ••  -  exp  (-2J/Q),  where  Q 
is  the  rotational  partition  function.  Our  c-stimated  values  for 
0  and  0,  obtained  using  the  sp>octroscopic  d  ita  contained  in  Ref. 
35,  also  are  contained  in  Table  1. 

Using  the  model  discussed  above,  the  ♦  ime  evolution  of  the 
chemical  species  and  laser  flux  in  a  cw  DF  laser  is  described 
by  the  following  equations: 


34.  Cohen,  "A  hiief  Review  of  Rate  Coe  f  f  i  c  i  t  ii  tr;  for  R('actions 
in  the  D2-F2  Chemical  System,"  TR- 00 7 4  (  4 5 30 ) - 9 ,  Aerospace 
Corp.,  January  1974. 

35.  Hmanuol,  G.,  "Numerical  Model  inn  of  Chtcmical  Lasers," 

Ihindbook  of  Chemical _ op.  cit..  Chapter  8. 
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TABLE  1.  CONSTANTS  FOR  SIMPLIFIED  DF  MODEL 


Constant 

^Fu 

^F1 

^Du 


Best  Estimate 


8.3  X  10  cm  /sec 

9.4  X  10  ^^cm^/sec 

1.5  X  10  ^^cm^/sGC 


Probable  Range  Heat  of 

5.9x10"^^-!. 2x10"^^  -7.3 

6.6xl0"^^-l. 3x10“^^  -23.4 

3.0xl0”^^-4.5xl0”^^  -0.0 


D1 


8.6  X  10  ^^cm^/sec 


1.7xl0"^^-2.6xl0"^^ 


-8.3 


a 


5.0  X  lO'^’^cm^ 


2.5xl0"^®-2. 5x10"^^ 


0  0.5 


0.4  -  0.6 


Reaction 


kcal/molc 


kcal/mole 


kcal/mole 


kcal/molo 


.  AVCO  I  vi.;nr_T-i 


d[DF  ]  , 

- ^  =  k  [FllD^]  -  (DF^llM)  -  -"-J  {[DF^j]  -  OIDF^)) 

d  t  u  u 


^  fief’ll  n,‘, 

-^=  kp^lFJlD^]  +  k^^[DFj(M]  -  k^^[DF^][M]  +  j-  (  [  DF  J  -  0  [  DF  ^  ]  ) 


d[DF  ] 

- 2_  =  k  [DF,  ]  [M] 

dt  1 


II  =  ca([DF^]  -  0[DF^] )« 

where  c  is  the  speed  of  light,  hv  is  the  laser  photoii  energy, 
and  ())  is  the  intracavity  flux. 

C.  ANALYTIC  MODEL 

The  laser  cavity  may  be  described  by  a  set  of  one-diriiensi  ona  ] 
flow  equations  plus  gain  and  flux  c'quation.  The  flowfield  is 
assumed  to  be  inviscid.  Thus  the  continuity  equations  are: 


Continuity  ^  ^  5  (ou)  ..  „ 
3t  3x 


(212) 


Momentum 


9  u  .  ,  3  u 

P  vr  +  PUs-a 


(213) 


Energy 


where 


+  K,,  Pk’  X, 

ml 


9  9 

9t  9x 


.  ;  Avco  f  VI  nf  11 


In  addition,  tho  gases  arc  assumed  perfect  giving 

P  =  pRgT  (215) 

Tho  reacting  species  may  bo  represented  Ijy  tlic  three  DF 
states  to  give  the  species  continuity  equations 


Upper  State 

DX 

u 

Dt 

K,.,^o.Xp  Kj,^PVu  -  r;  'b 

-  OX^) 

(216) 

Lower  State 

DX^ 

Dt 

bybdiy  "  RJ  "u  - 

rx^) 

(217) 

- 

K,s  OXdi 

Dj^  mi 

Ground  State 

DX 

= 

Dt 

K_  pX  X, 

D,^  m  1 

(21  8) 

T)ie  lasing  may  be  described  l.iy  tl'.e  flux  coupliiui  equation 
described  in  Section  II  and  tho  gain  ec]uation,  such  that 


Flux 


Gain 


St 


<Ij' 

2L 


[/cipu'dx  +  /u/j'dx] 


0  0 


g  =  (X  -  OX  ) 

DF 


(219) 

(220) 


by 


The  magnitudes  of  tho  coefficients  in  Eg.  (219)  are  given 


a 


R 


2nP 


_ ^0 _ 

^ATH 


b 

A 


Qj.  1.0 

These  cquatioiu:.  fully  des.cribe  Ih"  f  Ujw  inr.eU'  i  l.i’  cMvjty. 
For  tho  puriiosfv;  of  a  stability  analyr-is  va  •  a  le  i  n  1  i  ■  r '  t  i  ■'!  in 
tho  perturbation  of  the  flow  from  a  s  1  iM' iv- s  t  a  t  <  ■  c' t  i  r -n  .  Thu:'. 
re-writinc[  each  v.i  r  iab  le_a  s  the  sum  of  a  r.le.aly  st  it.,  and  a  fln.'- 
tuating  p)art  (i.e.,  u  u  +  u')  we  can  deiive  p'  t  t  ui  !  a.  t  i  eii 
tions  from  th('  <iljovo,  ni-alecl  ing  tiigfi'i'  (ad,  i  f,  a:-,:,  ii,  flue;  u  it  i:. 
quantities.  For  t  Ik'  j 'Urpic^;-..  t.  of  this  a  r.a  1  y  s.  i  s.  it  b"e'i:’,e:\  in  e,  s.,a 
to  consider  the'  me.in-valuc  rd  ('ady- s.f  a  t  ,  j  I'op,  i  t  ■  ■  ■: .  <'1  t  h"  f  1  mw 
variables  to  be  const.ant  tliroughaut  the  cavity.  Tlnre  i (  ne 
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exception  to  this  case;  the  ground  state  mass  fraction  obviously 
increases  with  time  whereas  the  other  equations  have  both  source 
and  sink  terms.  The  set  of  perturbation  equations  becomes: 


Contin  uity 


3p'  ^  -  3o’  ^  -  Du' 

at  ^  ^  a-^  +  P  ^ 


(221) 


Momentum 


-  au'  ^  —  Du'  Dp' 

Pat  PP  ax  =  -  a  X 


(222) 


Energy 


(Kp  AHp  +  Kp  AHp  )  {X.,  +  Xp  ■” 


)}(X'  t  X>  . 


^  P^i  ^  P’‘V 


(223) 


where  and  Xp^  arc  the  initial  mass  fractions  of  oxidizer  and 

fuel  respectively. 


State 


P’ 


pR  T'  + 
s 


P ' 


(224) 


r  ; 
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Spccios 


Upper 

DX' 
_ u 

Dt 


n  ^  O 


-  5(X„  F  i  K,  -(X  F  X^F  X  ))(X'„F  X'^F  X'^) 

u  o  o  “ 

DX'  , 


Dt 


Lower 


DX' 


Dir'  “LWh  -  2<V  -'V  ^1+  Sg”F' 

o  ^  o 


'  *  2  >^F  =‘'g> 

loo 


^m>»'  *  X" 


rs  "V  e>;'u-5>^'i)+  3-'--2”:t' 


\^W’''-  ''Dj'”'g.  "'l 


Ground 


DX'  „  _ 

=  ^'D,  [X-  X,o'  +  X  p  X'l  ]- 
Dt  1ml  m  l  Jx 


Gain 


(  (X^j-  0X^)p’  +  pX’^j  -  fOX'j  ^ 


pX^O  In*' 


T  ’  ] 


DF 
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(22! 


(22 


(22 


(22 


r  ,s  I  I 


1 


Flux 


ii' 

dt 


2L 


[/'  Uj^p'dx 

o 


/  a^g-dx] 


(229) 


As  wo  are  considering  supersonic  flow,  we  need  to  only  con¬ 
sider  the  flowficld  downstream  of  the  nozzle  bank;  thus  a  Laplace 
Transform  treataicnt  appears  tractable.  We  apply  a  double  trans¬ 
form  in  both  space  and  time;  i.e., 

f(s,k)  ~  f  j  e  e  f(t,.x)d.xdt 


This  produces  a  set  of  linear  equations  from  the  above 
partial  differential  ones.  Thus  Eqs.  (221)- (229)  give 


Continuity 


Momentum 


(s+uk)  p+pku=0 


p (s  +  uk) u  +  kp  -  0 


(230) 

(231) 


Energy 


1 


-{(K,  AH,  +  K,  AH,  ]{X„  -  J  (X^+  X^+  -  (X^+  X^+  X^)) 


F  F  F,  F,  '  ^'H 
u  u  11  o 


+  [K^  A  ,  +  K  Ali  )  X  +  X,]X  }~  p 
D  D  D,  D,  u  1'  m  c^  ^ 

u  u  1  1  P 


UlKp  AH  +K  All  HX  +  2  X  -  (V  ^1*  IK  AH  +  K  AH  1  X,p 

uull  o  o  uull 


X  {X  +  X.  } 
c  u  1 
p 


(232) 


+  (s  +  uk)T  -  k u/CpP  p 


+  ^r'S-  ‘'"r  *  h-  ''"p  1»H  *  b'  ^g>> 

puullo  o 


=  0 


/3/ 


,  .'AVerj  i  vr.FiL  :  r 


state 


—  p+pR  T-p=0 

Y  g 


(233) 


Species 


Upper 

State 


1,5 


U  O  O 


1  X. 


+  [(s+uk)  +  Kp  pCXp,  +2 


F  "  “‘F 
o  o 


(2  34) 


+  [(s+uk)  +  Kp  p(Xj,  +  I  Xp 


u  o 


o 


*i'f  *  -^F 

u  o  o 


g  g 


Lov/er 

State 


-IK, 


f/\-  I‘V  Xj+  x^)H- 


+  (K  P(X  +  2  S  '"''''rn  hv  ^u 

^1  ^o  o  ^ 


(235: 


+  [^'f  P(^h  +  I  ^F  ^1'"'  ■"  hv  ^  ‘''I 

1^0.0  1 


X  1  X. 


^'ifv  "f-'  ^ 


F  K  MX,,  +  2  X,,  -IX„x  K^-.  X  ))  X 
1  o  o 


Ground 
St  ate 


-(s4.2uk)Fp^X^Xjr  -  Kp^X^^c  X^uk  ^  ukdUu;)  ^  0.  (^it 


I  7n 
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r 

n 

i 

« 

i 

i 

i 

i 

i 

« 

i 

« 

M 

i 


Ga  in 


g  -  j.r-  [(X^  -  OX^)n  +  P  n 


(  In' 


^  ^  (  2  37) 


Flux 


sC  -  {o’j^  0  +  n^gl/k  ^  0 


(  2  3  y ) 


Fl  iini  nat  imi  n  from  Fqs.  (230)  and  (231)  aiul  g  f  ro;:i  Kq.u. 
(237)  and  (238)  wo  obtain 


Con  tin  ui  t\'-Monien  turn 


( s  +  uk ) ^  p  -  k^p  =  0 


(239) 


Gain-Flux 


«jO(X^-  0Xj_))p  +  -  ksl 


-  cl  -X  C  ,  “  cl  a^oXl''  lr.‘-  . 

2L  ®  2L  — - ^  ^ 


(240) 


whore  a  =  . 

The  seven  equations  in  variables,  p,  T,  p,  1,  X^,  X^., 

can  be  written  in  matrix  form  ■' 


[A]  (xl  =  0 


(241) 
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Thus  v;c  have 


Transf oiT.ieO  lAiual'-ioi^s 


a 

11 

0 

0 

0 

a 

15 

"21 

"22 

-ks 

"24 

Ctor 

2  D 

(s+uk) 

2  0 

0 

0 

0 

"41 

"42 

0 

44 

(s+uk) 

"51 

^52 

"53 

"54^ (s+uk) 

"55 

"61 

„+ (s+uk) 

0  *i 

0 

a64+ (s+uk) 

0 

"71{s+2uk)  ° 

0 

“74“’' 

0 

0  0 


0 


46  47 


0 


0  uk  (s-^uk) 


(242) 


where 


11  "  “l5  ~  “16  “ 


=  El 

21  2L  ‘”R'  "I"  '"u 


u,  =  ^  (a.+  a,o(X„-  0X^)1  ,  CX 


=  Ell  a  p 
22  2L 


24 


2L  '  "25  ^  “l  iL  «  — -— 


T 
1  /n 


41 


=  -  1-{IK  Mip  +  K^,  &Uj,  KXjj  -  ■^-(x^+  x^+  -(x^+  X^+  X^)} 

'■p  U  U  1  i  O  ^ 


u  u  11 
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42 


“44  "  ^  ^ 


F  “‘‘F  "I  , 
u  u  F 


+  K  '‘''1  ]{X  +  '2 

i  t  ft 


1  "o 


46 


u 

pc 


a  =  ~-{“u  +  4  ^  ^’’'F  ''”f  +^'F,'^‘”f,^ 

47  ”o  ^  o  “  ^  ^  u  u  1  1 


P  47  Cp  Q 


-iKp  (Xjj  -  -2-  (X^+  X^+  Xg))(Xp  -{a^+  X^+  Xg)+  Kj  X^X^ 


51 


1  o 


“52  "  ^F  I  ^F  “  ^m  hv 

I  o  o  ‘-’ 


“  (X  -  ex,) 


t__  —  —  r'—  \J\  \j  /\  ■>  j 

53  hv  u  1 


=>'f  5(>:„  "  l>=F  -<V  V  ^>>  *  ^ 


91n0 

“55  “  ~  hv  ;j 

“57  "«f  *  i\-  "'V  %’> 

^  1  o  o 


»61  =  -■'p  <>  '^.  -  I  ‘V  h*  h*  ^g>> 

”  u  o  ^ 


■62  =  “64"“67=%'^<V^V‘V  V  V> 


“71  "  '  ”74  ' 


.  :  AVCO  FVf  fvr  T  1 


The  stiibility  of  the;  sys;teni  can  l)c  examined  in  terms  of  tiie 
characteristic  equation  derived  from  He.  (241);  i.e., 


det(A]  -  0 


(24  3) 


An  analysis  stability  of  the  system  is  obtained  from  t)ie 
Routli  criterion  for  Ilq .  (  243).  (36) 

The  characteristic  equation  may  be  detained  from  Eq .  (242) 

After  tedious  algebra  this  reduces  to 


a^s 

6  ,  5 

S  -r  a  s  s 

+  a^s^  + 

a3s3 

+ 

2 

a2S 

+  a^s^ 

'''  '"'o 

where 

«7 

^6 

6riAg  + 

15n^A^ 

(d 

+  SriAg  + 

A4  +  Bg 

^  = 

20n^As 

+  lOn^Ag 

+  4  t)A^  + 

A3  + 

SnBg  + 

^4 

^3  = 

+  lOii^Ag 

+  6nA^  + 

3r,A3 

+ 

A2  t 

lOn^Bg 

= 

6r,\ 

+  Sn^^Ag  + 

4n^A^  + 

3.,S,3 

+ 

2qA2 

+  A^^  + 

lOn^Bg 

+  6n  +  3nD^  +  B2 


6  5  4  3  2 

-  h  As  +  n  A5  +  n  A.  +  n  A3  +  n  'A  +  ha,  +  a 


0 


+  Si/ds  +  4n^B^  +  Sn^B^  +  2nB2  + 


5  4  3  2 

Uq  —  ri  B^  +  n  +  n  B^  +  n  B2  +  nB^  +  Bq 


36.  DiStefano,  J.J.,  F.tubberud,  A.R.,  William';,  l.d.,  "Feedback 
c\nd  Control  .Systems,"  McGraw-Hill  (1  967). 
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r 

n 

1 

f 

f 

i 

« 

i 

i 

I  ■ 

ir 

i 

t 

« 

« 


and 


Ae  =  -ct^gkn 


^  =  “15^46’"  '1  -“16^62’^^  +“16^“52  " 


^5  “l6'^53  ^“22  °‘24'^ 


A4  n  '^15«46^“52  ~  '*‘“16^62  ^“52  ^54 

+“16^57^74^^ 


®4  “53^'^22  “  *^24^  ^  '*'  ^^62^ 


^3  ~  ^11*^  ''^^°‘52  ®54^~^62^‘''”l5''^46”57"^74^ 

■^“l5“47“71^  -  cc^^a^^nk  +«i6“57“62“74^'"^ 


““  16*^4  7*^5  5“  7  4^'^  ”l6“52”6  2”74'''" 


®3  “ll”53^“22  "”24^'^  '*'  ‘^15“46“53^62  ^'^'‘22  “«24*^’^ 


"“l6^2  5^4  7^53^74'^  ‘^’‘16“2  2”5  3”g2"7  4'^ 


)kr, 
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^“l5“46“57“62'‘74^^^  ■"  ^15^47"62“71^^ 


~“l5“47“71  ^“52  "  “54^'^  “l5“47'^7i’"'  "“l5“42°‘62 


‘15  47-^71' 


■*'“l5^'*4l''52  "'"4l"‘54^'"^'^  «1  S«4fi«52«fi2^74^^’' 


'15“46‘^52“62  ‘74' 


-“l5“42“62“7l'"^  ^  “l5“42"6l’^^'l 


■*'“l6“55°‘62“74  ^“42  “  047)^^'^ 


“ll“53°‘62^®22  ~  ^24^^  ~  “l5^47^53”71  ^^22  "^‘24' 


+'=^15^41^53  ^«22  "  '^24^'^  ^  «15“22“46«53“74’^^ 


■*'°‘l6®25°'53°‘62°‘74  ^“42  “47^^^ 


®ll‘^47“55“74^^'^  “  ‘^il“52“62“74^ 


■‘■“l5“47“62”71^^'^  r  ^15'^47^62“71^“52  "  “54^^ 


-“15«47“71<«52  ■  “54>'^  ^  r  “is"?  4 '^'4  l"5  7^5 1'M  7  ^ 


+“15“62«41^”52  ■  "54^’^^^  '  “l5'=^4 2^'g 2^'? 


■^“l5°‘42”62''‘71  ^"52  ~  “54^^'^  ~  "l5"42'^61 


1  “11^25  47  53  74  '  11" 22  53  G2  74 

“°‘l5°‘47°‘53“71  ^*^22  ”  “24^^  ^  '  "l5‘^47‘^‘53^G2“71  ^^22  ~  “24^*"' 

““l5“2l“4' ”53*^74’^"  ^  “i5^41'^53”G2  ^“22  “  ''' 

■*’“l5^42“53”61  ^^22  “  *^‘24^^  ^  ■*■  “i5“42"^53”G2"'‘71  ^^22  ~  °24^^ 


0  ==  ■“ii“^7“55“G2“74^^^  '  '‘l5“47“62«71  '"'52  “  “54^^'^’^ 


■^“l5“62“74  ^“41*^57  “  “51°‘47^^ 


■’'“ll“42“55“62^74^  ^  '^15'’61°‘74  ^°‘42'^57 


’'“l."'^62“74  ^'^4i“52  “  “5l“42^^" 


■^“l5“42®62'^71  ^“52  ~  “54^^ 


“'^ll"25“47”53®G2*^74’^  “  °‘i5"2l"47”53”G2“74^" 

~“15“47“53"G2‘^71  ^“22  '  “24^^‘  ”li'^2  5*^4  2*^5  3”g  2”7  4’^  ^ 

■*'“l5“53“74  ^”22®47“61  “g 2  ^*^2 1*^*4 2  "”22^41^^^  ’^' 


■*'“l5“42“53”62“71  ^"22  “  “24^^’‘ 


Tlic  value  of  the  spatial  wavenumber,  k,  is  taken  to  be 
given  by  the  cavity  dimension  in  the  flow  direction,  ■.  Thus, 

k  ^  1  /  f.  .  . 

Thus,  v.e  can  nov/  numerically  specify  all  the  cijef  f  ic  i  ents  f)f  tliL 
character i s t ic  equa t i on . 
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D.  STABILITY  ANALYSIS  OF  DF  LASER  MODEL 


The  values  of  the  laser  parameters  can  bo  obtained  from 
the  "gain  equals  loss"  condition  and  Eqs.  (21G)  and  (220).  The 
"gain  equals  loss"  condition  specifies  the  gain  in  terms  of  cav¬ 
ity  loss  and  cavity  length  such  that 


e2gL 


^2 


where  rj  and  r2  are  the  mirror  reflectivities.  Thus  wo  obtain 
the  population  density  of  upper  and  lower  states  from 

(X^  -  OX^)  =  g/cp 


The  flux  can  be  obtained  from  Eg.  (216)  where  we  assume 


DX 

u 

Dt 


^  0 


ex^) 


Thus  4i  =  Kp  P^X„  Xp,  ohv/(og) 

^o  °2,0 

Using  these  criteria,  the  parameters  for  a  wide  range  of 
flow  conditions  were  investigated.  The  laser  conditions  arc 
listed  in  Table  2.  Over  the  entire  range  of  mixing  rates,  flow 
velocities,  cavity  lengths  and  pressures,  use  of  the  Routh  cri¬ 
terion  specifies  that  the  laser  is  stable.  This  can  bo  simply 
described  by  reference  to  the  coefficients  of  Eq.  (244).  The 
only  matrix  coefficient  containing  contribution  from  (the 
density  fluctuation  induced  flux  variation)  is  ''t2i  wliich  appears 
in  coefficients  and  Bq;  those  affect  only  coefficients  a]  and 
ao-  In  general,  Bi  and  Bg  are  significantly  smaller  than  the 
other  terms  in  ai  and  ag,  resulting  in  both  a)  and  ag  being  j>osi- 
tive  (if  either  of  these  coefficients  were  negative  an  unstable 
root  would  occur) . 

The  response  of  the  laser  to  upstream- induced  noise  dis¬ 
turbances  may  be  described  in  terms  of  transfer  functions  derived 
from 


lA] 

[X]  - 

IF] 

where 

[F]  is  the 

fore 

i  n  g 

f unct i on 

vector . 

The  transfer  function 

between  parameter 

^n 

and 

J'n 

is  given  Ijy 

det. 

(A,,l 

f” 

det” 

r^r 

where 

[Ap]  is  tile 

ma  t. 

rix 

[A] 

wi  til 

the  nfh 

column  replaced  by  [ E j 
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TABLE  2. 


LASER  CONDITIONS  INVESTIGATED 


Cavity  length  (L)  0.1  -  3.0  m 
Mixing  length  £  1.0  -  10.0  cnis 
Flow  velocity  10^  -  3  x  10^  m/sec 
Cavity  pressure  1-10  torr 
Cavity  density  5xl0"^-5xl0^  kg/m^ 
Laser  mix  10  He/5  D2/O. 8  F2 
M  =  5.45 
Y  =  1.59 

Mirror  Reflectivities 

=  0.98  r2  ~  0.88 
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Tht;  rcspon.'U'  oi  ''.uch  a  .syr.tc-in  i:;  t-lnia  ch.na  f  t.  f  ■  r- i /Lt!  ) 
l/dot[A].  The  norn\a  1  i /< 'tl  macjn i  teacJcs  ol  thJ^;  fuj.el  i'':!  arc'  r  1  1  u:; 
trated  in  Fiqure  5S.  The'  r  esn-'onne  in  i-ne.can  1 1  a  1  I  y  Liiiifor:;i  uj)  t.o 
f  requeirc  ion  of  lO'j  I'ad/nec,  and  t  hera/a  i  l  cr  r  iredly  c; ,  >f  ra 'a  r.  ■  in. 
magnitude  with  a  ;;loi.>o  ariyinptoting  to  -  ‘  .  The  resjeonno  reduce 

in  magnitude  as  cavity  pressure  is  inerc-as'-'d . 

Obviously,  this  is  only  part  of  tiie  overall  tr'ansfer  f unc 
tion,  and  further  analysis  is  necessary  if  tlie  detailed  system 
response  is  to  be  determined. 

E.  DETAILED  SYSTEM  RESPONSE 

The  model  developed  here  is  a  simplification  t'f  tlie  re.i] 
laser  behavior.  A  more  detailed  and  cor  rc.'siH';.d  i  nal  y  diffn-ult 
approach  is  required  if  the  full  tiiree-d.  ■ '’i-.'n;;  i  c'na  1  iiafi>-(.-  of  tii 
mode  medium  interaction  is  to  be  inve.st  i.ejat o-d .  ;))<'vio 

analysis  a  simple  acoustic  or  gain  fluctuation  inn.tibility  does 
not  appear  to  occur  in  the  chemical  la;;er.  Hov.'evo!' ,  tire  ica! 
flowfield  may  result  in  other  significant  mocli.anir'.ias  occur!  mj. 

The  chemical  laser  consists  of  a  series  of  'tuas  i 1  an.ar 
flame  sheets  (Figure  54).  The  jets  from  the  norr'e  Ln.nr’.k  n.ro 
normally  in  the  transition  range  of  L.'.minar  to  turl.'V.lent  flow, 
'..'here  laminar  instabilities  of  the  f  1  aaic'.nhoef  .iro  .i:  o.i  rent . 
These  instabilities  arc  akin  to  the  Hay  loigr.-Tay  !  or  i  nr,  t.d; :  1  1 1.  y 
in  mixing  layers  (see  for  example,  Kintura(37)  .ipd  Toong.(3P) 
Fluctuations  in  tlie  local  gain  and  flow  pro'portios  will  rc.'-ult 
in  modulation  of  laser  output  flux.  'whether  such  a  flcwfield 
instablity  can  lead  to  an  output  instability  is  not  clear.  The 
coupling  between  the  array  of  mixing  layers  is  of  i  mport  a;ice  ; 
if  the  flowfield  is  driven  by  an  up.stream  acoustic  dislurijance 
which  is  identical  for  all  the  rihoar  layers,  this  ni.r/  re.su:  t  in 
significant  optical  coujiling.  Ilowevc';-,  a  random  ph.:i  r,  in.. :  of 
oscillation  of  the  array  of  flame  sheets  may  result  in  a  cjuite 
different  behavior. 

F.  SUMMARY 

A  simplific'd  model  of  a  raipers.on  i  c:  chi.-m  i  '  la.s.'T  he  I"'.’ 
derived  and  a  .stability  analysis,  pi'rform.ed.  This,  ind.ie.t.:,  t  i,  i 
the  laser  is  stable  >:K.:ross  a  wide  raiifU'  of  o['i,'r.i  t  i  nq  oo'ud  1  t  ion.: 


37.  Kimui'a,  I.,  'font  h  Eympos.ium  (Inl.)  or.  Co’.nL'as,  t  i  ( -r. ,  p.  lii'tf- 
1  300  ,  The  ComfiUL'. t  ion  Ins.titute  (l‘'()5). 

38.  Toon<i,  T.Y.,  .“'-^ilant,  R.F.,  Stopford,  J.M.,  And.  i.-on,  'MY., 
Tenth  ilyinp;  >;;  i  ur:;  (  I  r;  ‘  .  )  on  Conibus  t  i  o,;)  ,  p.  1301-1  ■:  1  3  ,  Tl;'' 
Combu.sLion  Ins.titut.  (1 '>•!'>). 
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However,  the  s  iinn  1  i  fy  i  nej  aso  t  i  ons  ri'niov-’  the  yor-.;  i i  1  i  tv  of 
a  coupling  occin'riiiq  ar.  a  con.-.i-ciuenco  o!  f  1 .1’ne- ;:)hm  ■  t  i  astcib  j  1  i  ty 
and  as  sucli  would  not  ai'j.ioar  in  Lht>  solutions  obtained  in  tins 
work . 


The  above  analysis  provides  a  useful  startinq  ooint  in 
obtaining  a  description  of  possible  niodo-inod  i  ui;i  iiu.cq;act  ions  in 
supersonic  mixing  lasers.  However,  it  is  far  from  complete  and 
it  is  reconimonded  that  further  explorati  on  of  tlio  t>o;,sible  be¬ 
havior  of  flame-sheet  instabilities  upon  optical  per ,■  ormaj'icc  i.s 
necessary . 


G.  GLOSSARY  OF  TERMS 
a  coefficients  in  Eq.  (244) 

A  coefficients  in  Eq.  (244) 

B  coefficients  in  Eq .  (244) 

c  speed  of  light 

Cg  sound  speed 

Cp  specific  heat 

g  gain 

h  enthalpy 

k  space  transform  variable 

k  kinetic  rate  coefficient 

K  effective  rate  coefficient  (=  Ak) 

1  mixing  length  scale 

L  cavity  length 

molecular  mass 
DF 

r  mirror  reflectivities 

p  c;as  pressure 


u  gas  velocity 

t  time 
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s 


time  transform  variable 


T  temperature 

X  axial  coordinate 


coupling  coefficient  in  Eq ,  (219) 

a  coefficient  in  matrix  Eq .  (242) 

3  Gladstone  Dale  Constant 

Y  ratio  of  specific  heats 

6  degeneracy  coefficient 

X  mixing  parameter 

X  laser  wavelength 

p  gas  density 

*^ATM  density  at  a  pressure  of  1  atm 
0  lasing  cross  section 

<{>  intracavity  flux 

T  reaction  time  constant 

X  mass  fraction 

AH  enthalpy  of  reaction 

0 


Subscripts 
u  upper  state 

1  lower  state 

g  ground  state 

F  formation  rate 

D  dciict  i  va  tion  rate 
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SECTION  VII 

SUGGESTIONS  I'OK  FUTURE  RESEARCH 


The  dynemiic  models  considc-rc-d  in  Sections  II,  1\'  and  VI  Inive 
made  simplifyincj  assurr.!-''tions  on  *:!'.<■  nature  of  the  la.'-.er  cavlvy. 

In  particular,  uniforr.i  valuers  of  ejain,  cicnF;ity  and  f  i 'a;-;  iuw';-  L, 
assumed  in  the  cavity.  .More  qvuintit  a.t  Lvo  ros.ult.s  may  Lcf  c;btainva 
by  perturbing  the  actual  steady-state  motias. 

The  theory  contained  in  Sections  Il-IV,  and  VI  of  th.is  rc'- 
port  deals  not  only  with  ;x5ssiblG  system  i  nst  ab  1 1  i  t  ic-r. ,  but  also 
in  the  system  tran.sfer  fvinctions.  Accord  i  p.!:  1  y  .statistics  of  tb^' 
systc.m  output  can  bo  [predicted  it  e>:t(?rnal  jawturbat  ions  can  be 
modeled.  In  particular,  it  would  be  useful  to  mo’-.V.  the  turbu¬ 
lent  nature  of  the  nonuniforni  flow  field  of  tite  f  1  cv.si  r.^ ;  ga:;. 

Once  the  correlation  (or  distribut  i  r’n)  functi^'n  of  l  be  d(-n;:-,ity 
(i.e.,  refractive  irsdex)  function  i.s  kno'.-.'n,  the  sta  t  j  :;i  t  ics  of 
the  noise  driver  are  known.  T!io  theory  o:  this  re:  art  tlmun 
allov/-s,  via  the  transfer  function,  l.hc  calculation  o:'  tlue  sta¬ 
tistics  of  the  ia.ser  flux  outioat.  In  cxldition,  the  r;'-t  ;';oy;s  rV' 
the  theory  of  light  proR'agation  tlir'ouch  t  uriru  1  ■ 'nee  can  be  used 
to  find  the  statisitics  of  wa  f  ro::  t  di  o,  tort  :  cm  of  t!;o  o\;tr-ut 
v;ave.  This  problem  is  especially  acute  for  th.e  chemical  laser, 
where  the  imperfect  mixing  of  the  reactant  cases  can  bo  moclolc'd 
to  predict  the  statistics  of  th.o  system  inhomocjoneitie.s . 


calculalii^g  the  sc  I  f -oc;wm  s  ton  t  modes  o:  a  ;e‘.  rur.m  e  u  ,  u;.  e  u;.;.  . 
resonator.  This  is  in  fart  the  only  r''aoonably  :ast  i  .-.cf.o.  :ctu'-- 
that  wo  know  of  for  making  .such  a  ca leu  1  at  i cri .  Althc'uu'.i  wo  i.ooec 
only  applied  it  to  a  str  j  p-rco,onator  cioomot.  rv ,  it  appears  t;.i  bio 
genera  1  i  xabl  e  to  a  circular  tu-omiotry  as;  well.  Cti  1  c:u  1  a  t  i  or.  r;  cif 
this  type  can  b('  used,  for  example,  to  detor’olne  the  eoupl.in.e 
between  modc.s  and  thc^-  modification  of  the  coup  linn  of  :n- 

dividu.il  mode.s  induced  by  an  externa]  ;k  r turba  1  ic.'n .  vk'  pri-d  ’  sf 
that  thi.s  techi'.i<iuc  will  find  wi  de- I'ang  i  ng  api.d  i  cat  i  c'pa;  lo  las.v'r 
re.sonator  probleiiis  in  tlie  future. 


(Sm 
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APPENDIX  A 


THE  PERTURBED  RESONATOR 


In  this  appendix  we  will  treat  the  pertnrbation  ex.p'an :: ion 
of  Section  II  somewhat  more  riciorously,  there!,''.’  c  1  e.  r  :  f'j- i  ne  tb.',.- 
pliysical  mcanimj  of  the  acoustically  induced  perturbation  oL  t  iie 
elcctromagiiotic  wave  in  the  resonator. 

Writing  the  part  of  the  electromagnetic  wave  vdiich  propa g-i  t  e ; 
to  the  right  as 


E(x,t) 


u(>:,t)  e 


iik^y.-rot) 


+  c .  c 


(A-1  ) 


we  find  that  u  satisfies 


u(x,t  +  T^)  =  .'■K(x,x',t)  u(k'  ,t)  dx' 


(A- 2) 


where  unperturbed  prepagator  (i.e.,  t'ne 

propagator  in  tlio  absence  of  tlio  acoustic  vaave  I  ;  i  .a  tiie  per¬ 
turbation  due  to  fluctuations  in  tlie  density  (acoustic  v.'ave)  and 
gain  of  the  medium  (see  Section  II);  x  and  x'  represent.  (v<.'Ctor) 
positions  in  an  x-y  plane  between  the  mirrors;  and  2L/c  is 

i^he  round-trip  time.  In  the  following  wo  write  x  to  denote 
X  =  (x,y)  . 


The  time  dependence  in  is  duo  to  an 
whose  frequency  is  very  low  compared  to  I/^q 
means  that 


acoustic  oscillation 
-  5  X  lo”^  sec-1.  Thi 


9K. 


9t 


o 


(A- 3) 


and  wo  therefore  may  construct  an  equation  for  riu/'t  (x,t)  as 
f ol lows ; 


?u 

3t 


( u  { X  ,  t  +  •!  ^ ) 
^  o 


u  (x  ,  t)  ]  = 


+  /  V  ^  ( X  ,  X  '  ,  t ) 


u (x '  , t)  dx  '  ) 


i  l/V^(x,x’)  u(x’,t)  dx' 

(A-l) 


1^ 
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where 


V 

o 

-  (K  (x,x’)  - 

^  (x-x' ) 1 

(A-5) 

and 

V,  .  f  K  (x,x' 
o 

'  ,t) 

(A-6) 

Equation  (A-4) 

has  been  written 

in  a  form  which  makes 

it 

identical  with  the  Schrodinger  equation.  Thus,  we  nay  j;olve  this 
equation  by  the  well  known  per turbation-theory  techniques  used  in 
quantum  mechanics. 

The  integral  operator  represents  an  unperturbed,  timo- 
indopendont  "Efainil tonian "  of  the  system,  while  is  tac  perturba- 
'tion.  To  calculate  the  behavior  of  u(x,t)  (analoyou.^  to  the  waw 
function  of  quantujn  mechanics)  W’e  may  expand  u(x,t)  '  any  complete 
set  space-dependent  functions: 


u  (x,t) 


a^(t) 


V 

■  n 


(X) 


(A-  7 ) 


As  the  set  '-'^c  choose  the  eigenfunctions  of  the  unpc''turbcd 

Hamiltonian  V^,  so  that 


A  1'  (x)  =  .'•  V  (x,x’)  T  (x')  dx'  (A- 8) 

nn  o  '  n 

i.e.,  is  the  eigenf  unction  of  eigenvalue  .  Equation  (A-8) 
may  be  shown  to  be  equivalent  to 


1  (x)  =  (1  +  iA  T  )  /  K  (x,x’)  T  (x')  dx’ 
n  n  o  o  n 


(A-8) 


or,  for  1 


-iA  T 

e  Oy  (x) 

n 


/  K  (x,x'  )  i'  (x’  )  dx' 
o  n 


(A- 10) 


is-i 
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which  is  the  usual  intcqral  form  of  the  oigonvaluv'  c.-quation  cjf 
a  laser  resonator.  It  will  be  shown  below  that  tlie  natural  optical 
frequency  of  mode  n  is 


0)  =  CO  +  Re  A 

no  n 

i.e.,  Re  represents  the  frequency  offset  of  mode  n  from  that 
of  the  pure  longitudinal  mode  of  frequency  q.  On  tlie  other  hand, 
-2  ImAj^  represents  the  fraction^il  energy  loss  per  second  of  mode  n. 

Substituting  Eqs.  (A-7)  and  (A-8)  into  Eq .  (A-4)  gives 

‘  s  fe 

m 


1 A  a 
m  m 


T  (x) 
m 


E 


m 


a  (t) 
m 


/  V,  (;<,x' 


,t)i^^(x 


)  dx' 


We  now  assume  that  the  modes  {i'li}  form  a  com;ilc'tc  ortho¬ 
normal  sot.  This  assumption  could  be  rigorously  justified  if 
the  unperturbed  Hamiltonian  were  Hcrmitian;  tlic  eigenvalue' 

Ap  would  then  bo  real,  corresponding  to  an  unperturbed  cavity 
in  wliich  all  modes  have  the  same  value  of  not  loss  (loss  minur^ 
gain).  This  will  clearly  not  bo  the  situation  in  an  a.ctual  un¬ 
stable  resonator.  ilowever,  thc^;  only  modes  w'nich  are  in. tercstinq 
are  the  initially  oscillating  mode  and  th.ose  models  who.'-.c'  net  losses 
arc  nearly  equal  to  that  of  this  initial  mode;  oth'sr  modes  (wit's 
significantly  higher  nc>t  losses)  will  bo  stroncly  damped  and  'will 
therefore  never  be  excited  by  th.e  acoustic  coupling.  Thus,  it 
should  be  a  good  approximation  to  assume  that  for  tlie  modes  of 
interest 


1mA  H  0 
n 


(A-12) 


and  therefore 


/  T  * 
n 


4’  dx 
m 


=  6 

nm 


(A-13) 


Equation  (A-12)  implies  that  multimodi^  osc  i  1  la  t  i  on  is 
model  las.er.  In  fact,  sucli  mul  t- iir.ode  beh.avic'r 


in  our 


oljsc'rved  in  cw  EPI,  Hascu's 
by  4  *(x)  and  integrii  t  i  nc} 
Eq.  (A-13) 


,(U 

over 


Mu  1  t  i  [1 1  y  ii'.g  belli  sides 
all  X  tlien  yields,  with 


:  'OSS  i  b  1  f' 
lias.  b.e'''n 
el  Eq.  (A-11) 
the  aid  of 


1.  R.  Patrick,  AERE,  private  commun ica I i on . 


ir:x 
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(A-14  ) 


da 
_ n 

dt 

where 

W^ni(t)  =  If  4'^*(x)  Vj^(x,x',t)  dxdx’  (A-15) 


i  A  a 
n  n 


-iEw 

nm 

m 


(t) 


m 


It  will  be  noted  that  no  use  has  been  made  so  far  of  per t  ur):)a t  i  on 
theory;  within  the  limits  of  our  assumption  tint  the  modes  of  in¬ 
terest  are  orthonormal  Eq.  (A-14)  is  exact . 


The  meaning  of  Eq .  (A-14)  may  be  elucidated  by  first  ex¬ 
amining  its  solution  in  tin  absence  of  the  perturbation  (i.o., 
with  =  0).  We  then  find  that  each  amplitude  ajj  satisfies 


da 

=  -  iA  a 
dt  n  n 


(A-16) 


with  the  solution 


-i A  t 

a  =  a  (0)e 
n  n 


(A-17) 


so  that  the  electromagnetic  wave  propagating  to  the  right  in  the 
cavity  is,  from  Eqs.  (A-1)  and  (A-7) 


E(x,t) 


(0) 


c 


i  (k  z-u'  t ) 
o  n 


^^(x) 


+  c .  c 


(A-18) 


where  ,  the  natural  frequency  of  mode  n,  is  given  by 


a)  =  (ii  +  I<o  ( A  )  =  ui  +  A  (A- 19) 

n  o  n  o  n 

This  i  r>  clearly  v;hat  one  would  exi''ect;  each  mode  v;hieh  is, 
present  initially  (at  t  =  0)  cont  imK.'s.  to  os.ci]la*o  freely  it  its 
natrual  frecpiency.  fU’cainsc  we  have  not  allowed  for  any  damj  vi  no , 
none  of  thc^se  modes  decays  with  time;  s  i  lure  we  have  not  allowod 
for  any  coupling  between  them,  tliey  do  not  affect  one  another, 
citlicr . 
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Con.sidor  tho  direct  elfc'ct  of  l  ho  pc'r  turhal  i  on  on  n.n  init¬ 
ially  oncillatinq  mode.  Wc  denote  tho  iniria]  niocie  by  t!io  nuii- 
script  0,  and  take  \q  =  0.  (Thi.n  i:;  ec]u  i  v.:i  1  ( 'ii  t  to  liolininn 
to  bo  tile  frequency  of  the  initi.al  mode.)  bi  rr.  t-ordLT  o.  ■  I'l  in  h .  ; ; 
theory  allows  us  to  innore  the  amp.]  i  tudes  of  all  the  other  mouen 
to  a  first  approximation.  We  thus  have,  from  Eq .  (A-14) 


da 

o 

dt 


i  h’  ( t ) 
oo 


a 

o 
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Tho  real  part  of  (defined  in  liq .  (A-lh))  er'r  r- ■5q>on  ;  n 

to  a  shift  in  tho  freqiU'acy  of  tiie  initie.!  r.iad  ■ ,  wi:’!,  the  in,- 
aginary  part  corresponds  to  a  ch.ange  in  th.-.  ;u't  cavity  loss 
gain)  per  unit  time  of  this  mode.  .sucli  a  U'ss  clian  ’  ecu !  d  co:.,.- 
about  as  a  standing  transverse  acoust  j,c  wave  osciliat.s  in  ampli¬ 
tude;  wo  would  then  expect  the  cavity  loss  to  be  mcshalated  at’ 
twice  the  frequency  of  tile  acou'itic  wave.  Lev..';  modal  at  ion  of  this 
type  is,  in  fact,  used  to  mode  lock  lasers..  !t  could  ,i  n  priin.-i;  e 
load  to  an  instability  if  tho  initial  mod^,'  is  not  as.itorm  in  in'.;.- 
sity  within  tho  cavity,  .since  thr.  heating  a.ssociatcsl  v.-ith  t  h"  'M'.’ity 
loss  modulation  is  then  nonuniform  a.n  i  can  on.  its.  an  acou.'tic  .iis- 
turbanco.  Indeed,  our  calculations  of  t  h«-.  ma.i,.;;  of  an  unst.ibm 
resonator  (Section  V)  siiov-;  explicitly  tint  a  ncnuniid.rm  mods.  i:.w 
sity  distribution  is  to  lie  oxp.:^cted,  so  tins,  is  a  possibK'  candi¬ 
date  for  an  acoustic  instability. 


Lot  us  now  consider  tho  effect  of  the  acoustic  ;  •  g- 1.  urb  1 1  i  on 
on  the  other  modes  of  tho  cavity.  Again,  we  may  assume  to  liisi 
order  that  the  amiilitude  a^  of  the  initial  mode  is  mu.-h  lai-  a.  i 
than  that  of  any  other  mode.  Wo  may  tlius-  icinore  all  amj.ditudi  :; 
a^n  other  than  on  the  right  of  Eej .  (A-14),  and  s.et  e:r.  il  te 
unity.  This  approximation  is,  of  coutmae,  only  .s.'lf  nsis*,.  :,t 
as  long  as  tho  oquatiens  predict  thar  1  f'.  r  all  n  ^  0, 

i.e.,  in  tho  initial  stages  of  any  buildup'  of  otiier  e.c.'ity  ms.,;...  . 
The  aj.j '  s  then  satisfy 


da 


I 


n 

dt 


+ 


iA  a 
n 


n 


i  W 


no 


(t)  a^(t) 


(A-2  1  ) 


It  is  useful  at  this,  point  to  cf.ns.id..  r  I  v.-o  ty;.' 
which  v.’o  label  a  and  M  for  dc  f  i  n  i  I  < ’in  ;;s  .  i.'  h  is  >: 

frequency  v;ith  the  initially  present  mcvle  (■,.  0),  v. 

is  offset  in  frequency  by  an  amount  /  0.  The  mcaii. 
satisfy 


'dt 


o !  ra  '  h -s  , 
a,  .,'1.  ■  1  1  !  .  .  , 

1  1  '■  :m  ■  ■  V 

.  iini  lilt  IS  !•  ■ . 


(A- 2..') 
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Associated  with  a^^  v/ill  be  the  pcrturbatioi-:  in  t!;o  ficlci, 
a^^  ( t )  ijg  (\ )  ,  which  we  denote  by  Uj(>:,t).  It  is  seen  iin;r.( -ciiate  1  y 
that  Uj  -■■..I  t  is  1' i  es 


TT 


(>^  (  t  )  i  (  t  )  >  J>I  (^  ) 


(A-24) 


Tlie  perturbation  V,Y,Q(t)  v;ill  qcnorally  oscillate  at  an 
acoustic  troquency  ,  so  Eq.  (A-24j  prodict.s  that  u,  will 
oscillate  at  this  frociuoncy  as  well.  Tiiis  mans  that^^ti'ie  per¬ 
turbation  in  the  actual  field  E  v/ill  oscillate  at  :.'q  t  ,  ^ ,  even 
though  the  cavity  mode  being  excited  would  os.cillate  naturally 
(i.e.,  in  the  absence  of  the  perturbation)  at  the  frequency  . 

(since  it  is  doqenoratc  in  frequency  with  t'le  initially  prcs.ent 
mode).  This  result  is  expected:  tlie  mode  N  is.  in  a  state  of 
forced  oscillation,  and  it  oscillate/:  at  the  ilrivu/n  fieqviency 

Wq  ± 

Equation  (7\-24)  is  in  the  form  of  the  .‘.t  u  tiiri  cquatifin  fr  r 
the  flux  perturbation  in  our  treatment  in  tlectiin  II.  !l  is  t  h\;  ■ 
seen  that  this  treatment  pertains  lo  th.e  flux  d.  i  s  t  urlwince  a  s.sot;  i.  a  t  ed 
with  the  buildup  c-f  a  cavity  mode  whos.e  urn  r's.  f;r  e.iency  is. 

degenerate  with  that  of  th.e  initial  mode.  si.'  w  in  bL/ction  IJ 

that  this  buildup  does  indeed  load  to  an  in;; I  abi  1  i  ty . 

On  the  other  hand,  we  now  consider  (lie  Iniildup  of  a  mode 
of  typo  M,  satisfying  Eq .  (.h-23).  .qince  osoillates,  at 

frequency  Eq .  (A-23)  predicts  tliat  the"Luildup  of  a.,,  will 

be  resonantly  cnlianced  wlK'n 


We  show  in  Sectie:i  IV  that  t).-'  i.uilaup  '.t  a  ted.-  c.  f  this 
type  can  also  lead  to  an  ins.tabi  1  ity ,  at  l<'as.t  in  the  1  ow- f  1  lumee 
limit. 


mode  f  IS  I 1  ari.  rs-t  W' 
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realist.  1C  las-'r  cav' ity.  .c 
or  nearly  rh'.q 'n' a  .i  t  <■  m'..!.;-, 
hand,  calcula.t  inn  t!'.'. 
pcctcd  to  ]jc  an  c,i:;y  t  is.k. 
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APPENDIX  B 


PROPERTIES  OF 
FOR 


degeni:rate  godixs  /TN'd 

IT  i  i  :  I R  RE  I N  l  0  RC 1 ;  M  E  N  T 


corjDiTior;s 


We  wish  to  derive  several  results  staLt.d  v.’ithuut  proof  in 
Section  II. 


First,  we  calculate  the  conditions  under  v.-hic)i  two  resonator 
modes  may  be  degenerate  in  their  ciejenf  rcquencios .  Consider  foi' 
simplicity  light  modes  contained  in  a  box  of  length  L  and  widtli  'f. 
The  eigenf requrcncies  (measured  in  cm"^)  satisfy 


(B-l  ) 


where  m  and  n  are  integers.  For  v  =  m/2L  n/2W  v.’c  may  write 
this  as 


V 

mn 


V  +  — 

mo  2  V 


mo 


(B-2) 


In  order  for  v  ,  _  to  equal  v  n  must  satisfy 
m-1 , 0  ^  mn  ^ 


(B-3) 


(B--1) 


where  1  =  v  ^  is  the  wavelength  of  the  light, 
mo 

Eejuation  P-l  suggerdr,  .nn  int'-e.'r  r.  ori  th,.'  orhi  i  ef  Id  v;  i  1  1 
permit  the  frequenc'/  of  a  tx’anrv-oT.e  i-.u-i..  t_o  F"  a  ■  's  ■  ;m  t  ( 

with  that  of  .i  longitudinal  jnode.  For  e.wir.ijd  e ,  let  W  a  em, 

L  =  400  cm,  and  A  =  10"^  cm.  We  obtain 


n  A-  1 G  , 


(B-'d 


Avea  ! 


From  this  vaJuo  wo  may  also  dotormiro  lln.'  na.'iibcr  n'  of 
transvorso  oscillatioiis  v;hich  are  to  bo  o:-;i><-o I od  in  an  aoour.t  ic 
wave  of  wave  vector  which  satisfies  the  ‘rlwase-iaa  tch  i  nq  con¬ 
dition  ‘ 


^  1  r, 

m-1 , 0 


+ 


•  ► 
k 

m ,  n 


(B-6) 


For  an  acoustic  mode  contained  in  the  transverse  dimension  of 
the  idealized  resonator  of  this  appendix  we  tlien  have 


=  Aa  =  W 

so  that 


(B-7) 


n'  ^  n/2  ^  8.  (B-8) 

This  is,  of  course,  only  a  ejua  I  i  fat  ivo  result,  mee'mt  to 
indicate  the  approximate  number  of  fringes  tc;  be  expef.-t  od .  If 
the  approximate  co:iditions  of  the  model  (uniform  steady  state 
properties  for  flux  density  and  gain)  wore  ritjorously  true,  n' 
would  liave  to  be;  a  integer  to  match  t!u>  acous-tic  bou'xla. ry  condi¬ 
tions  c'l  t  tlie  walls.  In  the  actual  case  Eq .  (B-8)  is  only 

approx  imiate . 

Next  we  show  that  the  wave  scattered  by  thc>  idiaso  grating 
(standing  acoustic  wave)  (.u'oducos  a  flux  perturbation  of  the 
same  shape  as  the  acoustic  wave. 

The  theory  of  scattering  of  light  by  an  ac'oustic  wave 
(Debye-Scars)  predicts  scattered  or  diffracted  waves  at  angles 
of 


(B-9) 


with  respect  to  the  initial  wave.  Tliesc  waves.,  .it  leas.l  '"or  a 
confocal  unstable  resonator,  will  after  one  round  Lrij',  return 
at  the  :  XU  lie-  aillUi’ii' 


TliP  in  ter  f  C'roJice  between  the;;e  waves  and  the  unpc’i' t  u!  l,ed. 
waves  will  be  of  the  form 


(B-lu) 


if  U]  lin.'’.  the  s.iiTK'  un;  ■  r  t  u  id '■ 'd  f  r  <  •.  j  n.  •!;>  •  y  (i  .e.  ,  if  I.)  1.)  ,\s. 

tile  original  wave.  Th"  s.ituation  is  i  1  1  u;.  1 1  a  t  .  in  I'luuie  p.- 1  . 
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cjuro  B- 1  .  1 1 1  Vist  ra  t  i  on  of  F''j>ro.!no  i  hi  1  i  ty  of  tlu;  Rc'attOT'od 

Wav'’o.  1'he  un;;o  t  t  urho;!  ■.••■.i'.v  focusor;  (\’i  r  tua  1  1  y )  at 
the  point  O  in  t  ho  oontmon  fcical  piano  of  t  fio  t  v.'o 
mirrorn.  'I’fio  Kca '  t  oinai  ’.vavo;.  a  1  :;o  fooun  i  thir. 
piano  hut  lit  si  i()tir  ly  d  i  .'-;p  1  a.tnnl  [Kunt/,  A  and  A'  . 

The  rf'’(.i!rri  'vaaor.  aft<'r  ono  cfnand  trip  '.■.■'11  t  n,  h.- 
in  lino  s.iiao  ion  as  t  ho  ori^linal  s.o.i  1  t  o  r  •  d. 
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Finally,  wo  justify  our  asrortion  that  tl-.o  donsity  pertur¬ 
bation  is  independent  of  z,  tlio  liyht  projiayatioii  direction.  If 
in  the  notation  of  Section  II  uj  ,  and  hene*'  ,  ha;;,  in  fact,  a 

different  k  ,|  than  fo/  then,  after  the  coninion  factor  exp  i  (kz  -  .ot 
is  factored  out,  j'lj,  can  be  written 

'i'u  ~  cos[k^x  -  (B-ll) 

Now  the  gain  and  hoatinc)  equations  v/ill  be  iiiconsir.  t  ent  unlo!'.;: 
gil  and  pjp,  have  this  sane  additional  z  behavior.  Hf)v.'ovc'r ,  nine- 
the  thickness  of  the  medium  is  of  the  order  of  half  tion  setruat'.i  : 
of  the  mirrors  (or  less)  in  a  typical  case,  the  z-vciriation  v.-p,  1 
be  in.suf f iciont  to  cancel  out  the  effoct.s  of  the  mode-medium  i;'.t,  . 
action  and  we  are  justified  in  neglecting  it. 
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In  this  ai^pondix ,  v;c  present  a  semi-h.ourii^t:  c  justifica¬ 
tion  for  Eq .  (GO)  of  ;  ;;c  text. 

Let  us  assume  th''  I'lcatinci  i  y>  onf^  d  i  mti'.s  i  ona  1  (flev:  direc¬ 
tion)  and  uniform  in  tn.e  reuion  (-  b/2)  ■"  (b/2).  In  thj  s 

case,  Fourier  and  Laplace  ti'ansf  orniiuq  Eq .  (57)  Iccids  to 


(k ,  s) 


-  2  0  -  1)  (s)  k  sin  k  ‘I 

(s  +  ivk)  [  (s  1-  ivk)^  +  c^^k‘ 


(C-1  ) 


where  F2(s)  represents  flic  tonii>ora]  h.isto:-y  of  tlie  lieatina:  and 
wc  have  us.ed  tiic  F- transf  orr.i  of  a  stes,.  1  u;u- 1  i  cn; .  Tr,en  ,  by  ccni- 
tour  integration  v;c  obtain 


P,(x=.0,s)  =  e 

C  S 


sj2  _  .sl5  J,  _ 

2v  2c. ' I  "’2 


sb  (  v/c 


cosh 


(C-2) 


Expandinq  al)f>ut  s  -  0  wo  find  that  f:.'-'  f  i  rr,  t  two  teriiu;  in  the 
brackets  arc  the  same  as  the  expansion  of 


(x  =  0 , s] 


(  y  -  1)  F(r.) _ 

2  s  +  ' 


(C-3) 


where 


1  _  _  V 

c  2 

s  c 


(C-a) 


li.0 
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\\’c  mny  arrive  at  a  result  similar  to  L’q.  (C-3)  by  m-akur; 
substitution 


s  +  ivk  '>•  s  + 


For  v/Cg  ~  0.5,  wo  have 


s  +  v_  SiS  +  1.25  p 
2  b 


In  t!ic  text  wo  have  taken,  somewhat  arbitrary  s  +  v/b  as 
dampinq  term  on  intuitive  qrounds  for  all  value's  of 
we  c'xpect  the  co.’ivcctivo  dampinq  to  bo  of  th.c;  order  of  (f 
time) . 
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RESONATOl^S 


1.  EMPTY  RESONATOR  (STABl.l-,) 

TJiero  aro  many  v;ays  in  wliich  tP-c'  oquat  icnn  sati;-!' i.od  by 
the  eigcnf unctions  of  an  eniuty,  stab]e  i’e5'-onci tor  can  be  formu¬ 
lated.  In  this  appendix  a  derivation  based  on  cbar^sjcal  poten¬ 
tial  theory  will  be  presented,  since  that  provides  t.he  most  lucid 
extension  to  the  analysis  to  be  presented  later  on. 

The  steady-state  wave  equation  for  a  tran.svorsc  electro¬ 
magnetic  wave  (Helmholtz  equation)  is  given  by  +  k^E  =  0 

where  k2  =  m  is  the  wave  f  rcqucr.cy ,  c  is  the  speed  of 

light  and  E  is  the;  Laplace-transformed  transverse  cc:'.r.,cn‘:nt  of 
the  electric  field.  The  solution  to  this  equation  can  be  writtcai 
in  terms  of  the  Huygons-Frcsnol-Kirchhof f  integral  by  moans  of 
Green's  identities ( 1 ) 


E (x ,y , z) 


(D-1) 


where  G  is  the  appropriate  Green's  function,  A  is  the  area  of 
interest.  Eg  and  r7E,_y?'n  are  the  electric  field  and  its  normal 
dorivativ'o  on  the  surface  respectively.  Thus,  once  the  field 
and  its  normal  derivatives  on  the  bounding  surfaces  are  given, 
the  field  cvcryv.'hcre  can  be  determined.  Tliis  rijjproach  h.as  tlio 
disadvantage  chat  it  requires  the  spcci  f  j  cat  i  on  of  both  tin.' 
electric  field  and  its  normal  derivative  on  the  surfactm  In 
general,  either  the  function  or  its  noimal  dcrivatit’c  will  be 
specified.  Those  are  the  Diriciilct  and  Neumann  boundary  con¬ 
ditions,  respectively.  Note  that  inixc'd  conditions  are  a  1  1  fneed 
wlicu'c  Dirichlet  is  specified  ever  part  of  th.e  boundiiw;  siirfacc 
and  Neumann  over  the  rest  of  it  but  not.  l.)Oth  at  tlu':  sar.ie  time 
over  the;  s.ame  .surface'.  Thus,  it  is  in 'cer-.sa  ry  to  ro  f  orr.vj  ]  a  t  e 
the  solution  in  terms  of  so  called  dipole  and  s.  i  i-ej  1 C' t  charge' 
distribution.  For  the  Dirichle.t  problem  tliis  implie.s  that(--) 


(1)  Dorn,  M.  and  I'.'ol  f ,  E.,  Principles  of  Optics,  Per  g.imon 
Pre.ss,  NY  (  1959). 

(2)  Miklilin,  S.G.,  Inteqral  Equations,  Pernamon  Pres,;-.,  NY 
(19:57). 
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I" 


t;(x',y,z)  =  —  G(x,x')  dA 


(0-2) 


and  the  dipole  layer  ;i  satisfies 


E  —  n  + 


i-  dA 


(D-3) 


v;hcrG  E  is  the  value  of  the  electric  fie]d  on  tlic  boundary,  n 
is  the  normal  direction,  A  is  the  surfeace  aj'oa  iind  G  is  the  ap¬ 
propriate  Green's  function  for  the  problem  at  hand.  For  tliree- 
dimensional  problems  it  is 


G  (x  ,x  '  )  =  e 


ik  •  r 


(D-4  ) 


Consider  applying  this  procedure  to  the  stnl:>'!e  resonator 
sho\'m  in  Figure  D-1  vdiich  consists  of  the  tv;o  mirrors  v;it,h.  un¬ 
equal  radii  of  curvature  separated  by  a  distance  D.  In  tl'.o 
Fox-Li  calcualtions  for  the  sLablc  resonator  only  th.c  fiedds  on 
t))c  mirrors  are  determined,  that  is  only  n  is  solved  for.  Since 
fields  on  the  mirrors  vanish  Eq.  (D-3)  reduces  to 


‘l  J  '^1  3n 


2  3n . 


(D-5) 


^2  + 


f  ^‘1 


2  'in. 


(D-G) 


where  uj  ,  \i2r  ,  M2  cire  the  dipole  layers  and  minor 
mirrors  1,  2,  res[)cct i  ve  1  y ;  ^11'  ^22 

tween  the  mirrors  including  c.-urva.turc  t'arm.s  and  ih- 
of  different  points  on  tlie  s.ime  mirre^r.  b'ritiivi  Lh.o 
terms  in  IF].  (D-C)  a:;  ;i  ;i  i  and  the  first  and.  ]a:;t  t 
(D-6)  as  Y2  V2  these  ('Cjuations  can  bt;  rewritten  a;: 


a  rear 
s  t  a  n  cs ; 


of 


'*2  i,T. 


(n-7) 


/C.3 
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Figure  D-1.  Empty  Stable  Resonator 
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Substituting  Eq.  (0-8)  into  (D~7)  results  in  the  cosiventiona  1 
integral  equations  for  the  Fox-Li  prolilcm.O)  In  order  to  solve 
these  equations  the  paraxial  approximation  is  made,  th.at  i;^,  the 
term  proportional  to  S~2  in  Eqs.  (0-7)  and  (D-3)  is  neglected. 

The  question  of  what  happens  when  the  resonator  is  loaded  will 
now  be  addressed. 

2.  LOADED  RESONATOR  (STABLE) 

The  Maxv.'ell  equations  appropriate  to  the  loaded  resonator 

are , 

V  •  C  E  =  -  7  •  P  (D-9) 

o 


V  •  B  =  0  (D-10) 


V  X  E 


M} 

'dt 


(D-11) 


V  X  H 


9E 

c  + 

o  9t 


3P 

at 


(D-12) 


where  it  has  been  assumed  tliat  there  are  no  freo-charge  currents 
flowing,  but  the  median  responds  to  electromagnetic  waves  through 
the  polarization  term. 

Equations  (D-9)  and  (D-IU)  can  bo  combined  into  a  single 
wave  equation, 

,  2  2 

E  -  -irr  — T  E  =  y  P - ^  V  V  •  P  (D-13) 

c  ■  9t  9t  o 


For  transverse  electromagnetic  waves  the.'  la;.t  term  va.nislic'S. 
The  polarization  term  is  assumed  to  be  given  by 


(3)  Fox,  A.G.,  Li,  T.,  Bell  Syst.  Tecli.  J  '1 0 ,  (1901) 
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I'  =  (a  +  c.c) 


(1)-]  ) 


wlicro  X  is  the  complex  dielectric  susceptibility.  It;;  form  in 
term-s  of  more  fundamental  quantitie;:  i;;  qivc'n  in  Aj  p  n  lix  11. 
Substituting  Eq.  {D-14)  into  the  Laplace:  ti  an;;loir.i  of  Eq .  (l)-l  i) 

results  in 


E  +  E  =  O  (D-15) 


where 

=  —  (^1  +  Xj-  +  iXi  +  i  (b-lC) 

where  Oj,  denotes  a  phenomenological  loss  term  duo  to  output 
coupling  from  the  resonator  (x^  <  0  for  lar.ii'.n  medivim)  .  The 
solution  to  Eq .  (D-16)  can  be  v/ritten  in  tiio  su.iik!  form  as 
Eq.  (D-1).  However,  since  again  only  tlie  field  or  t  lu,;'  normal 
derivative  would  normally  be  specified  it  it,  iiecctsarv  to  re¬ 
cast  the  problem  in  terms  of  oiie  for  the  singlet  or  cioublet 
charge  layers.  However  because  of  the;  angjl  i  f  ica  t  ion  Iv/  the 
medium  it  is  no\-i  necessary  to  calculate  the  fields  betwe-en  the 
mirrors  as  well.  Tlic  only  difference  in  tlie  douiilet  onuntic/n 
(Eq.  [D-3])  v^ould  be  that  the  Green's  function  has  a  crewing 
(decaying)  part  to  it  because  of  tlie  complex  values  k. 

Hence  special  attention  has  to  be  paid  to  the  converc.,nce  {uop- 
ertics  of  the  various  functions  if  equations  of  tlie  form  r;f  l.'qs. 
(D-7)  and  (D-8)  arc  to  be  derived.  Mo.'.t  of  the  di  f  f  i  cvi  1  t  ie  s  ar,,' 
swept  into  finding  the  appropriate  Grec>n '  s  function. 

At  this  point  an  alternative  procedure  could  b--  used. 

This  is  the  method  of  multiple  scalin<).  Thi;;  tecl.n ;  ljis'  i;; 
used  to  essentially  eliminate  thc'  rajud  variation  in  t  lie 
direction  (direction  normal  to  the  mirrors)  !  i a  in  th;'  gqu  ttions. 

To  implement  this  the  fields  are  as.sumed  to  liave  e*''^'  \Mria- 
tions  in  tlic  d  i  root  ion  .  Substitutin':  ihis  fca  in  i  n  t 'a  Ecj .  (D-lf) 
and  averaging  over  tlie  fast  variation  results,  in 

2 

E  +  2ik  — -  K  +  E  -  -  i  -  >.  E  -  ias  E  (I)-17) 

T  z  dz  z  .1  o  r 
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whore 


k 

z 


2 


2 


2 


o» 


-f 


(D-  ]  H) 


2 

Vr^,  is  the  transverse  Laplacinn  operator. 

TFie  solution  to  Eq .  {D-17)  can  in  {.irinciplc  Ijl'  1  ound  by 

the  following  procedure: 


E  (r) 


1-;  (  r  ■  ) 


c  (1  ,r  '  ) 


(D-iy) 


(for  a  homogeneous  medium)  and  G  sati.sfi 


G  +  2ik  — -  G  ^  k'  ^G 
z  o  z 


(r-r-  ) 


(I)-2  0  ) 


Note  that  E  and  G  have  to  bo  found  i,n  Ih.,'  .■•I'.tiro  r*  :  i  f'n  L" 'tva,  ■ 
the  mirrors,  not  just  on  the  mirrors  tin -'i; ,e  1  viu'.  as  Isui  bi.;on  tii" 
case  for  the  empty  resonator. 

Even  assuming  that  the  Green '  r;  function  c'ci  ie  found 
Eq .  (D-19)  is  stil]  an  integral  cjuatioii  v):ich  his  b'  be 

solved,  oithc'r  numerically  or  in  a  pi-r  t  u:  :  s;  t  ;>  n  .■■■  n  :<  n  {'r  '■rvi 
appproximation )  .  Solving  the  coupled  ret.  ot  i::r,.  ('.-1'.')  and 

(D-20)  is  more  akin  to  a  scattcM'in'!  c:a  1  cu  1  a.  t  ;  on  i  r,  it  i  ■■ 
using  the  Green's  theorem  of  the  pireviour,  r'ct  ion. 

3.  LOADED  RESONATOR  (LUnISTABLE) 


The  unstable  resonator  is  evv-n  more  cer;  1  i e  ■  i.  Ennire 
D-2  illustrates  why  tliir.  is  the  cm.'O'.  T.he  v.nr.a  i;  j  .  (D-IV)  :1  jl1 

api'jlies.  However,  it  has  to  Ije  jjol'.n  d  i  a  ■  ;  '  ;  I'-moi:'.  : 

the  solution;;  imve  to  liv  matciied  ac.-i'or.r  tl:-  i  onn  :  .;  :■  1-5)  - 

gion  I  contcains  both  a  r  i  giit-rm:)  inn  an.;  a  1  e  ;  ; -tr  > ;  n  ...cre.  -or. 

II  contains  only  a  r  i  giit.-cjoi  n<j  wa\’e  witli  coin.  Iv-e-ion  lli  C'  n- 
tainr-J  a  ricil.t  coing  v.-.e/e  but  with  no  (,.  0).  in  t  t;e 


analy.sis  w’here 

the 

e(  !  U  It 

.  i  on 

<1  r  '  ‘ 

r ' o  1  ■  d;  in  , !  i  1  ;  1 

.  • '  :  '  a  1 '  ■  1 1 :  ■  , 

0  -■  0  Ijf'CM  ur.e 

fc;r  .solviiK.;  E| 

til.' 

lor.s 

i CM 

,  leul 

at  ed.  1  i  ,  •  i  1  1  c  . 

•1  '  pi  o 

.  (D- 

17)  1 

i  to 

::i-  ■  a  r  i  ■ :  ‘ 1  '  :m  i  a 

;  evl  K,  ;t  or, 

wave  and  a  ]t:f 

t  -  CIO  1 

nq  ;; 

'her  1  c' 

■al  w 

■  •  in  ia  :  ;  on  I  , 

M  '  .  i : ,  u::,- ' 

a  c;y  1  i  nd  i  i  ca  1 

v;a  v<* 

with 

d  1  f  1  . 

T'  ■!!  t 

I  1  op.  1  '■  1 1  ion  cr  : 

1 :  ■  t  .  1 ; ,  t  : ;  in 

Reqii'i'.'s  11  .in'! 

Til, 

.ilicl 

t  h‘  -n 

m.it  o- 

ii  1  r. ■:  t  :.e  v;  ive  i.i 

'lot  ion;;  .  ;e  l 

the  d  i  r  t  (  r  ell  t 

bound 

1  U  1  *  ■ : 

; .  Tb 

lur  , 

the  field  in  Re ' 

I  1  on  1  i  'ii 

i)’/  ( 1  oi  t  tie  1  o 

(  > !  <  i ‘  1 

■  mod' 

) 

(4)  Mo.te,  G. 

'I' .  and  Met 

’a  r  thy 

',  H. 

,I .  ,  d  .  ('pt  .  .hoc  . 

AM,  G7,  22 

(1977)  . 


E 


ik  (z  +  d  4-  D) 
f  ( r  ,  0  )  e  ” 


G  ( ,  i- ) 


v/here  f  and  G  arc  to  be  detorniined .  The  boundary  cojidition  E  =  O 
on  the  two  mirrors  lead  to 


f (o,r) 


G(d  +  11,1-/  ( [:>  4  d )  ) 
D  +  d 


(0-22) 


f {D,r) 


G(d,  r/d) 
d 


(0-23) 


where  the  magnification  is  given  by  M  (D  +  d) /d  and  tlv’  lonui- 
tudinal  mode  number,  k  =  nr/0.  Substituting  Kq.  (0-21)  into 
(0-17)  leads  to 


It 

<)  Z 


2 

_T_ 

2ik 


(0-24) 


V 

^G  ^i;_ 
?p  2ik 


(0-2 b) 


where  V'p2  the  transverse'  I.aplacian  in  cry  1  i  nd  r  i  c  a  1  and  sidrei'ical 
coordinates  in  the  tv;o  ecjuaticjns  anc:  e  -  -  (  ,  2 /i;  ,,2c2 )  i;  ,  ,  ^  is  tiio 
ejain.  Equations  (0-24),  and  (0-23)  werc'  solvc'd  liv  .“'iaore  a.nd 
McCarthy(4)  usinc;  an  as;;amed  eain  profile.  Tlio  so  1  f-con;;i  r.tC'Ut 
gain  profile  is  cjiven  by  tlic  Higrod  eejuation. 


V 


+  P 


|0-2(.) 


wlierc'  V  is  tlio  flow  velocrity  of  tin-  la.-dium  ;  ,  is  the  dt'.:o1  ivat  ’(n 
time,  i  ,  are  the  optical  and  sat’i.ral  ion  {  1  u;-:  and.  !'  i:  t  !u  pn,.,p 

incj  rate.  Since  i  ~  id'  this  implies  that  the  (ijU'l  ren;-.  i  or  t  h*  ■ 

fields  ((O-2'l),  (0-25)1  are  coupli-d  thtouih  the  ciin  lie.  (O-.di). 

Since  thi'  two  field;;  are  drea' r  ibed  in  d  i  !  1  .  ■  i  i.  :  1 1  c  'a  i;,, 
terns,  this  m.i);;  ;;  the  .se  1  f-cons  i  :.t  <  m  ;e-,ii  a.,;-!  p;M  icalai  ly  cen- 
filc'X.  A  comiiU'f.e  numerical  analysis,  of  tlu  se  hqs.  (O-.dl)  ,  (O-.’h) 

and  (0-2G)  is  described  in  S-'Ct  ion  V. 
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In  the  integral  Cornvalation  ot  the'  inodo-nied  i  u;n  inataijilit 
the  solution  to  the  equilibrium  problem  i.s  written  aa;  u^  ^  n 

where  is  the  appropriate  iluygens-rresne  1  operator. 
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APPEND  IX  i: 


DERIVATION  OF  THE  POLAR  I  ZAT  10r«'  TEI'.M 


In  the  ni.in[)ln  eipi^roach  taken  In.M'c  thi'  <■  1  ec t rc;n f,  aja,'  treat 
nemic  1  ii  ns  i  ca  1  ly  .  Fjinco  the  median  treated  is  a  ]  c/.v-cien  r;  i  ty  a  a.  ;i 
the  electric  field  -v'etMi  by  the  i  nd  i  \’ i  diua  1  r-.'.o  1  ■ 'C  n  1  es  is  the  Fanv 
as  the  e:-:tc!'nal  field.  (In  a  d(‘n:;er  median  lik.a  a  solid  ti'c  e.‘ 
tive  f  i e  1  l!  would  be  Eeff  --  E  +  Il'’/3  wliere  IP  i ;;  tlie  ;  ol  an  ;uit  j  or. 
shieldincj  .  )  ( 1 ) 

Thas  individual  electrons  c.ibey  ti;e  following  'C  j'aa t  i ''m , 

X  .  .X  .  _  e 


where  y  ,  '.<3  are  the  danijiin'i  (crewth)  r.i'e,  fir'  eri.'it  il  l.-iiriin  i 
freejuency  and  E  is  tiu'  s  ti  iiui  1  a  t  :  na  ebetne  fnld..  IRdvina  Ej, 
(E-1)  leads  to 

-  i  .  t 

E,  e 

x=  <'-• 

L  ‘  o  B  ^  '  o  -* 

where  c-c  is  the  complex  conjuqate.  Tlte  dipole  moment  of  the 
molecule  is  given  by 

p  =  -ex  (E- 

Finally  the  macroscopic.-  pol.ir'  ation  for  a  density  of  N  moK-ctf 
per  cubic  centimeter  is  give. 

a'  =-  Nex 


IP  -N 


-  i  .m  t 
o 


m  2  i 

L  ‘  'o  ‘  h 


The  (--ffects  of  a  spatially  v.'nyin'’;  w.i'.’»>  ai"  inclsl 
fashion.  Th"  Fouriei-  t  r.ans  f  o  I'mc  al  p'O  !.■  n  xa  t  1  on  hir. 
nary  parts  givop  by, 


;c  !  u  h  I  :  n  ,1  .■  in 


1  1  .ini  11. 


1.  Pannfsky,  W.K.H.,  Phillips,  M.,  n.i-ml  ibc 
M,i  (jn  <  t.  i  an  ,  Add  1  s.i  in-V.'i  s.  1  e-y  ,  M.iss.  (i'lan). 


AVCO  i  .'I  HI 


where  R  =  “  '■'b2)2  +  Thus  the  real  part  moclific'S  tlie 

index  of  refraction  while  the  iraaq  inary  !>art  i^rodnces  rd.  ir.ird  atod 
emission  of  absorption  depending  up^on  the  population  inversion.. 

In  Section  IV  II>  is  written  as  follows; 

IP  =  f^oXj,E  +  ie^  E  (E-7) 


Comparing  Eq.  (E-7)  with  Eqs.  (E-5)  and  (E-6)  leads  to 


X'r 


Ne^ 

O' 

o 

Ne^ 

Y'.i 

( 

t  m 
o 

R 

(E-8) 


{E~9) 


Tn  addition  I!'’was  defined  in  terms  of  a.  Equationr'  (E-P)  ami 
(E-9)  give  a  simple  definition  of  the  frequency  dei'iendenco  of  i. 


In  the  somic  ]  assica  1  trcwitnent  crov.'t)i  (sfinul.ited  i.  i:ii  .s  .s  i  on  ) 
i.s  oljtaincd  by  letting  >  be  negative.  (2)  in  adflitir:'.,  the  effect: 
of  saturation  enter  y  by  makiim  it  aiin.ililardc  d.r  p.,'ne.' 'U  t  .  Thus  a 
relationship  between  th<}  junn  ic  V  a  s.-;  i  oa  1  ‘reatmenf.  ..nvi  th.e  curantura 
treatment  is  obtained  by  identifying  witli  tlu.'  c;ain  of  tiie 
medium.  The  gam  is  given  by 


9 


01 

o 


2 


c 


2 


N 

o 


.  k 
1  o 


(E-l 0) 


Hence 


2  0-  "  y:.- 

e  _ o  _  o 

nn  ,22  K  "  ‘o 
o  k  c 
o 


(K-ll ) 


finally  this  allows  the  pc  d  a  i' i  ra  t  i  on  t''  f'o  written  ;:i  teims.  ol 
the  gains  as, 


2,  Yariv,  A.,  (.>Ua:P  i ,  1 . -e  '  ;  e: . ;  e 


(1  ^'(.7)  . 
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(K-1 2) 


c  NE, 
o  R  k 


.  2  2 
k  c 

O  Cl 

o  2  ■  k 
Cl)  o 
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This  result  is  derived  on  the  basis  of  assuminci  that  and  fi 
remained  constant  over  the  time  period  of  interest.  Tiio  ueiierali- 
zation  allows  both  of  these  to  c.hancje  and  is  treatc'd.  in  Si'ctioM 
IV.  Note  that  even  in  cejuilibrium  rj  i  s  al  lo'wed  to  '.’ary  bi  CLi use 
of  saturation  dependence  on  the  fields.  Thiat  is  v  drynends  on  tin- 
fields,  but  is  assumed  to  liave  a  weak  dependence  on  the  fields 
at  Lest.  For  single  freejuency  electric  fields  the  polarizaticm 
in  real  space  and  time  is  given  by 

IP  (r  ,  t)  =  (H-13) 

where  IP)^  is  given  by  Eq.  (E-12).  For  m.ore  complicated  elect  ric 
fields  the  Fourior-Laplace  inversion  of  Eq .  (E-12)  has  to  be  per¬ 

formed.  This  leads  to  a  convolution  integral  of  tl'ic  form, 

IP  (r  ,t)  =  Jdr 'df  x(r',t')  E(r-r’,t-t’)  (E-14) 

where  X  is  given  by  Eqs.  (E-8)  and  (E-9). 
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